A Group Theoretic Analysis of Mutual Interactions of Heat and Mass Transfer in a Thermally Slip Semi-Inﬁnite Domain

: Group theoretic analysis is performed to get a new Lie group of transformations for non-linear differential systems constructed against mass and heat transfer in the thermally magnetized non-Newtonian ﬂuid ﬂow towards a heated stretched porous surface. The energy equation is used with additional effects, namely heat sink and heat source. The chemical reaction is also considered by the use of the concentration equation. The symmetry analysis helps us in numerical computations of surface quantities for (i) permeable and non-permeable surfaces, (ii) thermal slip and non-thermal slip ﬂows, (iii) magnetized and non-magnetized ﬂows, (iv) chemically reactive and non-reactive ﬂows. For all these cases, the concerned emerging partial differential system is transformed into a reduced ordinary differential system and later solved numerically by using the shooting method along with the Runge-Kutta scheme. The observations are debated graphically, and numerical values are reported in tabular forms. It is noticed that the heat transfer rate increases for both the thermal slip and non-slip cases. The skin friction coefﬁcient declines towards the Weissenberg number in the magnetized ﬁeld.


Introduction
Materials processing, power generation, transportation, civil infrastructure, food production, automobiles, and hydroelectric power plants, to mention just a few, are the dominant parts of fluid mechanics applications.Owning such a vital role in fluid science, mathematicians, oceanographers, geologists, biologists, atmospheric scientists, physicists, and engineers pay attention to fluid flow fields.For example, the heat transmission characteristics of a visco-elastic fluid flow by way of a stretched surface carrying a heat sink and heat source were examined using mathematical analysis by Vajravelu and Rollins [1].In terms of Kummer's and parabolic cylinder functions, solutions for heat transfer rate and temperature towards the Prandtl number were obtained.For small Prandtl numbers, it was demonstrated that there was no boundary layer in the solution.A mathematical analysis of heat and momentum for flow towards a stretched surface was performed by Andersson et al. [2].Accurate similarity transformations were carried out to diminish the flow of time-dependent equations to ordinary differential equations.For certain iterations of an unsteadiness parameter and Prandtl number, the resulting problem was numerically solved.The temperature grew monotonically towards the ambient field.The visco-elastic heat and flow properties subject to the porous stretched medium were considered by Abel et al. [3].Temperature-dependent viscosity was considered.Such assumption results in non-linear equations.Such involved equations were solved numerically.The impact of visco-elastic permeability parameters and fluid viscosity towards numerous conditions were investigated for two separate scenarios, namely heat flux (PHF) and prescribed surface temperature (PST).In this examination, the most noteworthy discovery was that the skin friction admits an inverse relation towards the permeability parameter.The magnetized heat transmission over a stretched sheet in the visco-elastic fluid was investigated by Zakaria [4].The successive approximation approach was used to get the solution of flow equations.On the temperature and velocity, the impacts of relaxation time parameter, Prandtl number, Alfven velocity, surface mass transfer, and elastic velocity coefficients were debated graphically.The heat and flow aspects of a second-grade liquid towards a stretched surface were investigated by Cortell [5].The order reduction of the involved flow equations was conducted via similarity transformations.The energy equation with viscous dissipation was considered, and the variation in both temperature gradient and temperature was investigated.Abel and Mahesha [6] investigated magnetized a visco-elastic liquid towards a sheet with the assumptions of a heat source that is non-uniform with thermal radiation.The thermal conductivity was believed to change with temperature in a linear fashion.The primary flow equations were partial differential equations (PDEs), which were later reduced as ordinary differential equations (ODEs) by using the appropriate transformations.The regular perturbation scheme was used to solve the altered equations.The efficient shooting method also yielded a numerical solution, which agrees well with the analytical answer.In several plots, the various flow variables that determine temperature profiles, such as heat sink/source, thermal radiation, and visco-elastic parameters, Eckert, Prandtl, and Chandrasekhar numbers, were shown, and the rate of heat transfer coefficients was also measured against these parameters.By considering heat sink/source, viscous dissipation, and magnetic field, an analysis was performed to examine the convective flow having an exponentially stretched plate and temperature regime.The extremely non-linear momentum and energy equations have approximate analytical similarity solutions.The current findings admitted great accord with the existing work on a variety of special instances.For various values of the governing parameters, heat transfer rate and temperature regimes were examined.On the directional flow coordinate, numerical solutions were derived towards exponentially stretched velocity.The implications of numerous physical parameters on dimensionless heat transfer characteristics, such as Prandtl number, Hartman number, and Grashof number, were thoroughly examined.It was discovered that increasing the Prandtl number lowers the drag faced by fluid while increasing the magnetic field strength raises the local Nusselt number.By considering heat transfer, volumetric heating, magnetic field, and a variety of other factors, the flow of a visco-elastic fluid over a porous stretched surface was considered by Pal [7].It was demonstrated that their solutions do not appear to have emerged earlier and can be developed by making careful picks of the involved functions.On a stretching sheet, the pseudoplastic flow field was investigated by Ashrafi and Meysam [8].The flow equations were converted into lower-order equations via transformations.After that, the system was numerically integrated with the RK scheme.The impact of various flow variables, such as heat transfer coefficient, apparent viscosity, temperature, velocity, towards Prandtl number, pseudo-plasticity index, and unsteadiness parameter, were studied in depth.The stretched surface was used by Zhang et al. [9] to explore heat transfer aspects of Oldroyd-B fluid with suspended nanoparticles.Here, with Ag and Cu, the polyvinyl alcohol is treated as a base fluid.The stretched sheet velocity and temperature were supposed to vary.The flow equations were first constructed, and the similarity transformation was used to convert them to ODEs.The homotopy analysis method (HAM) was used to derive the analytical answers, which exhibit good agreement with earlier results.For analyzing the heat transmission of fractional visco-elastic magnetohydrodynamics (MHD) fluid flow field, the Lie group was introduced by Chen et al. [10].The Grünwald scheme approximation was used to minimize and solve fractional equations conjectured with Riemann-Liouville operators numerically.The results reveal that the wall stretching exponent, fractional derivative, and magnetic field all have a significant impact on skin friction and heat conductivity.For high fractional-order derivatives, visco-elastic fluids move faster and do not stick near the outer flow.As the magnetic field parameter is increased, skin friction increases, while heat transmission decreases.Refs.[11][12][13][14][15][16][17][18][19][20][21][22][23][24][25][26][27][28][29][30] can be used to examine past and recent trends.Besides this, heat transmission and transport phenomena in porous media are significant processes in many technical applications, including heat pipe technology, chemical catalytic reactors, electronic cooling, pack-sphere bed, and heat exchangers, to name a few.Heat transport in the porous fibrous medium is a well-known and complicated subject that has received a lot of attention.Many drying and heating uses of such materials require a detailed grasp of this challenge.Similarly, several textile materials and, more recently, ovens that are utilized as heat barrier materials are being investigated.As a result, due to its ever-increasing applications in industries and modern technology, the study of flow and heat transfer in porous media exposed to non-Newtonian and Newtonian fluid flows has gotten a lot of attention [31][32][33][34][35].
By considering mathematical formulations, their solutions, and the motivations carried by the researchers reported above, we propose a new symmetry transformation for non-Newtonian thermally magnetized flow subject to both non-permeable and permeable sheets.The effort in design is explained in the various sections.In Section 1, the limited literature survey on mathematical formulation subject to fluid flow fields is reported.In Section 2, the flow formulation is disclosed.The procedure to obtain a scaling transformation is debated in Section 3. Section 4 is devoted to sharing the procedure adopted for numerical solutions.The observations are debated in Section 5 with the parallel graphical outcomes.The important results are itemized in Section 6.We conclude that the new scaling group of transformations will help researchers implement such exercise on their problems for batter narration of flow problems by using symmetry analysis with mathematical modeling.

Flow Formulation
The flow of Williamson fluid (WF) is considered over a heated flat surface.The surface was taken as porous, and an external magnetic field was applied perpendicular to the flow.The heat transfer in fluid with heat source/sink effect was entertained, and concentration aspects of fluid were taken in the presence of a first-order chemical reaction.At the surface, both temperature and velocity slip effects were also assumed.The physical model of the flow problem is given in Figure 1.The present physical model was studied by using the continuity equation, Cauchy momentum equation, energy equation, and concentration equation.The said ultimate equations are: To obtain a system of differentials for our problem, we needed the particular rheology of the Williamson fluid model [36], which could be assessed using the following constitutive relation for the Williamson fluid model: Ω is given as: . To obtain a system of differentials for our problem, we needed the particular rheology of the Williamson fluid model [36], which could be assessed using the following constitutive relation for the Williamson fluid model: where  is given as: For the present problem we considered 0 and 1. gives: Carrying the binomial expansion in Equation (3), we get: ( ) The Williamson fluid velocity . Further, by considering assumptions of steady flow, the external magnetic field effect, heat sink/source effect, chemical reaction effect, and using the Williamson fluid constitutive relation (Equation ( 5)) into Equation (1), resulted in the following coupled non-linear flow narrating differential equations for our physical problem: For the present problem we considered µ ∞ = 0 and Γ .
Ω < 1.Hence, Equation (2) gives: Carrying the binomial expansion in Equation (3), we get: The Williamson fluid velocity U 1 is taken along the X 1 -axis and U 2 is towards X 2 -axis.Further, by considering assumptions of steady flow, the external magnetic field effect, heat sink/source effect, chemical reaction effect, and using the Williamson fluid constitutive relation (Equation ( 5)) into Equation (1), resulted in the following coupled non-linear flow narrating differential equations for our physical problem: along with boundary constraints: Further, we used: and obtained: ∂U ∂X while the corresponding equations are: The partial differential Equations ( 12) to ( 15) form a structure.These equations are strongly non-linear and coupled.So far, finding an exact solution was not possible.Therefore, we reduced the independent variables by using single parameter transformations for examination.Rather than continuing with transformations accessible in the literature, we found some specifically for the current problem.We needed the stream function relation: Equations ( 12)-( 16) under Equation ( 17) gets the forms: the concerned endpoint conditions are:

Symmetry Exploration
The Equations ( 6)-( 9) are a series of partial differential equations that are strongly non-linear and coupled.Our aim was to decompose the equations into independent variables.For this, we needed a series of scaling transformations.As a result, the following one-parameter transformation category can be considered: (23) here, the old coordinates (C, T, Γ, Ψ, Y, X) were replaced with (C * , T * , Γ * , Ψ * , Y 1 , X 1 ) due to Equation (23).In this case, the implementation results in: carrying invariant condition for Equations ( 24)-( 26) under G s , we get: further (boundary conditions) BCs give λ 4 = 0, and λ 6 = 0. Equation ( 27) results in: Equation ( 23) under Equation ( 28) offers: Equation ( 29) with the use of Taylor's expansion around ε = 0 up-to O(ε), we get: Therefore, we have: The use of Equation (31) results in: Equations ( 19)- (22) under Equation (32) gives: while the reduced endpoint conditions are: The Williamson fluid fills a semi-infinite domain.As a result, the Sherwood number (ShD), Skin Friction Coefficient (Skin-FC), and Nusselt number (Nm) are among the surface quantities.These surface quantities have the following mathematical expression: the corresponding dimensionless forms can be written as: and the flow parameters are given as: where Wb, Mg, Pr, Hr, Sc, Vp, Ts, Pm, and Rs denote the Wessinberg number, magnetic field parameter, Prandtl number, heat generation parameter, Schmidt number, velocity slip parameter, thermal slip parameter, porosity parameter, and chemical reaction parameter, respectively.

Solution Procedure
We took into account the Williamson fluid flow towards a stretched porous surface in this paper.In the existence of heat absorption, heat generation, and chemically reactive species, the slip flow field was thermally magnetized.The assumptions are mathematically represented as PDEs.The use of a Williamson fluid with many physical special effects resulted in highly non-linear PDEs, making an exact solution unattainable at this time.When an exact answer appears to be impossible, we always look for a numerical solution [37,38].Researchers in this subject, in particular, used a variety of approaches to report numerical solutions.The majority of methods for investigators seeking a numerical solution are by transforming the PDEs into ODEs, which are then solved by the numerical scheme.It is worth noting that the conversion from PDEs to ODEs is accomplished through a specific set of transformations [39,40].Many scholars consider such transformations straight from the literature rather than obtaining specific transformations of the flow problem.This does not allow for a better narration of the flow fields.In our situation, we used symmetry analysis to provide transformations for equations.We created coupled ODEs using these transformations and solved them using a shooting method combined with the RK arrangement.It is worth noting that we employed self-coding to implement the shooting approach by converting Equations ( 33)-( 36) into an initial value problem.The reduced initial value problem was then transformed into a system of seven first-order differential equations, which were solved by selecting three missing conditions as appropriate starting estimates.The computed solution converges if the boundary residuals are less than the tolerance error 10 −6 .If the calculated results do not meet this criterion, the starting estimates are adjusted using Newton's method, and the operation is repeated until the solution meets the chosen convergence threshold.Further, ShD, Nm, Skin-FC, WF concentration, WF velocity, and WF temperature are considered as the quantities of interest.Here, we have eight differential flow variables, namely Wb, Vp, Mg, Pm, Ts, Hs, Pr, and Hr.The effects of such flow variables are inspected on WF temperature, WF velocity, WF concentration, Skin-FC, Nm, and ShD.The range of flow variables was chosen in such a way that the method's convergence and stability are maintained.The final observations in this direction are shared by using line graphs and tables.The ShD, Nm, and Skin-FC were the surface quantities of concern.These quantities were estimated by taking into account a variety of physical frames.The variation in Skin-FC for non-slip and slip flows is shown in Table 1 with positive Wb values.When Wb = 0.1, 0.2, 0.3, and 0.4 were increased for Vp = 0, the Skin-FC decreased intensely.Furthermore, the Skin-FC was found to be a diminishing function of Wb when Vp = 0.5.For iteration in porosity parameter, Table 2 provides the Skin-FC outcomes for both magnetic and non-magnetized cases.When Mg = 0, the Skin-FC was found to be rising as a function of Pm in an absolute sense.Furthermore, when Pm = 0.1, 0.2, 0.3, and 0.4 were increased for Mg = 0.5, the Skin-FC displayed inciting values.The Skin-FC was found to be an increasing function of Pm for both fields.Skin-FC numerical values for non-permeable and permeable sheets are listed in Table 3.In this scenario, we looked at the magnetized flow on a case-by-case basis, namely permeable surface (Pm = 0) and non-permeable surface (Pm = 0.5).When we increased Wb = 0.1, 0.2, 0.3, and 0.4 for Pm = 0, the Skin-FC decreased considerably.Furthermore, at Pm = 0.5, Skin-FC was found to be a decreasing function of Wm.Physically, the resistance encountered by particles decreased as Wb grew.

Results Analysis
The WF was fitted above the permeable magnetized surface.In the existence of heat generation and absorption special effects, the thermal flow regime was carried out.Velocity and thermal slip were also taken into account.Consideration of the concentration equation, as well as the chemical reaction, adds to the novelty.The problem was investigated numerically.In detail, the impact of Pm, Mg, Vp, and Wb were investigated on WF velocity and shown with the help of Figures 2-5, respectively.In detail, the effect of Pm = 0.0, 0.3, 0.6 on WF velocity was investigated, and observations in this way are shown in Figure 2. It is worth noting that when Pm = 0.0, the fluid flow is across a non-permeable surface, and the velocity strength of WF is much higher than when Pm = 0.3 and Pm = 0.6.Overall, it is seen that velocity declines as Pm increases.The porosity parameter possesses an inverse relation towards permeability of the porous medium; therefore, increases in the porosity parameter cause a decline in the stretching rate of the flat surface.Since we considered the no-slip condition, a reduction in stretching rate resulted in lower values of the velocity of WF. Figure 3 offers the variations Mg = 0.0, 0.5, and 0.9, which are all positive values.Here, Mg = 0.0 denotes the absence of a magnetic field, and velocity has a large value because there is no Lorentz force.The WF velocity is retarded for Mg = 0.5 and Mg = 0.9 as compared to Mg = 0.0.When Mg = 0.9 was compared to Mg = 0.0 and 0.5, the fall in WF velocity magnitude was greater.It has been noticed that as Mg is increased, the WF velocity decreases.One should note that the non-zero value of the magnetic field parameter gives birth to Lorentz forces.Here, we simulated Mg = 0.0, 0.5, and 0.9.The positive values 0.5 and 0.9 enhance the strength of the Lorentz force, and due to the resistive nature of the Lorentz force, the WF particles met higher resistance; this led to a decline in WF velocity.For Vp = 0.0, 0.2, and 0.4, the velocity distributions subject to the permeable magnetic surface are shown in Figure 4. We established non-slip Williamson fluid flow with Vp = 0, and the WF velocity magnitude was larger than with Vp = 0.2 and 0.4.We discovered that WF velocity over a magnetic surface decreased as Vp increased.Variations in Wb = 0.1, 0.3, and 0.7 are shown in Figure 5  It is worth noting that when Pm = 0.0, the fluid flow is across a non-permeable surf and the velocity strength of WF is much higher than when Pm = 0.3 and Pm = 0.6.Ove it is seen that velocity declines as Pm increases.The porosity parameter possesses an verse relation towards permeability of the porous medium; therefore, increases in the rosity parameter cause a decline in the stretching rate of the flat surface.Since we con ered the no-slip condition, a reduction in stretching rate resulted in lower values of velocity of WF. Figure 3 offers the variations Mg = 0.0, 0.5, and 0.9, which are all posi values.Here, Mg = 0.0 denotes the absence of a magnetic field, and velocity has a la value because there is no Lorentz force.The WF velocity is retarded for Mg = 0.5 and = 0.9 as compared to Mg = 0.0.When Mg = 0.9 was compared to Mg = 0.0 and 0.5, the in WF velocity magnitude was greater.It has been noticed that as Mg is increased, the velocity decreases.One should note that the non-zero value of the magnetic field par eter gives birth to Lorentz forces.Here, we simulated Mg = 0.0, 0.5, and 0.9.The posi values 0.5 and 0.9 enhance the strength of the Lorentz force, and due to the resistive nat of the Lorentz force, the WF particles met higher resistance; this led to a decline in velocity.For Vp = 0.0, 0.2, and 0.4, the velocity distributions subject to the permeable m netic surface are shown in         For both magnetized and non-magnetized situations, Table 4 was used to calculate the Skin-FC changes towards Wb.To begin, we assumed that the WF fluid was positioned above a non-magnetic (Mg = 0) permeable surface.When Mg was set to 0, we saw a drop in Skin-FC in absolute terms as Wb was increased.For magnetized WF flow across a stretched surface, the Skin-FC was also determined.Therefore, for Mg = 0.5, we discovered that Skin-FC decreased the function of Wb.In terms of physics, Skin-FC implied that the surface exerts an opposing drag force on fluid particles.As a result, increasing Wb = 0.1, 0.2, 0.3, 0.4 reduced the resistance faced by particles in both frames.Skin-FC numerical values for non-permeable and permeable sheets are presented in Table 5 for Mg = 0.1, 0.2, 0.3, and 0.4, respectively.Table 5 shows that in both non-magnetic and magnetized fields, the resistance provided by the surface increased as Mg values increased.Table 6 shows Nm variations for both magnetic and non-magnetized frames as they approach greater Pr values.When Pr = 1.1, 1.2, 1.3, and 1.4 were increased for Mg = 0, which is the non-magnetized condition, the Nm increased considerably.Similarly, when Mg = 0.5, Nm rose in response to greater Pr values.For both magnetized and non-magnetized situations, Table 4 was used to calculate the Skin-FC changes towards Wb.To begin, we assumed that the WF fluid was positioned above a non-magnetic (Mg = 0) permeable surface.When Mg was set to 0, we saw a drop in Skin-FC in absolute terms as Wb was increased.For magnetized WF flow across a stretched surface, the Skin-FC was also determined.Therefore, for Mg = 0.5, we discovered that Skin-FC decreased the function of Wb.In terms of physics, Skin-FC implied that the surface exerts an opposing drag force on fluid particles.As a result, increasing Wb = 0.1, 0.2, 0.3, 0.4 reduced the resistance faced by particles in both frames.Skin-FC numerical values for non-permeable and permeable sheets are presented in Table 5 for Mg = 0.1, 0.2, 0.3, and 0.4, respectively.Table 5 shows that in both non-magnetic and magnetized fields, the resistance provided by the surface increased as Mg values increased.Table 6 shows Nm variations for both magnetic and non-magnetized frames as they approach greater Pr values.When Pr = 1.1, 1.2, 1.3, and 1.4 were increased for Mg = 0, which is the non-magnetized condition, the Nm increased considerably.Similarly, when Mg = 0.5, Nm rose in response to greater Pr values.The WF temperature dependency on Pr, Hs, Ts, and Hr was examined and offered in terms of line graphs, as shown in Figures 6-9, respectively.To be more specific, the influence of Pr on temperature is depicted in Figure 6.The WF temperature dependency on Pr, Hs, Ts, and Hr was examined and offered in terms of line graphs, as shown in Figures 6-9, respectively.To be more specific, the influence of Pr on temperature is depicted in Figure 6.The WF temperature dependency on Pr, Hs, Ts, and Hr was examined and offered in terms of line graphs, as shown in Figures 6-9, respectively.To be more specific, the influence of Pr on temperature is depicted in Figure 6.In this case, the positive variation in Pr = 1.5,1.8,and 2.1 was taken into account.We have seen that when we have such a higher iteration, the temperature drops.Because Pr has an inverse relation with the fluid thermal diffusivity, we witnessed a decline in thermal diffusivity and, as a result, a considerable drop in WF temperature when we iterated Pr = 1.5, 1.8, and 2.1.The effect of Hs on WF temperature is depicted in Figure 7. WF temperature decreased when we iterated Hs = 0.0, −0.3, and −0.5.Iterations with Hs = 0.0, −0.3, and −0.5 caused energy loss and a drop in the overall temperature.The influence of the thermal slip parameter on WF temperature is offered in Figure 8. Ts = 0.0, 0.2, and 0.3 were employed in this line graph investigation.It is worth noting that we got a non-thermal case for Ts = 0.0.In this situation, the temperature appeared to be higher than Ts = 0.2 and Ts=0.3.Overall, the WF temperature appeared to be falling as Ts increased.Figure 9 depicts the impact of heat generation on the thermal flow regime of a WF fluid towards a   In this case, the positive variation in Pr = 1.5,1.8,and 2.1 was taken into account.We have seen that when we have such a higher iteration, the temperature drops.Because Pr has an inverse relation with the fluid thermal diffusivity, we witnessed a decline in thermal diffusivity and, as a result, a considerable drop in WF temperature when we iterated Pr = 1.5, 1.8, and 2.1.The effect of Hs on WF temperature is depicted in Figure 7. WF temperature decreased when we iterated Hs = 0.0, −0.3, and −0.5.Iterations with Hs = 0.0, −0.3, and −0.5 caused energy loss and a drop in the overall temperature.The influence of the thermal slip parameter on WF temperature is offered in Figure 8. Ts = 0.0, 0.2, and 0.3 were employed in this line graph investigation.It is worth noting that we got a non-thermal case for Ts = 0.0.In this situation, the temperature appeared to be higher than Ts = 0.2 and Ts=0.3.Overall, the WF temperature appeared to be falling as Ts increased.Figure 9 depicts the impact of heat generation on the thermal flow regime of a WF fluid towards a In this case, the positive variation in Pr = 1.5, 1.8, and 2.1 was taken into account.We have seen that when we have such a higher iteration, the temperature drops.Because Pr has an inverse relation with the fluid thermal diffusivity, we witnessed a decline in thermal diffusivity and, as a result, a considerable drop in WF temperature when we iterated Pr = 1.5, 1.8, and 2.1.The effect of Hs on WF temperature is depicted in Figure 7. WF temperature decreased when we iterated Hs = 0.0, −0.3, and −0.5.Iterations with Hs = 0.0, −0.3, and −0.5 caused energy loss and a drop in the overall temperature.The influence of the thermal slip parameter on WF temperature is offered in Figure 8. Ts = 0.0, 0.2, and 0.3 were employed in this line graph investigation.It is worth noting that we got a non-thermal case for Ts = 0.0.In this situation, the temperature appeared to be higher than Ts = 0.2 and Ts = 0.3.Overall, the WF temperature appeared to be falling as Ts increased.Figure 9 depicts the impact of heat generation on the thermal flow regime of a WF fluid towards a porous magnetic surface.The WF temperature increased dramatically when we iterated Hr = 0, 0.2, and 0.3.When Hr = 0.0, the thermal flow regime has no heat-producing impact, and the temperature of WF is smaller in magnitude than when Hr = 0.2 and 0.3 are utilized.This is because energy is generated inside the flow regime when we iterate the Hr, resulting in a rise in WF temperature.The impact of Sc and Rs is examined on WF concentration, see Figures 10 and 11.Particularly, Figure 10 deals with the impact of Sc on WF concentration.Sc was carried out in the following way; Sc = 1.8, 2.2, and 2.6.The WF concentration decreased as we iterated Sc.This effect was comparable to the Pr effect on WF temperature.As the Sc value increased, the mass diffusivity reduced, and the concentration of WF fell.The flow was considered with chemical reaction, and the resultant flow variable was the chemical reaction parameter Rs.In this direction we iterated Rs = 0.0, 0.2, and 0.4.It was found that for Rs = 0.0, the WF flow was non-reactive, and one can note that here, the concentration magnitude was higher as compared to Rs = 0.2 and Rs = 0.4.Together, the WF concentration showed declining values towards Rs = 0.0, 0.2, and 0.4.
Appl.Sci.2021, 11, x FOR PEER REVIEW 15 of 20 porous magnetic surface.The WF temperature increased dramatically when we iterated Hr = 0, 0.2, and 0.3.When Hr = 0.0, the thermal flow regime has no heat-producing impact, and the temperature of WF is smaller in magnitude than when Hr = 0.2 and 0.3 are utilized.This is because energy is generated inside the flow regime when we iterate the Hr, resulting in a rise in WF temperature.The impact of Sc and Rs is examined on WF concentration, see Figures 10 and 11.Particularly, Figure 10 deals with the impact of Sc on WF concentration.Sc was carried out in the following way; Sc = 1.8, 2.2, and 2.6.The WF concentration decreased as we iterated Sc.This effect was comparable to the Pr effect on WF temperature.As the Sc value increased, the mass diffusivity reduced, and the concentration of WF fell.The flow was considered with chemical reaction, and the resultant flow variable was the chemical reaction parameter Rs.In this direction we iterated Rs = 0.0, 0.2, and 0.4.It was found that for Rs = 0.0, the WF flow was non-reactive, and one can note that here, the concentration magnitude was higher as compared to Rs = 0.2 and Rs = 0.4.Together, the WF concentration showed declining values towards Rs = 0.0, 0.2, and 0.4.porous magnetic surface.The WF temperature increased dramatically when we iterated Hr = 0, 0.2, and 0.3.When Hr = 0.0, the thermal flow regime has no heat-producing impact, and the temperature of WF is smaller in magnitude than when Hr = 0.2 and 0.3 are utilized.This is because energy is generated inside the flow regime when we iterate the Hr, resulting in a rise in WF temperature.The impact of Sc and Rs is examined on WF concentration, see Figures 10 and 11.Particularly, Figure 10 deals with the impact of Sc on WF concentration.Sc was carried out in the following way; Sc = 1.8, 2.2, and 2.6.The WF concentration decreased as we iterated Sc.This effect was comparable to the Pr effect on WF temperature.As the Sc value increased, the mass diffusivity reduced, and the concentration of WF fell.The flow was considered with chemical reaction, and the resultant flow variable was the chemical reaction parameter Rs.In this direction we iterated Rs = 0.0, 0.2, and 0.4.It was found that for Rs = 0.0, the WF flow was non-reactive, and one can note that here, the concentration magnitude was higher as compared to Rs = 0.2 and Rs = 0.4.Together, the WF concentration showed declining values towards Rs = 0.0, 0.2, and 0.4.Tables 7 and 8 offer the Nusselt number outcomes for various parameters in different frames of reference.The Nm variations in thermal/non-thermal slip flows are reported in Table 7.For the non-thermal case (Ts = 0), we confirmed that when Pr = 1.1, 1.2, 1.3, and 1.4 increased, the Nm upturns.The Nm increased as Pr increased when Ts = 0.5.In thermal/non-thermal cases, the transfer of heat normal to the permeable surface was directly related to Pr.The transfer of heat normal to both permeable and non-permeable surfaces was apparent, see Table 8.We found that for both surfaces, the rate was an increasing function of Pr.Tables 9-11 were constructed to examine the variation in ShD towards three different frames, namely permeable and non-permeable surfaces, magnetic and non-magnetic fields, chemically reactive and non-reactive flows.The variation in ShD was observed for both non-permeable and permeable sheets (Table 9).The Pm = 0.0 implies the porous surface assumption, and for this case, we observed that when we increased Sc, the ShD increased reasonably.Together in both cases, porous and non-porous mediums, the ShD had a direct relation with Sc.Table 10 offers the ShD variation in both magnetic and non-magnetic fields towards Sc.We noticed that for both Mg = 0.0 and Mg = 0.5, the ShD showed inciting values towards Sc = 2.1, 2.2, 2.3, and 2.4.The impact of iteration in Sc on ShD for both chemically reactive and non-reactive flows was investigated and offered in terms of Table 11.In detail, Rs = 0.0 implies the case of non-reactive flow, and in this case, we noticed that ShD showed inciting values towards Sc = 2.1, 2.2, 2.3, and 2.4.Further, it can be seen that this impact was the same for chemically reactive flow.The magnitude of ShD variation was higher in the case of reactive flow.

Key Outcomes
The group-theoretic analysis was performed to determine the specific scaling transformations of the heat transfer problem, and through these specific transformations, we narrate the whole description of chemically reactive thermally magnetized Williamson fluid toward a stretched heated porous surface.Owing to the numerical solution, we arrived at the following conclusions:

•
WF velocity had inverse relations with Pm, Mg, Vp, and Wb.

•
WF temperature showed a declining nature with Ts, Pr, and Hs.

•
Higher values of Hr resulted in inciting values of temperature.

•
WF concentration was found with decreasing functions of Sc and Rs.

•
Nm showed inciting values towards Pr for thermal no-slip and slip regimes.

•
For magnetized and non-magnetized flows, ShD showed higher values of Sc.

•
In both chemically reactive and non-reactive cases, ShD increased as Sc increased.

•
Skin-FC showed declining values towards Wb in both non-magnetic and magnetic fields.
as a line graph of WF velocity.WF velocity has an inverse relationship with positive Wb values, as shown in the graph.It is worth mentioning that raising Wb reduces WF velocity, and the drop is of average scale.Positive values of Wb increased the Williamson fluid relaxation time.The higher relaxation time increased the viscosity of fluid due to which fluid flow faced higher resistance, and as a result, the WF declined.Appl.Sci.2021, 11, x FOR PEER REVIEW 10 o 0.6 on WF velocity was investigated, and observations in this way are shown in Figur

Figure 4 .
We established non-slip Williamson fluid flow w Vp = 0, and the WF velocity magnitude was larger than with Vp = 0.2 and 0.4.We disc ered that WF velocity over a magnetic surface decreased as Vp increased.Variations in = 0.1, 0.3, and 0.7 are shown in Figure 5 as a line graph of WF velocity.WF velocity has inverse relationship with positive Wb values, as shown in the graph.It is worth menti ing that raising Wb reduces WF velocity, and the drop is of average scale.Positive val of Wb increased the Williamson fluid relaxation time.The higher relaxation time increa the viscosity of fluid due to which fluid flow faced higher resistance, and as a result, WF declined.

Table 6 .Table 6 .
Nm variation towards Pr in magnetic and non-magnetic fields.Nm variation towards Pr in magnetic and non-magnetic fields.

Figure 11 .
Figure 11.Rs versus WF concentration.Tables 7 and 8 offer the Nusselt number outcomes for various parameters in different frames of reference.The Nm variations in thermal/non-thermal slip flows are reported in

Figure 11 .
Figure 11.Rs versus WF concentration.Tables 7 and 8 offer the Nusselt number outcomes for various parameters in different frames of reference.The Nm variations in thermal/non-thermal slip flows are reported in

Table 1 .
Skin-FC variation towards Wb in slip and non-slip flows.

Table 2 .
Skin-FC variation towards Pm in magnetic and non-magnetic fields.

Table 3 .
Skin-FC variation towards Wb in porous and non-porous mediums.

Table 4 .
Skin-FC variation towards Wb in magnetic and non-magnetic fields.
w F

Table 5 .
Skin -FC variation towards Mg in porous and non-porous mediums.

Table 7 .
Nm variation towards Pr in thermal and non-thermal flow fields.

Table 8 .
Nm variation towards Pr in porous and non-porous mediums.

Table 9 .
ShD variation towards Sc in porous and non-porous mediums.

Table 10 .
ShD variation towards Sc in magnetized and non-magnetized flow fields.

Table 11 .
ShD variation towards Sc in reactive and non-reactive flow fields.

•
Skin-FC at both non-permeable and permeable sheets showed a direct relation with Mg.
F Skin friction coefficient