Heat and mass transfer analysis of assisting and opposing radiative flow conveying ternary hybrid nanofluid over an exponentially stretching surface

Access to dependable and environmentally friendly energy sources is critical to a country's economic growth and long-term development. As countries seek greener energy alternatives, the interaction of environmental elements, temperature, and sunlight becomes more critical in utilizing renewable energy sources such as wind and bioenergy. Solar power has received much attention due to extraordinary efficiency advances. under this context, the present work focus on solar radiation and chemical processes in the presence of modified ternary hybrid nanofluids (THNFs) circulating over an exponentially stretched surface in both aiding flow (A-F) and opposing flow (O-F) circumstances. The primary objective of this investigation is to dive into the complicated dynamics of these structures, which are distinguished by complex interactions involving radiation, chemical reactions, and the movement of fluids. We construct reduced ordinary differential equations from the governing equations using suitable similarity transformations, which allows for a more in-depth examination of the liquid's behavior. Numerical simulations using the Runge–Kutta Fehlberg (RKF) approach and shooting techniques are used to understand the underlying difficulties of these reduced equations. The results show that thermal radiation improves heat transmission substantially under O-F circumstances in contrast to A-F conditions. Furthermore, the reaction rate parameter has an exciting connection with concentration levels, with greater rates corresponding to lower concentrations. Furthermore, compared to the O-F scenario, the A-F scenario promotes higher heat transfer in the context of a modified nanofluid. Rising reaction rate and solid fraction volume enhanced mass transfer rate. The rate of thermal distribution in THNFs improves from 0.13 to 20.4% in A-F and 0.16 to 15.06% in O-F case when compared to HNFs. This study has real-world implications in several fields, including developing more efficient solar water heaters, solar thermal generating plants, and energy-saving air conditioners.


List of symbols
Science and advanced methods as well as recent technology have played a prominent role in the production of power 1 , refrigeration and heating 2 , production 3 , applications in medicine 4 , pharmaceutical industries 5 .To operate effectively and functionally of these systems will depends on excellent thermal management systems 6 .Further, improving the power efficiency and minimizing the components will challenge the existing techniques.
To overcome this a new concept of enhancing the thermal conductivity in working fluids by involving micro-scale nanoparticles was introduced by Hamilton and Crosser 7 .In 1995, Choi and Eastman 8 conducted pioneering research on using nanoparticles (NPs) to improve liquid thermal conductivity.The resulting liquid formed is coined as nanofluid (NFs).These fluid shows gradual improvement in its thermal conductivity due to this NFs applications are found in [9][10][11][12] .Based on the applications, we can find some of the advantages of nanofluids with relevant to base liquid.The rate at which heat is exchanged in thermal systems may be increased by utilising a nanofluid because of its better thermal conductivity relative to the base liquid.Therefore, by employing nanofluids to increase the temperature transfer rate, the size of a thermal system may be lowered, meaning a more compact system with savings in material weight and cost.Nanofluids are more stable and can improve the transmission of heat more than carrier fluids can.In view of this, numerous researchers have examined and evaluated the usage of nanofluids in thermal transfer applications.
Sandeep et al. 13 analyzed the flow characteristics of a chemically reactive Casson liquid over a convectively heated curved area.The results demonstrate that non-Newtonian fluids, in contrast to Newtonian fluids, have a much higher heat transfer rate under unequal heat sources and viscous dissipation conditions.Yasir et al. 14 looked at the dynamics of ethylene glycol transporting copper and titania NPs over a stretchable/shrinkable curved structure, including a stability in their study.Because copper enhances porosity and titania functions as a photocatalyst, the study also emphasises the significance of these nanocomposites.Prasannakumara et al. 15 explored the TPD (Thermophoretic particle deposition) in a bioconvective NFs flowing over an exponentially extended surface.The findings reveal that adding NPs in a carrier fluid will increases the heat transformation rate.Poojari et al. 16 evaluated improved conveyance of heat in unsteady magneto-NFs circulation caused by a stochastic expanding surface with convective boundaries.They employed the Maxwell and Xue models of NFs in their research.According to the results, the Maxwell nano model predicts a far higher rate of thermal circulation than the Xue nano model predicts.Khan et al. 17 evaluated the effect of mixed convection and radiation on the temperature transfer of a nanofluid in a slip flow through a bending sheet subject to activation energy and a binary reaction.Results reveal that the addition of nanoparticles slows the flow rate for both upper and lower branch solutions.Some of the works with different kinds of nanoparticles in base liquds are listed in [18][19][20][21] .
A specialised group of NFs known as HNFs (hybrid nanofluids) has arisen in recent years.An innovative kind of nanofluid, HNFs are made by dispersing many NP types in a single solvent.A specific material may not contain all desirable features necessary for a certain application; it may have excellent thermal or rheological capabilities.An owing to the synergistic effect, the HNFs is projected to have higher thermal conductivity than separate nanofluids.Recently, Nanda et al. 22 studied the nonlinear/linear expandable surface containing joule impact, radiation with three dimensional movement of tangent hyperbolic liquid containing aluminium alloys.The influence of nanoparticles, the heat transfer rate in the nonlinear extended scenario is double that of the linear stretched case.Ramesh et al. 23 examined the HNFs circulation across a nonlinear/linear stretched surface by considering TPD.In the study, linear situation, including nanoparticles augments heat propagation but decreases concentration while improving axial velocity in the nonlinear case.Yasir et al. 24,25 explored the thermal performance of different kind of nanoparticles (SWCNT, CuO, MgO, Ag) with H 2 O and EG carrier liquids.Sulochana et al. 26 examined the tangent hyperbolic liquid containing Mgo and Cuo nano sized particles in the presence of magnetic field.In the nanoparticle volume factor, the energy transmission rate in the nonlinear extending case is greater than the linear case.
Researchers are also investigating the inclusion of different NPs into HNFs, resulting in modified nanofluids.This novel class of functional liquids shows promise and is the topic of much investigation.Trihybrid nanofluids (TNFs), made up of three solid NPs combined with a carrier liquid, have been studied for their improved thermal conductivity and thermal expansion.Recently, Yook et al. 27 used a multi-linear regression model to examine ternary hybrid nanoparticles' thermal and momentum transmission in a channel with varying permeabilities and porous walls.In their study they used three different cases of nanoparticle combination to examine the thermal performance.The results show that more rapid heat transmission is shown in mixtures of paraffin wax, sand, and AA7072.Madhukesh et al. 28,29 examined the TNFs flow to examine the thermal and mass transfer analysis.The TNFs shows greater performance in these two aspects than HNFs and NFs.Animasaun et al. 30 studied thermal examination of magnitudes of different kinds of nanoparticles with water as a base liquid.More thermal significant performance is achieved with smaller densities of the nanoparticles.
Aerodynamics, plastic sheets preparation, processing of materials, wire drawing and environmental flows are some of the examples which draw the attention of exponential stretching sheet (E-SS).The study also draws an attention the impacts of this geometry and the movement of the surface on motion of the fluid and thermal distribution.Due to this many works are carried out on E-SS geometries.Prasannakumara et al. 15 investigated the TPD, H-SS and movement of microorganisms over an E-SS using NFs.Alqahtani et al. 31 examined a change in energy and mass due to the motion of a Casson hybrid nanofluid over an elongated stretching sheet.In their study thermal and slippage of velocity circumstances, the absorption of heat, viscous dissipation, thermal radiation, the Darcy impact, and thermophoresis diffusion impacts are studied.Carreau flow of Cu-water nanofluids through an exponentially permeable stretched thin surface with an MHD thermal boundary layer is quantitatively studied by Yousif et al. 32 .Souayeh and Ramesh 33 studied the mobility of metallic ternary nanofluids (Ag-Au-Cu/ H 2 O), taking into account a wide variety of phenomena like gyrotactic organisms, energy of activation, buoyancy forces, and thermal radiation, as they enacted through an exponentially extending sheet.
The mechanics of fluids and the transfer of heat combine at the interface, where thermal radiation plays a crucial role in determining liquid behavior.Without a medium, heat may be transferred from one place to another by means of electromagnetic waves; this process is governed by thermal radiation (T-R).Intricate studies are required to effectively characterize and forecast liquid behavior when T-R is included in the governing equations.To create effective thermal systems, such as solar panels, burning chambers, or spacecraft re-entry scenarios, it is necessary to have a firm grasp of these interactions.Yasir et al. 34 looked into the effects of nonlinear T-R and non-uniform H-SS on the movement of ethylene glycol in the presence of hybrid nanomaterials like Silicon dioxide and Titanium oxide.The results indicate that in both variants of the optimal solution, heat transmission improves as a result of an increase in the radiation effect.The effect of electromagnetic radiation and convective slippery conditions on second-grade nanofluids with high viscosity via permeable medium was studied by Jamshed et al. 35 .The study reveals that increasing radiation levels modify its temperature distribution and the Nusselt number.The consequences of heat production and absorption in the T-R mixed convective circulation of a hybrid nanofluid through an inclinically contracting interface were studied by Yasir et al. 36 .the study reveals that boosting thermal dispersion through a rise in T-R and Eckert numbers.Some recent works on T-R with TNFs are given in 37,38 .
In a chemical reaction, two or more chemicals undergo a transformation to produce a new compound.Chemical reactions (C-R) form the backbone of the discipline and perform an essential role in many other areas of www.nature.com/scientificreports/science, industry, and daily life.They play a role in everything from combustion to metabolism to photosynthesis to the creation of new materials.The areas of health, materials research, energy generation, and environmental protection all stand to benefit greatly from a deeper understanding of, and ability to manipulate, chemical interactions.Khan et al. 39 examined the time-dependent movement of non-Newtonian fluid, emphasizing thermal and solutal movement.The results demonstrate that an upsurge in the surface concentration of the catalyst enhanced the efficiency of both homogeneous and heterogeneous techniques.The impact of endothermic/exothermic C-Rs with activation energy on a wedge-shaped ternary hybrid nanofluid was studied by Sajid et al. 40 .The study reveals that C-R will upsurges rate of thermal distribution.Using a squeezed parallel infinite plate as a flow channel, Bilal et al. 41 studied the characteristics of a fluctuating electroviscous TNFs.The study shows that as C-R values rise, the pace at which mass is allocated rises.
In all the above served literatures, the studies individually examined the impacts of T-R, chemical reaction, and porous medium over different geometries.No work is found to be studied by integrating all the abovementioned impacts over E-SS in the presence of TNFs.Analyzing the resultant complicated system of ordinary differential equations requires the use of innovative numerical techniques like the Runge-Kutta Fehlberg (RKF) method and shooting algorithms.The major engineering coefficients are also covered in this study.The results have real-world applications in areas such as thermal engineering, manufacturing processes, and renewable energy infrastructure.Finally, the aim of the article provides answers to the following pertinent research questions: • How do several important variables affect both A-F and O-F cases?
• What effect do changes in various parameters have on the rate of heat transfer and skin friction?• How change in chemical reaction constraint will interact with mass transfer rate?

Mathematical modelling
The steady, two-dimensional, laminar and incompressible flow of THNF (combination of Al 2 O 3 + TiO 2 + Ag nanoparticles and H 2 O base fluid), see Fig. 1.The sheet is stretching along x-direction and y-axis normal to it.The scenario of mutual assistance between the stretching-generated flow, thermally buoyant flow, and as well as the scenario of mutual antagonism between these two flows are studied.The sheet uniform velocity is represented by U = U * e x l in the corresponding x− axis direction.Further, the sheets variable temperature is provided by T w = T ∞ + T 0 e x 2l .Here, T w ,T ∞ and T 0 denotes the wall, ambient and reference temperature of the sheet, respectively.It is also assumed that the variable concentration of the sheet is In which; C w ,C ∞ and C 0 denotes the respective wall, ambient and reference concentration of the exponential surface of the sheet.However, δ r is the variable chemical rate introduced into concentration equation 42 and it is defined as δ r = 0.5 × δ 0 e x l .Meanwhile, at y = 0 , the temperature and concentration at the surface of E-SS is denoted by T w and C w , respectively.The far field boundaries as y → ∞ are presented by the respective T ∞ and C ∞ .Based on these assumptions with boundary layer and Boussinesq approximations, the governing equations for the present problem are as follows: (see [43][44][45][46] ): (1) The last term in Eq. ( 2) specifies the effect of buoyancy force which has a negative and positive sign, where this positive sign refers the buoyancy assisting flow and the negative sign refers the case of buoyancy opposing flow.Moreover, in Eq. ( 3) the term q r is given by q r = − 4σ * 3k * ∂T 4 ∂y (Rosseland approximation), see 47 .Additionally, expressing the term T 4 as a linear function of T and ignoring the higher order terms, we get Finally, the equation become as q r = − 16σ * 3k * T 3 ∞ ∂T ∂y .By using the below suitable similarity (see 44 ): Using Eq. ( 7) in Eqs.(2-4), one obtains the form as: with boundary conditions are: From the above equations the non-dimensionless parameters are tabulated in Table 1.Nanofluid and its correlations are discussed by many researchers.The first model of thermal conductivity of the nanofluids proposed by Maxwell model 48 .Later, Yu and Choi 49 considered the nanoliquid layer and proposed new thermal conductivity.By considering Brownian motion and aggregation, Xuan 50 proposed new (

Sl. No
Symbol Name of the parameter Special case thermal conductivity of the nanoliquids.Some of the other thermophysical correlations of the nanofluids are given in the works of 4,[51][52][53] .Thermophysical properties of TNFs used in the present study are given as follows (see 54 ) In the above expressions when 3 = 0 , the properties reduce to HNFs and in the absence of 3 and 2 , the properties reduce to requisite posited NFs.

Gradients.
The important engineering factors such as Cf x , and Nu x are described as follows: and The expression for Sh x (Sherwood number) is given as: In the above equations, τ w , q w and q m are expressed as, By using (15) in (13-14), we get where Re = lU * e x l 2ν f is called the local Reynolds number.

Methodology
The employing computational approaches, where the simplified governing Eqs.(8-10) and BCs are addressed using the RKF-45 55 approach and shooting procedures 56 .We turn the revised equations into a first-order system via introducing new variables.In order to transform the system of equations into first order, we will select, f = r 1 , f ′ = r 2 , f ′′ = r 3 , θ = r 4 , θ ′ = r 5 , χ = r 6 and χ ′ = r 7 .Hence, the equations become as: To solve the IVP defined by Eqs. ( 20)- (22) and BCs in (23), a quantitative approach called the RKF-45 order method is applied.The shooting procedure is utilized with carefully selected error tolerance of 1/10 6 and step size value of h 1 = 1/100 correspondingly, to ensure that BCs at infinity are satisfied.For obtaining mathematical estimates, the built-in MATLAB program called bvp4c solver is implemented.It involves substituting parameters K 1 = 0.1, Rd = 1, Sc = 0.8, Rr = 0.1 and = ±0.3, considering properties mentioned in Table 2. Graphical rep- resentation of the results generated for each constraint by varying each parameter while keeping the remaining values constant.In order to validate the numerical code for −θ ′ (0) , a comparison is made between the findings of present study and formerly published research.The findings, as presented in Table 3, demonstrate a satisfactory level of agreement between the outcomes of the two datasets.
The algorithm of the RKF -45 is given below (see 55 ): The current methodology can be implemented by the following six steps: RK 4 th order : Table 2.The NPs and H 2 O thermophysical properties taken from the works of 54 .

Results and discussion
The flow of THNF across an ESS in the presence of a permeable medium is examined in this work.The modelling takes into account the impact of chemical reactions and renewable (solar) radiation.The dimensionless equations are formed from a collection of suitable similarity variables to simplify the consider task.Several critical elements are uncovered using this technique, and their impact on the profiles is visually depicted and briefly described.
The variation of K 1 on f ′ , θ , and χ in both A-F and O-F situations are shown in Figs. 2, 3, 4. In both flow circumstances, the velocity f ′ falls as K 1 rises (see Fig. 2).This is happening due to the reason that the greater frictional force caused by higher porosity values, the system is more resistive and fluid flow is constrained.Potential drop in the fluid velocity is observed, and this results in higher resistance.Additionally, when K 1 rises, the fall in f ′ (η) for O-F is more noticeable than the case of A-F.In both situations, higher porosity values also lead to higher θ profiles.Temperature is improved by increasing K 1 because it thickens the thermal boundary layer (TBL).This finding is consistent with the reasonable assumption that higher K 1 improves heat retention and transmission within the liquid.In addition, when K 1 is present, heat transmission in circumstances with O-F is greater than in situations with A-F.Rise in the values of K 1 will improves the concentration for both A-F and O-F cases as shown in Fig. 4. Improvement in K 1 will improves the frictional force which results improvement in the concentration boundary layer (CBL).Consequently, the level of concentration shows significant improvement for the O-F case than the A-F case. Figure 5 reveals the variation of Rd on thermal profile.An increase in Rd will enhances thermal distribution for both A-F and O-F cases.Rd has an inverse relationship with the mean absorption coefficient, which drops as the value of Rd increase.The buoyancy force is reduced, as a result, of conductive heat transmission, which proves to be more efficient compared to radiative heat transfer.In fact, a higher thermal dispersion is achieved with increased heat transfer to operating fluids through a higher Rd value.Furthermore, gradual improvement in temperature of fluid is observed as the Rd values are set to be high, and the fluid becomes more heated.Following that, heat transmission significantly increases.In the presence of Rd , the O-F shows high thermal distribution than A-F case.
The impact of Sc on the χ profile is observed in Fig. 6.In both scenarios, increasing Sc lowers the concentra- tion.Increased values of Sc improves the mass diffusion coefficient, which lead to a reduction in concentration.This implies that the fluid's capacity to convey velocity is comparatively more effective than its capacity to convey mass or concentration.Consequently, the predominance of fluid motion supersedes mass transfer, leading to less efficient dispersion of concentration gradients arising from the origin.In comparison to the O-F situation, concentration is reduced in the A-F case.
Figure 7 shows how the variance of χ changes with various Rr values.In both situations, increasing values cause the concentration to drop.Higher Rr indicates that the transformation of reactants into products is taking     www.nature.com/scientificreports/place at an accelerated rate causes lower CBL.When Rr is present, the O-F case concentrates less than the A-F case.
Figures 8 and 9 indicate how the key engineering interests of Cf x and Nu x vary with different dimensionless restrictions.Figure 8 depicts the fluctuation of Cf x on K 1 for various values of 3 for both flow situations.The graph shows that increasing the value of 3 decreases the friction drag force in both flows.The increase in K 1 values will resist the flow of the liquid.The addition of K 1 to the scale of 3 enhances the MBL thickness.This results in a reduction in surface drag force.In the O-F situation, Cf x is much lower than in the A-F example.Figure 9 depicts the fluctuation of Nu x on Rd as 3 values vary.The rise in the value of 3 and Rd will increase the thermal transfer rate.
The elevation in the values of Rd improves the rate of heat transfer due to presence of k * .Increment in 3 also improves the TBL which results in improvement of Nu x .Nu x is more in A-F case than O-F case.
Figure 10 displays the variation of Sh x with change in the values of 3 and Rr for both A-F and O-F scenarios.As the values of Rr escalates, the rate of transformation of reacting substances into outputs is accelerated further.On the other hand, 3 will also improves the thickness of the CBL due to improvement in the surface area.A-F case shows greater rate of mass distribution than O-F.
Further, Table 4 shows the computational values of Cf x ,Nu x and Sh x for both A-F and O-F cases in the pres- ence of various dimensionless constraints.Table 5 displays the change in the percentages of Cf x ,Nu x and Sh x for both A-F and O-F cases in comparison with TNFs and HNFs.The tabulated outcome shows that TNFs shows gradual improved performance in all the constraints in comparison with HNFs.
The rate of change in the percentage of skin friction, Nusselt number and Sherwood number in the presence of TNF and HNF are presented in Table 5.To obtain Cf x %,Nu x % and Sh x % the formulas are given below:

Final remarks
The idea of this study is to investigates the transfer of heat and mass on ESS by considering the influence of A-F and O-F with ternary hybrid nanofluid.It also incorporates T-R, chemical reaction, and porous medium.Utilizing RKF-45 the numerical solutions are found for reduced ODEs.Graphs are utilized to illustrate the important dimensionless limitations.The key findings of this work are as follows: • The velocity profile declines with higher impacts of the porosity parameter while the temperature and con- centration profiles boosted up.• The augmentation of the solar radiation parameter decelerates the velocity, thermal distribution, and tem- perature gradients.This indicates that reducing thermal radiation flux decelerates the thickness of TBL while increasing the thickness of the MBL.• For increasing thermal radiation and 3 , the heat transfer rate augmented.In general, the radiation parameter enriches the thermal conductivity, as a consequence, the heat transfer improves.• The shear stress and mass transfer rate decelerate with superior influences of 3 but the heat transfer rate remarkably uplifts.• When compared to HNFs, the rate of heat distribution in THNFs goes from 0.13% to 20.4% and from 0.16% to 15.06%, respectively, for the case of A-F and as well as for the case of O-F.
The current research is relevant in several applications like capacitors, biofuel, batteries, nanomaterials, and power storage etc.By captivating the following distinct impact such as the ternary nanofluid flow over an inclined stretching/shrinking surface, mass suction/injection, convective boundary conditions, Two-phase model and Newtonian heating is also possible to consider in the extended work of this problem.

1 .
Figure 1.Geometric representation of the problem.

Figure 8 .
Figure 8. Variation of Cf x on K 1 for numerous values of 3 .

Figure 9 .
Figure 9. Variation of Nu x on Rd for numerous values of 3 .

Figure 10 .
Figure 10.Variation of Sh x on 3 for numerous values of Rr.

Table 3 .
Validation is performed by comparing the results with −θ ′ (0) values for a selection of the limiting cases.

Table 4 .
Computed numerical results of Cf x , Nu x and Sh x for several constraints which are dimensionless with 1 = 2 = 0.01.

Table 5 .
Comparative analysis of change in Cf x %,Nu x % and Sh x % for dimensionless constraints in the case of A-F and O-F.