Impacts of entropy generation in second-grade fuzzy hybrid nanofluids on exponentially permeable stretching/shrinking surface

The present investigation aims to use entropy analysis to analyze the unsteady Magnetohydrodynamic (MHD) flow in a second-grade fuzzy hybrid \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left( {{\text{Al}}_{{2}} {\text{O}}_{{3}} - {\text{Cu/SA}}} \right)$$\end{document}Al2O3-Cu/SA nanofluid over an exponentially shrinking/stretching surface. The model for hybridization of the mixture of alumina \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left( {{\text{Al}}_{{2}} {\text{O}}_{{3}} } \right)$$\end{document}Al2O3 and copper (Cu) nanoparticles in the sodium alginate (SA) base fluid under heat source/sink, nonlinear thermal radiation, and viscous dissipation. The fundamental partial differential equations (PDEs) are simplified using an appropriate similarity conversion to generate the ordinary differential equations (ODEs). The analytical computation occurs in the MATHEMATICA program implementing the homotopy analysis method (HAM). In terms of code validity, our results are preferable to previous findings. The features of several parameters against the velocity, surface friction coefficient, entropy, temperature, and Nusselt number are described through graphs. According to our findings, the rise in the Brinkman and Reynolds numbers enhanced the total entropy of the system. Furthermore, the nanoparticle volume fraction and viscus dissipation magnifies the fluid temperature while retards the flow profile throughout the domain. Fluid velocity declined due to the Lorentz force using magnetic impact applications. The imprecision of nanofluid and hybrid nanofluid volume fractions was modelled as a triangular fuzzy number (TFN) [0%, 1%, 2%] for comparison. The double parametric approach was applied to deal with the fuzziness of the associated fuzzy parameters. The nonlinear ODEs convert into fuzzy differential equations (FDEs) and use HAM for the fuzzy solution. From our observation, the hybrid nanofluid displays the maximum heat transfer compared to nanofluids. This important contribution will support industrial growth, particularly in the processing and manufacturing sectors. The percentage increase in skin friction factor is 18.3 and 15.0 when \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha$$\end{document}α and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta$$\end{document}β take input in the ranges of 0 ≤ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha$$\end{document}α ≤ 0.8 and 0 ≤ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta$$\end{document}β ≤ 1, respectively.


Impacts
) and copper (Cu) nanoparticles in the sodium alginate (SA) base fluid under heat source/sink, nonlinear thermal radiation, and viscous dissipation.The fundamental partial differential equations (PDEs) are simplified using an appropriate similarity conversion to generate the ordinary differential equations (ODEs).The analytical computation occurs in the MATHEMATICA program implementing the homotopy analysis method (HAM).In terms of code validity, our results are preferable to previous findings.The features of several parameters against the velocity, surface friction coefficient, entropy, temperature, and Nusselt number are described through graphs.According to our findings, the rise in the Brinkman and Reynolds numbers enhanced the total entropy of the system.Furthermore, the nanoparticle volume fraction and viscus dissipation magnifies the fluid temperature while retards the flow profile throughout the domain.Fluid velocity declined due to the Lorentz force using magnetic impact applications.The imprecision of nanofluid and hybrid nanofluid volume fractions was modelled as a triangular fuzzy number (TFN) [0%, 1%, 2%] for comparison.The double parametric approach was applied to deal with the fuzziness of the associated fuzzy parameters.The nonlinear ODEs convert into fuzzy differential equations (FDEs) and use HAM for the fuzzy solution.From our observation, the hybrid nanofluid displays the maximum heat transfer compared to nanofluids.This important contribution will support industrial growth, particularly in the processing and manufacturing sectors.The percentage increase in skin friction factor is 18.3 and 15.0 when α and β take input in the ranges of 0 ≤ α ≤ 0.8 and 0 ≤ β ≤ 1, respectively.Non-Newtonian materials have attracted the attention of investigators for over a century because of their interdisciplinary characteristics and stimulating texture.Such fluids have been used in various industries, including manufacturing, metallic substances, chemical processing, materials for plastics, and foodstuffs.Glassware exhaling, biologic fluids, skin care products, synthetic fibers, healthcare, food, metallic twirling, shampoos and conditioners, and numerous other sectors use non-Newtonian fluids.These fluids possess a broad range of features and are categorized as dilatant, shear-thickening, shear-thinning fluids, or thixotropic.Rheologists have identified and studied various fluid simulations: Casson, Williamson, Maxwell, third grade, Burgers, micropolar, Oldroyd-B, Sutterby Cross, Jeffrey, and Sisko.But second-grade fluid behaves differently under different conditions, which explains the characteristics of shear-thickening/thinning and Newtonian effects.According to its dynamical properties, second-grade fluid is acclaimed and esteemed among researchers 1-4 .

Abbreviations
Stretching a plastic sheet, on the other hand, is not always linear.An exponentially stretched sheet's heat transport characteristics have a wider range of technical applicability.The outcomes of annealing and thinning wires made of copper are measured by the efficiency of heat transfer on the developing consistent surface, which grows exponentially with evolving temperature and velocity variation.During such methods, the kinematics of stretching and the concurrent cooling or heating considerably impact the eminence of the final products.One real-world application of an exponentially shrinking or stretching sheet is in materials engineering and manufacturing, particularly in the production of micro and nanoscale devices.By manipulating the dimensions of a sheet exponentially, it is possible to create structures with precise and controlled variations in size.Overall, the use of exponentially shrinking or stretching sheets provides a versatile approach for fabricating structures with precise and controlled variations in size.These techniques find applications in various fields, including microelectronics, optics, biotechnology, and advanced materials research.Khan and Sanjayanand 5 studied a second-grade fluid's steady flow and heat conductivity over increasing exponentially surfaces employing the Runge-Kutta fourth-order (RK4) technique.Rehman et al. 6 explored the stable flow of a second-grade fluid across an exponentially stretched sheet employing the Keller-box technique and the HAM.Nadeem et al. 7 examined the flow and transmission of the heat of second-grade (viscoelastic) fluid through an exponentially expanding sheet in an environment of thermal radiation.Ramzan and Bilal 8 studied time-varying MHD, thermal radiation, and mixed convection and conduction flow on the stretched surface in a second-grade nanofluid.Through perturbations assessment, Pakdemirli et al. 9 evaluated the features of second-grade fluid.Over the past few decades, there has been a lot of investigation on the second-grade nanofluid over an exponentially extending surface [10][11][12][13][14][15] .
Scientists and technicians appreciate a chaotic flow in several industries because it gives them more significant influence over their decisions 16 .Moreover, even in the most effective flow scenarios, an undesirable insecure impact can appear encircling a structure.The behaviour of unstable BL flow is different from that of regular flow because of the inclusion of additional time-dependent aspects of the mathematical models that influenced the boundary layer (BL) segregation and the movement of the fluid pattern.Still, modern engineering techniques that enable improved reliability, effectiveness, and economics for multiple dynamic components are feasible via an adequate assessment of the unstable flow of fluid implications for production processes 17 .Zaib et al. 18 described the computational investigation of a time-dependent flow with heat flux past an exponentially contracting surface.
Thermal radiation is a phenomenon that causes heat transmission through electromagnetic waves.It develops since there is significant variation in temperature across the two mediums.Radiative influences are significant in the engineering profession and sciences.In the polymer manufacturing sectors, where heat-controlling variables influence the ultimate product quality to some extent, thermal radiation impacts are essential in controlling heat transfer.Moreover, weapons of mass destruction, planes, combustion engines, solar radiations, spacecraft, metallic liquid fluids, MHD accelerators, and nuclear plants have significant radiative implications.Pantokratoras and Fang 19 were the inaugural researchers to look at the impacts of nonlinear radiation from nonlinear sources on the Sakiadis flow.Dogonchi and Ganji 20 tested the influence of radiant heat on the MHD flow over divergent/convergent channels.Khan et al. 21studied the electromagnetic radiation propagation of hybrid nanofluid technology throughout an opaque surface.Non-linear thermal radiation has captured the attention of multiple investigators [22][23][24][25] .
Recognizing the need for improved thermal conductivity in traditional fluids, a new type of nanofluid called "hybrid nanofluid" is presented to provide highly industrialized heat conductivity.Two or more semiconductor materials are mixed with a base fluid to make a hybrid nanofluid.Different nanomaterials include carbon nanotubes 26 , metals, metal oxides, and carbides.Numerous investigators are now interested in hybrid nanofluid due to its significance for the betterment of thermodynamic characteristics in real-world applications 27 , as a result of Choi and Eastman's 28 outstanding findings that gave the unique notion of nano liquid.Hybrid nanofluids are also used in various applications, including electrical gadget cooling 29 , cooling of domestic refrigerators 30 , automobile braking fluid, transformers, heat exchangers, and solar water heating 31 .Suresh et al. 32 explored the effects of a hybrid (Al 2 O + Cu/Water) nanofluid in a circular tube that was uniformly heated.Momin 33 investigated the thermal act of a hybrid nanofluid in a spherical tube and demonstrated that the hybrid (Al 2 O + Cu/Water) nanofluid improves thermal conductivity compared to a conventional working liquid.Waini et al. 34 scrutinized the influence of buoyancy on hybrid (Al 2 O + Cu/Water) nanofluid toward the stagnation point of an exponen- tially stretching/shrinking vertical sheet.They determined that (Al 2 O + Cu/Water) had a higher heat transfer rate than Cu/water.Khan 35 numerically examined the convection of copper (Cu + Water) and (Al 2 O/Water) nano-liquid across a spinning disc in a porous media.Cu-water has a faster heat transfer rate than Al2O3-water, and the presence of porous media raised the thickness of the thermal BL.Takabi and Salehi 36 analyzed the heat transfers of Al 2 O/Water and (Al 2 O + Cu/Water) nanofluids with a heat source.The literature is well-stocked with further information on this topic [37][38][39][40][41] .
To reduce wasted energy, the research community has made extraordinary efforts.They were going through energy-converting devices again, coming up with new strategies and approaches to use better what they had.The amount of energy that does not accomplish productive work is called entropy generation.The ratio of recovered heat to temperature is defined as entropy generation.They measured the performance of thermal types of machinery such as heat pumps, electronics, power plants, air conditioners, and energy gadgets employing entropy production.Bejan 42 utilized entropy optimization to show how thermal performance in fluids behaves.Butta et al. 43 conducted a theoretical study of the entropy implications of a Casson nanofluid flow across an unsteady stretching surface.They observed that unpredictability aided the creation of entropy effects.The entropy generation in an MHD viscoelastic (Second-grade) fluid past a linearly stretched sheet was investigated by Aiboud and Saouli 44 .Afridi and Qasim 45 investigated the entropy effects of a 3D Viscous liquid across an exponential sheet.Butta et al. 46 investigated the entropy generated within a mixed convective magnetohydrodynamic flow of a viscoelastic fluid over a stretching sheet.They observed that unpredictability aided the creation of entropy effects.Dalir et al. 47 demonstrated entropy generation by a viscoelastic nano-liquid against a stretchable layer.Sithole et al. 48use a stretchable surface with viscous dissipation to examine irreversibility in the electromagnetic flow of second-grade nanoparticles.Some relevant entropy optimization efforts are described in [49][50][51] .
Fuzzy set theory (FST) 52 has proved to be a valuable technique for modelling uncertainties in recent decades, providing models with a more accurate view of reality and allowing them to express themselves with a wider perspective 53,54 .After modelling real-world problems, they convert into PDEs or ODEs.Uncertainty issues may arise during the development of a dynamic model.Researchers must deal with inaccurate data, parameters, dynamical variability, and complex relationships.As a result, many scientists use fuzzy models to depict dynamical systems to prevent artificial data accuracy and produce more realistic results.The FDE is critical in overcoming these challenges.Chang and Zadeh 55 proposed the basic idea of fuzzy derivatives.Dubois and Prade 56 proposed the idea of fuzzy numbers (FNs) for solving an FDE.Kaleva 57 introduced the concept of FDEs in a fuzzy environment.Recently, FDEs played a significant part in the research field of fluid dynamics, such as when Nadeem et al. 58 systematically studied the influence of MHD and gravitation on the third-grade fluid across an inclined channel in an environment that is fuzzy.The triangular membership function (MF) was applied to analyze the degree of indeterminacy.Siddiqui et al. 59 explored the heat transmission of SWCNTs MWCNTS on third-grade nanofluid down an inclined channel in a fuzzy environment.They used the volume fraction of nanoparticles as TFN for comparison and uncertainty.The effect of fuzzy volume fraction on Jeffery-Hamel nanofluid across inclined plates was studied by Biswal et al. 60 .Some current studies relating to the FDEs are included in [61][62][63][64] .
On the other hand, a comprehensive examination of the previously cited literature exposes substantial discrepancies and drawbacks.No preceding works have investigated entropy generation due to unsteady MHD flow in second-grade hybrid nanofluid over the exponentially stretching/shrinking sheet with non-linear thermal radiation and fluid dissipation in their investigation structure, according to the authors' understanding.Furthermore, the nanoparticle volume percentages for each nanofluid and a hybrid nanofluid are determined as fuzzy triangular numbers utilizing the dual parametric concept for comparisons and instability.The homotopy analysis technique was used to tackle the problem under consideration.The impact of essential parameters on heat and flow properties and entropy generation are graphed and briefly discussed.This creative approach can help in the development of industrial production, particularly in the areas of engineering and fabrication.We will be able to answer the following research questions once we have solved this problem: • What causes a decline in the nanofluid and hybrid nanofluid velocity profile when using the magnetic effect accelerates Lorentz forces?• What are the differences in transport behaviours between non-Newtonian (second-grade fluid) and hybrid nanofluid?• How does the entropy generation behave concerning different flow parameters?• How do the fuzzy phenomena affect the temperature of nanofluids and hybrid nanofluids?
• Why does the rise in the magnetic field, non-linear thermal radiation and the nanoparticle volume percentage cause the nanofluid temperature to elevate? • How do the Nusselt number and drag force respond to different flow parameters?

Mathematical model
This study takes into account the time-dependent, 2D incompressible flow of MHD viscoelastic (Second-grade) hybrid Al 2 O 3 + Cu/SA nanofluid over the exponentially shrinking/stretching surface, as shown in Fig. 1. u w (x, t) = e x / L (a/(1 − ct)) signifies the stretching/shrinking velocity, where represents a constant that relates to stretching > 0 and shrinking < 0 cases of the velocity rate, L indicates the characteristic length, and c denotes the unsteadiness.v w (x, t) = ν o e x / 2L / √ 1 − ct represents the mass flux velocity, where v o is constant.The variation in temperature near the surface is determined by T w = T ∞ + e x / 2L (T o /(1 − ct)) , the ambient and reference temperatures, denoted as T ∞ and T w , respectively.Also, the fluid at the surface is thought to contain a wall temperature higher than the ambient level.The magnetic field is supposed to be uniform , with B o , denoting a magnetic field, and there is no pressure gradient on the surface.The vicious, source/sink, and nonlinear thermal radiation effects are also considered.
where α 1 , and α 2 are material constants, where p stands for pressure, Ŝ is Cauchy tensor, µ is the dynamic vis- cosity, I is an identity matrix, d dt is material time derivative T indicates the matrix transpose, A 1 , and A 2 are the first two Revilin-Erickson tensors 1-4 .
When using the BL approximation, the following assumptions are used in calculating the governing equations for continuity, momentum, and heat: 12,48 the boundary conditions are: (1) where u and v are the velocity components along the x-axis and y-axis, respectively, and T shows the fluid tem- perature.The sheet surface's suction and injection velocities are assumed to be distributed.The dimensionless wall mass transfer coefficient is denoted by v w here.It displays the mass injection process when we select a velocity of v w > 0 across the boundary wall's surface, and the mass suction process when we select a velocity of mass v w < 0. According to the Roseland approximation for radiation in Eq. ( 7), q r is the radiative heat flow is defined as follows: in which the k * * and δ * * are the delegate of the mean absorption coefficient and Stefan-Boltzman constant, respectively.
The thermophysical properties of the hybrid nanofluids are 63 : To clarify Eqs. ( 3)-( 6), the subsequent similarity transformations are presented 10 in a considerably simpler method, where the stream's function ω is adjustable as v = −∂ω ∂x, u = ∂ω ∂y, and its similarity variable is η : Considering Eq. ( 9), the subsequent nonlinear ODEs Eqs.(3)-( 6) are capable of being condensed in the framework of the abovementioned links: e a t g e n e r a t i o n / a b s o r p t i o n p a r a m e t e r i s , the suction parameter is s = −v o 2L aν f and the stretching/shrinking parameter is .The coefficient of skin friction C fx and the local Nusselt number (Nu x ) are thus demarcated as: Then use ( 9) into ( 13) and ( 14), which produces an apparent connection: where Re x = u e x ν f is the x-axis local Reynolds number.

Entropy analysis
Entropy generation is a concept used in thermodynamics and fluid dynamics to quantify the irreversibility and losses that occur within a system.It measures the inefficiency or degradation of energy within a process.Expert's concerns about energy management in any experiment or study were sufficient to prevent valuable energy from vanishing.Entropy formation analysis might be more efficient for carrying out this crucial activity.This allows doctors to discover issues such as MHD impacting its irreversibility.While implementing the boundary layer approximation, the equation for the local volumetric rate of entropy development in the presence of a magnetic field for the second-grade hybrid nanofluid is 48 : here viscous dissipation, heat transmission, and magnetic field are all crucial elements in the formation of entropy effects, as seen in the above formula.
Equation ( 9) can be used to construct the dimensionless form of Eq. ( 16)'s formula for local entropy generation: here, the dimensionless temperature gradient is A = T w −T ∞ T ∞ , the local entropy generation number is , and the Brinkman number is Br = Pr Ec.

Homotopy analysis method (HAM)
The HAM is a sophisticated analytical technique for solving nonlinear equations with many variables.The HAM technique is used to determine the consequent Eqs. ( 10) and (11) with the help of the boundary conditions (12).
To start the process of employing this technique, linear operators and preliminary predictions are essential.Consequently, we utilized the previously mentioned technique to resolve the energy and motion transformation models employing � f , � θ as linear operators and f 0 (η), θ 0 (η) the initial estimations.See 9-12 for more information on this technique.
where h f and h θ signify non-zero auxiliary parameters and q ∈ [0, 1] represents an embedding parameter while F and θ representing the mapping occupations for f and θ.
The boundary conditions become the Eqs.( 10) to ( 12) convert in the non-linear operators like Eqs. (22-24), then the series solution becomes:

For fuzzification
The governing FDE is converted into a double parametric form.In this case, Eq. ( 8) can be converted into interval form using the χ-cut concept.Here χ and ω are parameters that range from 0 to 1, controlling the fuzziness of the uncertain parameters.The named problem was solved using the HAM as well.See for further details 53 .TFNs come in especially handy when you have higher confidence in the center estimate while allowing for some variance and the range of probable values is poorly understood.They offer a straightforward and natural method of modelling and representing uncertainty in decision-making.Using fuzzy concepts, comparing nanofluid and hybrid nanofluid is also explored in this study.The nonlinear ODEs convert into an FDE, and the nanoparticle's volume percentage is a TFN.The slight variation in the volume percentage of nanoparticles (19) F 0; q = s, F′ 0; q = , θ 0; q = 1, impacts the flow rate and heat.These parameters alone determine the nanofluid's flow rate and heat transfer because some researchers estimate that the volume percentage of nanoparticles falls within the [1-4%] range.It is preferable to address a challenging situation in a fuzzy atmosphere by getting both volume fractions as a TFN, since φ 1 and φ 2 signify the volume fraction of Al 2 O 3 /SA and Cu/SA, separately as shown in Table 2.The volume fractions of nanoparticles used in this study are classified as TFNs, with the TFNs being transformed into χ-cut methods, and the fuzziness of the TFNs is controlled by χ-cut [61][62][63][64] .
Let φ 1 = φ 2 = [ 0, 0.05, 0.1] be a TFN that is described by the three values highlighted in Fig. 2: 0 (lower bound), 0.05 (most belief value), and 0.1 (upper bound).The membership function can be represented by the TFN as follows: The σ -cut technique is employed to convert TFNs into interval form and is represented as The FDEs are renewed into lower θ 1 (η, χ) and upper bounds.θ 2 (η, χ). to handle this scenario.

Results and discussion
An unsteady flow analysis was carried out on the second-grade hybrid nanofluid over an exponential surface.In this situation, a fluid that conducts electricity is studied because it combines viscoelastic and nonlinear radiative heat transmission.Examine the effects of dynamic factors on our system's velocity, temperature, and entropy generation contours.Moreover, an approximation analytical method has been used to control the framework's modified equations.We have acknowledged the most significant characteristics of the simulation we run as , Pr = 12 and Ec = 0.3.Table 3 was generated to confirm the numerical results of θ ′ (0) Haider et al. 51 .The numerical results generated by the HAM in this research have been determined to exist in outstanding agreement with previous research.
Figure 3a illustrates how the magnetic (M) parameter controls the velocity field.For greater values of M, the flow rate dropped in both cases.Lorentz pressure is responsible for this phenomenon, which arises from the cooperation of electric and magnetic fields during an electrically conducted fluid flow.So, the fluid velocity in the BL is controlled by the generated Lorentz force.As a result, as M rises, the velocity of the fluid and hybrid nanofluids falls.The interaction of magnetic fields is significant in different technical and industrial applications such as crude oil extraction, geothermal systems, groundwater hydrology, and so on.Figure 3b illustrates the variation of the second-grade parameter ( α ) on mobility.The technical meaning of it is the square of the shear stress divided by the first normal stress difference.It is found that α increasing generates a boost in the veloc- ity of both liquid and nanofluid hybrids.This is because viscous forces and fluids' density decrease when alpha 0, otherwise.
Table 2. φ 1 and φ 2 transforms into TFN [61][62][63][64] .www.nature.com/scientificreports/grows.Figures 4 and 5a,b demonstrate the impact of an unsteady parameter ( β ) and the suction parameter (s) on the flow and temperature fields.The respective temperature and profiles of velocity drop as β and s rise.This demonstrates that with higher levels of β , the rate of cooling is substantially faster.This is because rising both and s decreases momentum and thermal BL.This is consistent with the concept that suction reduces the thickness of the outermost layer of fluid.Figure 6a,b shows the impression of stretching/shrinking parameters ( ) on velocity Table 3.Comparison of current results of θ ′ (0) Haider et al. 64 for variation in the Pr and M when H = 0.0, Nr = 0.0, β = 0.0, Ec = 0.0.www.nature.com/scientificreports/and temperature distributions.When increased, the velocity also upsurges while the temperature diminishes.

Fuzzy numbers
Because the stretching parameters are set to higher levels, the temperature and thickness of the BL are reduced.Due to the exposure of the cooler to the ambient fluid, the BL thickness reduces with growing values of stretching parameters.Figure 7a,b shows the variation of nanoparticle volume fraction (φ 1 ) on the velocity and temperature distributions.When φ 1 has raised, the velocity declines while the temperature boosts up.The variability of volume fractional of nanoparticles (φ 2 ) on velocity and temperature gradients is exposed in Fig. 8a,b.When φ 2 progress, speed drops while temperature upsurges.The main reason for the decay in the velocity is that as the values of volume fractional of nanoparticle growths, the resistive force also increases, reducing the fluid flow speed.Physically, the energy is discharged from the exponential sheet due to the resistive force of the nanoparticles.More energy is generated when more nanoparticles are added, causing the temperature to rise.Furthermore, the optimum temperature may be achieved because hybrid nanofluids have higher thermal conductivity than nanofluids.The characteristics of the radiation parameter (Nr) and the liquid temperature are seen in Fig. 9a,b.Electromagnetic radiation produced by the heat-driven movement of particles within materials is known as thermal radiation.
Higher predictions of the Nr help in the unpredictable mobility of nanoparticles.Consequently, there have been greater collisions among particles and more energy released.As a consequence, the fluid heat rises.The thermal profile extends as the temperature ratio parameter (θ w ) grows, as shown in Fig. 9b.The findings illustrate that as θw increases, the temperature disparity (T w − T 0 ) increases and causes the fluid temperature to rise. Figure 10 pierces the heat generation parameter (H) impressions on the temperature field.The non-dimensional heat generation/absorption parameter is based on the heat produced or absorbed by the fluid.It is noticed that as the    www.nature.com/scientificreports/H > 0 grows, the temperature distribution improves.A physically higher thermal production value indicates that more heat is generated within the boundary layer, which produces higher temperatures in the field.Figure 11 demonstrates the upshot of the Reynolds number (Re) on the entropy rate (Ns).The ratio of viscous forces to inertial forces in a fluid that experiences relative internal movement due to varying fluid velocities is described as the Reynolds number.A rise in Re causes the entropy rate to rise.Physically, since greater values of Re lead to larger inertial force and decay of the viscous force, the entropy rate hikes near the exponentially stretching surface.Figure 12 shows the temperature difference parameter (A) 's contribution to the entropy rate near the exponentially stretching sheet.When A is raised, then the entropy rate also boosts.Figure 13 influences the effect of α on Ns.It is seen that Ns upsurge due to the rise in α.Physical, Entropy formation is encouraged by stronger elastic effects.Figure 14 outlines the impact of the Brinkman number (Br) on the entropy rate (Ns).It is the proportion of heat transferred by molecular conduction to heat generated by viscous dissipation.The entropy rate is elevated via Br.Because Br and thermal conduction are inversely proportional.A larger value of Br causes low thermal conduction, consequently enhancing the entropy rate near the surface.Also, the hybrid nanoparticles enhance the entropy rate compared to a regular fluid.Hence, the quantity of heat energy is carried out.
Figure 15 the impression of varying values of M on Ns is shown.It was shown that the entropy rate improves as a result of M. The higher value of M leads to an upsurge in the resistance to the movement of liquid in the exponentially stretching sheet.Consequently, heat boosts and the entropy rate also rises.
As demonstrated in Fig. 16, f′′(0) increases as M increases, whereas it declines as α increases.Rising M gener- ates substantial obstacles to fluid flow because of the Lorentz drag force, which decreases the velocity of the fluid and momentum BL width, raising velocity and, consequently, shear stress at the exponentially stretching sheet.The behavior of f ′′ (0) the unsteady parameter (β) and suction/injection parameter (s) is revealed in Fig. 17.It can     be detected the drag force declines with the rise in β and s. Figure 18 impact of Nr and H on Nu x .It is observed that the Nu x reductions with an increase in Nr and H.The Nu x lower has φ 1 and φ 2 increased, as shown in Fig. 19.Physically, the heat emitted from the exponential sheet when enhancing φ 1 and φ 2 ., is displayed on the vertical axis, while the fuzzy temperature profile for different values of η, is displayed on the horizontal axis.The derived fuzzy temperatures are TFN, nevertheless not symmetrical, while the derived fuzzy volume fraction contains both TFN and symmetrical.These variations might be due to the nonlinearity of the governing FDE.It was also revealed that hybrid nanofluids had a wider width than nanofluids.As a result, the hybrid nanofluid is uncertain according to the TFN.On the other hand, Fig. 20 shows the comparison of nanofluids Al 2 O 3 /SA (φ 1 ), Cu/SA (φ 2 ), and Al 2 O 3 + Cu/SA hybrid nanofluids through MF for numerous values of η.In these figures, we evaluated three scenarios.Blue lines represent the situation in which φ 1 is considered as TFN and φ 2 = 0 .Whenever the φ 2 is considered TFN and φ 1 = 0 , it is represented by red lines.In the third case, black lines show the hybrid nanofluid that has both φ 1 and φ 2 non-zero.It has been determined that the hybrid nanofluid operates well since its thermal variation is greater than that of the other nanofluids.Physically, the combined thermal conductivities of Al

Conclusions
This study analyzed the unsteady MHD flow caused by exponentially stretching/shrinking surfaces in a secondgrade hybrid (Al 2 O 3 + Cu/SA) nanofluid.Alumina (Al 2 O 3 ) and copper (Cu) nanoparticles are mixed with sodium alginate to create a hybrid (Al 2 O 3 + Cu/SA) nanofluid.Entropy generation, viscous dissipation, heat scores/sink, and non-linear thermal radiation are also considered.An analytical approach such as HAM is implemented for the outcome of the suggested problem.For validity, present results were found to be good when comparing existing results.The impacts of non-dimensional regulating parameters on velocity and temperature profiles for second-grade fluid and hybrid nanofluid are examined and discussed via graphs.Further, φ 1 and φ 2 are said to be TFNs using the χ-cut technique.Comparison and uncertainty are studied through triangular fuzzy graphs.The foremost goals of this study are presented below: • Due to strong Lorentz effects, the fluid velocity is dropped with the magnetic parameter, which is boosted with the second-grade fluid parameter.However, this retardation to the flow of hybrid nanofluid is higher than the retardation experienced by the flow of second-grade fluid.• The fluid temperature is boosted while the velocity profile declines with the improvement of φ 1 and φ 2 .
• The fluid temperature boosts against higher θ w , Nr, and H, values, whereas the reverse holds for the unsteady parameter, suction parameter, and Prandtl number.However, the radiations emitted by the hybrid nanofluid are more potent than those emitted by a second-grade fluid.• The fluid velocity grows with an upsurge in the stretching/shirking parameter, whereas the fluid temperature declines.• The entropy profiles enhance as the second-grade fluid, Brinkman number, magnetic parameter, and Reynolds number rise.• Heat transmission entropy effects dominate the exponentially stretching/shrinking sheet, whereas viscous dissipation and magnetic field entropy impacts dominate heat transmission entropy impacts in the distant region.Finally, we conclude that the entropy effects of heat transmission are the primary source of entropy production.The hybrid nanofluid shows maximum heat transfer in entropy production compared to the second-grade fluid.• Skin friction coefficients are reduced with increased unsteady and second-grade parameters while growing with magnetic parameters.• For higher values of Nr, H, φ 1 and φ 2 , the rate of heat transfer at the surface decreases.
• The variation among the lower and upper limits of fuzzy temperature profiles of hybrid nanofluids is highest for fuzzy analysis employing the triangle membership function.Hence, the degree of fuzziness is highest when compared to conventional nanofluids.
of entropy generation in second-grade fuzzy hybrid nanofluids on exponentially permeable stretching/shrinking surface Rana Muhammad Zulqarnain 1 , Muhammad Nadeem 2 , Imran Siddique 3 , Adeel Mansha 4* , Abdullatif Saleh Ghallab 5* & Mahvish Samar 2 The present investigation aims to use entropy analysis to analyze the unsteady Magnetohydrodynamic (MHD) flow in a second-grade fuzzy hybrid (Al 2 O 3 − Cu/SA) nanofluid over an exponentially shrinking/stretching surface.The model for hybridization of the mixture of alumina (Al 2 O 3

Figure 11 .
Figure 11.Impression of Ns for Re.

Figure 12 .
Figure 12.Impression of Ns for A.

Figure 14 .
Figure 14.Impression of Ns for Br.

Figure 20
Figure20portrays the calculated fuzzy temperature using volume fractions of φ 1 and φ 2 as the TFN [0%, 10%, 20%].Four subdivisions for triangular MFs demonstrate the fuzzy temperatures for various thresholds of η, the following: 1, 2, 3, and 4. The MF of the fuzzy temperature profile ζ -cut (0 ≤ ζ -cut ≤ 1), is displayed on the vertical axis, while the fuzzy temperature profile for different values of η, is displayed on the horizontal axis.The derived fuzzy temperatures are TFN, nevertheless not symmetrical, while the derived fuzzy volume fraction contains both TFN and symmetrical.These variations might be due to the nonlinearity of the governing FDE.It was also revealed that hybrid nanofluids had a wider width than nanofluids.As a result, the hybrid nanofluid is uncertain according to the TFN.On the other hand, Fig.20shows the comparison of nanofluids Al 2 O 3 /SA (φ 1 ), Cu/SA (φ 2 ), and Al 2 O 3 + Cu/SA hybrid nanofluids through MF for numerous values of η.In these figures, we evaluated three scenarios.Blue lines represent the situation in which φ 1 is considered as TFN and φ 2 = 0 .Whenever the φ 2 is considered TFN and φ 1 = 0 , it is represented by red lines.In the third case, black lines show the hybrid nanofluid that has both φ 1 and φ 2 non-zero.It has been determined that the hybrid nanofluid operates well since its thermal variation is greater than that of the other nanofluids.Physically, the combined thermal conductivities of Al 2 O 3 and Cu are incorporated in a nanofluid with a hybrid structure to offer maximal heat transmission.When discussing Al 2 O 3 /SA and Cu/SA nanofluids have Al 2 O 3 /SA an enhanced heat transfer rate while Al 2 O 3 having higher temperature conductivity than Cu.
Figure20portrays the calculated fuzzy temperature using volume fractions of φ 1 and φ 2 as the TFN [0%, 10%, 20%].Four subdivisions for triangular MFs demonstrate the fuzzy temperatures for various thresholds of η, the following: 1, 2, 3, and 4. The MF of the fuzzy temperature profile ζ -cut (0 ≤ ζ -cut ≤ 1), is displayed on the vertical axis, while the fuzzy temperature profile for different values of η, is displayed on the horizontal axis.The derived fuzzy temperatures are TFN, nevertheless not symmetrical, while the derived fuzzy volume fraction contains both TFN and symmetrical.These variations might be due to the nonlinearity of the governing FDE.It was also revealed that hybrid nanofluids had a wider width than nanofluids.As a result, the hybrid nanofluid is uncertain according to the TFN.On the other hand, Fig.20shows the comparison of nanofluids Al 2 O 3 /SA (φ 1 ), Cu/SA (φ 2 ), and Al 2 O 3 + Cu/SA hybrid nanofluids through MF for numerous values of η.In these figures, we evaluated three scenarios.Blue lines represent the situation in which φ 1 is considered as TFN and φ 2 = 0 .Whenever the φ 2 is considered TFN and φ 1 = 0 , it is represented by red lines.In the third case, black lines show the hybrid nanofluid that has both φ 1 and φ 2 non-zero.It has been determined that the hybrid nanofluid operates well since its thermal variation is greater than that of the other nanofluids.Physically, the combined thermal conductivities of Al 2 O 3 and Cu are incorporated in a nanofluid with a hybrid structure to offer maximal heat transmission.When discussing Al 2 O 3 /SA and Cu/SA nanofluids have Al 2 O 3 /SA an enhanced heat transfer rate while Al 2 O 3 having higher temperature conductivity than Cu.

Figure 15 .
Figure 15.Impression of Ns for M.

Figure 16 .
Figure 16.Impression of M and α on C f .

Figure 17 .
Figure 17.Impression of s and β on C f .

Figure 18 .
Figure 18.Impression of Nr and H on Nu.

•
Based on triangle fuzzy MF, Al 2 O 3 + Cu/SA nanofluids that are hybrids are particularly able to enhance the temperature flow rate compared to Al 2 O 3 /SA and Cu/SA nanofluids in fuzzy heat transfer examinations.Moreover, when examined, the Cu/SA nanofluid outperformed the nanofluid Al 2 O 3 /SA.
The dynamic viscosity of Al 2 O 3 + Cu/SA is µ hnf , ρCp hnf is the Al 2 O 3 + Cu/SA heat capacity, k hnf is the Al 2 O 3 + Cu/SA thermal/heat conductivity, ρ hnf the density of Al 2 O 3 + Cu/SA , and δ hnf the electrical conductivity Al 2 O 3 + Cu/SA .The aluminium oxide (Al 2 O 3 ) thermophysical properties, along with copper (Cu) and Sodium alginate (SA) nanoparticles, are revealed in Table 1.Equation 8 contains the thermophysical properties of Al 2 O 3 + Cu/SA .Here φ 1 and φ 2 are volume fractions of Al 2 O 3 and Cu respectively.Further, if α 1 = 0 , then, Eq. 4 is converted into a Newtonian fluid.