Non-Newtonian Slippery Nanofluid Flow Due to a Stretching Sheet Through a Porous Medium with Heat Generation and Thermal Slip

The domains of engineering, electrical, and medicine all have a significant demand for nanofluids. Applications for nanofluid flow include electronic device storage, industrial cooling and heating frameworks, and associated medicinal management information systems. Nanofluids are utilized generally as coolants in heat exchangers such as thermostats, electronic cooling systems, and radiators due to their enhanced thermal characteristics. This study aims to explain the mixed convection phenomenon’s applications on the thermal impact of Maxwell nanofluid. The mass diffusivity is supposed to be a function of concentration, whereas the thermal conductivity and viscosity of Maxwell nanofluid are assumed to be functions of temperature. It is recommended to consider the additional thermal effects of thermal slip, magnetic fields, and heat generation phenomena. The fluid flow motion was caused by the vertically stretched sheet. The dimensionless formulation of the suggested physical model is shown by the suitable variables interacting. The shooting approach is used in the numerical simulations, and it is based on lowering higher-order nonlinear differential equations to first-order. The slip velocity and the magnetic parameters have a direct impact on the local skin friction coefficient and velocity, as indicated by the research findings. Also, the increase in values of the Maxwell parameter, porous parameter, and viscosity parameter leads to the enhancement of temperature distribution, while the decline in velocity distribution can be attributed to the same factors. A comparison is also made with the results described in the literature that is currently available, and a superb agreement is discovered.


Introduction
Many implementations are found in conducted over the past twenty years nanotechnology has advanced to become the leading technology, science, and biomechanics, chemical of the twenty-first century [1].Due to its many industrial, electronic, and medicinal uses, such as high-tech nuclear systems, solar collection, heat exchangers, engine cooling, and biological materials [2,3].The term of nanofluid was first introduced by Choi [4], he found that when added small nanoparticles in a base fluid increase the thermal conductivity.Earlier several researchers have carried out the studied the growth of nanofluid flow with heat and mass transfer under different physical conditions.Goyal and Bhargava [5] studied the non-newtonian nanofluid double diffusive boundary layer flow over a stretching sheet.Das et al. [6] investigated the heat and mass transfer of magnetohydrodynamics nanofluid flow over a vertical stretching sheet with non-uniform heat generation/absorption effect.Sulochana et al. [7] studied the effect of thermophoresis and Brownian motion on the heat and mass transmission over a stretching surface.Pal et al. [8] analyzed the Soret and Dufour effects on the mixed convection heat and mass transport of nanofluids over non-linear stretching and contracting sheets.The conjugate mixed convection heat and mass transfer with electrical magnetohydrodynamic and heat source/sink effects of nanofluid flow over a slip boundary stretching sheet surface was studied by Hsiao [9].Khan et al. [10] investigated numerically the changing viscosity and angled Lorentz force effects on Williamson nanofluid flow across a stretching sheet.Pal et al. [11] discussed the combined impacts of non-linear thermal radiation, non-uniform heat source, and sink on the heat transfer of Casson nanofluid flow over a stretching surface.The flow and heat transfer for unsteady nanofluid past a moving rotating plate with Brownian, and thermophoresis diffusion impacts was studied by Abbas and Magdy [12].The three-dimensional continuous thin film flow of tangent hyperbolic fluid past a stretching surface accompanied by nonlinear mixed convection flow and entropy formation was studied by Ibrahim and Gizewu [13].Loganathan et al. [14] investigated analytically the effects of heat-absorbing viscoelastic nanofluidic flow through heated porous Riga plate with Cattaneo-Christov double flux.Shao et al. [15] investigated the significance of visco-elastic time dependent materials and electrically conductive motion over gravitationally affected cylindrical surface.
Several researchers have accomplished various types of study the characteristics of flow and heat transfer over stretched surfaces due to its significance in numerous technological, industrial, and scientific applications, including the production of fibre, the condensation process, textile machines, and plastic sheets [16].Turkyilmazoglu [17] analyzed the heat and mass transfer characteristics of the magnetohydrodynamic viscous flow through a permeable stretched surface.The effects of velocity slip and heat generation/absorption on magnetohydrodynamic stagnation point flow and heat transfer over a contracting/expanding surface was investigated by Nandy et al. [18].Gireesha et al. [19] examined the impact of nanoparticles on a steady magnetohydrodynamic of Eyring-Powell

3
Journal of Nonlinear Mathematical Physics (2023) 30:1221-1238 fluid flow with heat and mass transfer over a stretching sheet.The impacts of viscous dissipation, Joule heating and convective boundary condition of magnetohydrodynamics micropolar liquid flow past nonlinear stretched surface were presented by Waqas et al. [20].Jena et al. [21] explored the Soret and Dufour effects on magnetohydrodynamics viscoelastic nanofluid flow over a porous vertical stretching sheet.The boundary layer viscous flow of nanofluids and heat transfer over a non-linearly stretching sheet in the presence of a magnetic field was demonstrated by Ramya et al. [22].Mabood et al. [23] investigated numerically the thermal radiation, suction/injection, and multiple slip effects on MHD unsteady nanofluid flow with heat and mass transfer over a stretching sheet.The effect of the aligned magnetic field on Williamson's nanofluid flow on a stretching surface with convective boundary conditions was investigated numerically by Srinivasulu et al. [24].Megahed [25] addressed the steady MHD of non-Newtonian Maxwell fluid flow over an imbedded stretched sheet in a porous medium with convective boundary condition.The MHD free convective stagnation-point flow toward an inclined nonlinearly stretching sheet embedded in a porous medium was presented by Biswal et al. [26].Alali et al. [27] investigated numerically the thin film flow and heat mass transfer for Casson nanofluid on a horizontal elastic sheet in presence of viscous dissipation, thermophoresis, and the slip velocity effects.The flow of a non-Newtonian Williamson fluid over stretched sheet embedded in a porous medium and combining with the Cattaneo-Christov model was studied numerically by Abbas et al. [28].The literature review reveals several studies that discuss the feasibility of obtaining exact solutions for certain nonlinear models, as demonstrated in works by Yu-Hang Yin et al. [29] and Si-Jia Chen et al. [30].However, our proposed model exhibits high nonlinearity, necessitating the utilization of numerical methods for solution determination.
Based on the research above and motivated by the possible technological applications regarding the nanofluid topic, the novelty of the proposed study can be emphasized by examining the combined effects of mixed convection and viscous dissipation phenomena on the flow of a Maxwell nanofluid over a stretching sheet embedded in a porous medium, considering the presence of slip velocity.The unique contribution and significance of this study lie in its unprecedented treatment of Maxwell nanofluids, offering novel insights that have not been explored previously.Further, the motivation for this study stems from several factors.Firstly, Maxwell fluids exhibit distinct rheological properties that differ from conventional Newtonian fluids.Gaining an understanding of the flow characteristics and heat transfer properties of these non-Newtonian fluids is essential for diverse industrial and engineering applications.Secondly, the inclusion of nanoparticles in slippery nanofluids offers the potential for significant enhancements in heat transfer performance.Exploring the flow and heat transfer behavior of these nanofluids is highly valuable in improving energy efficiency and thermal management across various systems.The findings reported here, according to the authors, have not been previously published and may fill in any gaps.

Formulation and Description for the Physical Model
Let's examine the two-dimensional, steady, mixed convective boundary layer flow caused by a vertically stretched sheet under the influence of non-uniform heat generation, as shown in Fig. 1, where the sheet velocity is U w = ax and a is a constant.
The assumption is that the nanoparticles are homogenous in size and shape.Choose a cartesian coordinate system where the x−axis is generally parallel to the sheet and the y−axis is orthogonal to the sheet.The sheet is also considered to be immersed in a porous medium.With regard to the fact that the current work adopts Buongiorno's nanofluid model, which incorporates thermophoresis and Brownian motion impacts.It is assumed that there is a magnetic field with a strength of B 0 in the y−direction.Due to the low Reynolds number, the induced magnetic field is disregarded.
In this research, we assume that the suggested model is subject to an internal heat generation or absorption phenomenon that can be expressed by the term q ′′′ and may be mathematically modelled as follows [31]: where ∞ is the kinematic viscosity of the nanofluid away from the sheet, T w is the sheet temperature, T ∞ is the ambient temperature.Also, in our model, T w is assumed to correlate with T ∞ according to the relation T w = T ∞ + Ax , where A is a constant.Likewise, is the nanofluid thermal conductivity which obey the following relation [25]: (1) where ∞ is the ambient thermal conductivity and 1 is the thermal conductivity parameter.Additionally, the following relationship is considered to affect the Maxwell nanofluid's viscosity [25]: where is the viscosity parameter and ∞ is the ambient viscosity of the nanofluid.Also, the base fluid's nanoparticle concentration is denoted as C w .Additionally, as y approaches infinity C ∞ exhibits a concentration of nanoparticles.Further, we assume that the motion of the nanoparticles in the base fluid is caused according to the Brownian diffusion coefficient D B which obey the following relation [32]: where 2 is the mass diffusivity parameter and d 1 is the reference mass diffusivity.Moreover, the movement of nanofluids exhibits a thermophoretic behavior with a thermophoretic diffusion coefficient D T , where the reactivity of the nanoparticles to temperature gradients is distinct.In fact, the suitable partial differential equations governing the aforementioned hypothesized MHD Maxwell nanofluid flow model are provided in terms of the following equations for mass, momentum, energy, and nanoparticles' molar concentrations: where u and v, respectively, stand for the velocity elements in the x and y directions.Maxwell's parameter is 1 , the density of the ambient nanofluid is ∞ , the acceleration due to gravity is g, T is the nanofluid's temperature, C is the concentration expansion coefficient, K is the permeability of the porous medium and the coefficient of thermal expansion is T .Further,the boundary conditions for the temperature, velocity, and concentration distributions are as follows: ( where a, c are constants, h T is the coefficient of heat transfer and 0 is the factor of the slip velocity.Now, we begin the dimensionless variables in the formats listed below before creating the solution algorithm: where is the non-dimensional fluid temperature, f is the non-dimensional stream function, is the dimensionless concentration and is the dimensionless similarity variable.We now see that the continuity equation Eq. ( 5) is quickly and easily satisfied once these transformations (11)(12) have been applied.Additionally, employing the same non-dimensional transformations that were suggested earlier, we can obtain the following governing ordinary differential equations: Additionally, with equivalent boundary conditions, the following is stated: Journal of Nonlinear Mathematical Physics (2023) 30:1221-1238 where M = B 2 0 ∞ a is the magnetic parameter, = ∞ aK is the porous parameter, = a 1 is the Maxwell parameter, is the thermal slip parameter, is the Lewis parameter.Here in Eq. ( 16), it is important to note that the condition will physically represent the thermal slip phenomena, which is independent of thermal conductivity, if 1 = 0 .Additionally, when = 0 , our non-Newtonian nanofluid model can be transformed into a Newtonian one.The following are other significant physical factors for the drag force Cf x , heat transfer rate or the local Nusselt number Nu x , and the rate of mass transfer or the local Sherwood number Sh x : where is the local Reynolds number.

Code Validation
Comparisons with the findings reported in the literature by Hayat et al. [33] for the Newtonian fluid situation for various values of slip velocity parameter 1 when = = = M = = 0 are done in order to confirm the numerical results.
Here, it is important to note that Hayat et al. [33] used a homotopy analysis method to examine the problem of viscous liquid flow past a stretching sheet under conditions of maintaining velocity and thermal slip as well as the presence of thermal radiation effects.Table 1 provides this comparison.This table demonstrates the ( 18) great concordance between our findings and those reported in the literature.So we may certainly say that the data provided here is authentic.

Discussion of the Results
Using graphical and tabular representations, relevant numerical results are presented in this section.Numerous physical parameter values, including the magnetic parameter M, the porous parameter , the Maxwell parameter , the slip velocity parameter 1 , the mixed convection parameter , the thermal slip parameter 2 , the thermophoresis parameter Nt, the Lewis parameter Le and the Brownian motion parameter Nb, were computed.Figure 2 depicts how altering the Maxwell parameter affects the nanofluid flow's velocity f � ( ) , temperature ( ) and concentration ( ) fields.The velocity profile clearly exhibits a magnifying effect caused by the Maxwell parameter, which is in opposition to the temperature and concentration fields.In physical terms, a higher value of the Maxwell parameter indicates an increased resistance affecting the flow behavior.This heightened resistance exerts a stronger drag force on the fluid particles, leading to a decrease in velocity.
Figure 3 illustrates the relationship between the porous parameter and the temperature ( ) , velocity f � ( ) , and concentration ( ) of the nanofluid.Because of the increased flow impedance caused by the presence of a porous medium, the flow velocity falls and the temperature and concentration of solutes rise.Furthermore, from a physical interpretation standpoint, this finding aligns with the observation that reducing the permeability of the porous medium decreases the drag force.As a consequence, this reduction in drag force leads to a decrease in the velocity of the flow.
Figure 4 shows how M varies with f � ( ) , ( ) and ( ) .As M grows, the tempera- ture and concentration profiles both rise, however the velocity behavior shows the opposite tendency.The Lorentz force causes the frictional drag to increase, which causes the temperature and concentration boundary layer thickness to boost and the velocity field to diminish.A magnetic parameter indicates a stronger magnetic field's impact on the flow, leading to an increased Lorentz force that hampers the velocity.Additionally, the magnetic field promotes heat and mass transfer processes, resulting in elevated temperature and concentration.The magnetic forces induce enhanced mixing and facilitate the transportation of thermal energy and solute particles, consequently raising the temperature and concentration levels.
The velocity f � ( ) of the nanofluid, the temperature ( ) and the concentration of the nanoparticles ( ) are plotted against various viscosity parameter in Fig. 5.While a large values of the viscosity parameter causes the concentration and temperature fields to be amplified and the velocity field to be reduced, it is also seen that the parameter increases the thickness of the thermal boundary.Further, higher viscosity results in less heat being transferred to the surroundings, warming the flow and increasing both temperature and the concentration of nanoparticles.
As the slip velocity parameter 1 was increased, Fig. 6 showed how the motion of nanofluids f � ( ) , temperature distribution ( ) , and the concentration ( ) of nanoparticles changed.This graph displayed the positive slip parameter values, which cause the nanofluids' concentration and temperature to increase while slowing the flow of the fluid.Physically, a high slip velocity parameter signifies a substantial slip between the and solid surface, resulting in increased momentum and mass transport.This heightened transport fosters improved heat and mass transfer processes.The amplified slip velocity encourages greater mixing and interaction among fluid particles, facilitating the movement of thermal energy and solute particles.Consequently, the temperature and concentration levels escalate within the flow region, indicative of the intensified heat and mass transfer arising from the larger slip velocity parameter.
In Figure 7, while keeping the other parameters unchanged, the influence of the mixed convection parameter is examined in relation to the distribution of velocity f � ( ) , temperature ( ) , and nanoparticle concentration ( ) .It is obvi- ous that the mixed convection parameter only appears in the equation for the momentum boundary layer.As a result, it has a direct impact on the velocity field and a secondary impact on the temperature and concentration fields.In contrast to the temperature and concentration fields, which show a tendency in the opposite direction as the mixed convection parameter is increased, the velocity distribution is found to increase.
In Figure 8a, the impact of the thermal slip parameter 2 on the temperature distribution ( ) is examined while maintaining the other parameters at their constant values.It is discovered that as the thermal slip parameter's increases up, both the sheet temperature (0) and the nanofluid temperature within the Journal of Nonlinear Mathematical Physics (2023) 30:1221-1238 boundary layer decrease.As a result, the thermal boundary layer will be thinner and the rate of heat transmission will be increased.The temperature field ( ) results are shown in Fig. 8b for the thermal conductivity parameter 1 .Due to the increased thermal conductivity of solid particles compared to base fluid, the volume of nanoparticles enhances the overall thermal conductivity of nanofluids.The thickening of the thermal boundary layer is caused by an increase in thermal conductivity.Therefore, the temperature as well as the thickness of the boundary layer grow as the overall thermal conductivity of nanofluids improves.
The performance of the temperature profile ( ) in response to the space dependent heat source parameter A * is depicted in Fig. 9a.The findings demon- strated that as the space based heat source parameter was altered, both the temperature distributions ( ) and the sheet temperature (0) significantly increased.Further, Fig. 9b showed the Maxwell nanofluid temperature distribution characteristics as the thermally heat source parameter B * was increased.It has been found that raising the values of the temperature-dependent heat source parameter raises both the temperature of the sheet (0) and the temperature distribution ( ) , which has an significant impact on the thermal thickness.
Figure 10a show how the thermophoresis parameter Nt impacts the thermal field ( ) .The temperature and concentration boundary layer equations can be observed to contain the thermophoresis parameter.As we observe, it is connected with the temperature field and significantly affects the boundary layer's heat transport and nanoparticle concentration.The increase in the thermophoresis parameter thus leads to an improvement in the thermal boundary layer thickness and temperature distributions for the flow of Maxwell nanofluids.Also, the properties of the thermal thickness of the boundary layer and the temperature distribution ( ) are affected by the Brownian motion parameter Nb, as shown in Fig. 10b.The presence of nanoparticles in nanofluid media enhances heat transmission capabilities through Brownian movement, which promotes heat conduction and increases the heat transfer surface area.The temperature distribution in the nanofluid flow, represented by ( ) , is influenced by both the Brownian motion parameter and the sheet temperature, (0) .The Brownian motion parameter affects the dispersion and movement of nanoparticles, while the sheet temperature determines the initial conditions for the temperature distribution.The concentration patterns ( ) for various values of the mass diffusivity param- eter 2 are shown in Fig. 11a.When the mass diffusivity parameter is set to a significant value, the amount of nanoparticles present increases in concentration.Additionally, Fig. 11b depicts the behavior of the concentration distribution ( ) for the investigated range of Lewis number Le.In addition to describing the diffusivity relationship between heat and mass, Lewis number is a ratio of Schmidt number and Prandtl number.
It makes sense that when Le rises, kinematic viscosity falls to facilitate quicker mass transfer to the surroundings.The nanoparticle concentration decreases as you move farther away from the sheet, as seen by the flows in Fig. 11b.
When it comes to engineering and application, the forms of the velocity, temperature, and concentration profiles are frequently less important than the skin-friction coefficient  2 shows that an increase in the Maxwell parameter, the porous parameter, the magnetic number, the viscosity parameter, and the velocity slip parameter causes a reduction in the Nusselt number or the rate of heat transfer at the boundary.Shear stress and the local Nusselt number grow along with the thermal slip parameter values, but they diminish when the slip velocity parameter values are elevated.The local Sherwood number is shown to have a falling function for the Maxwell parameter, porous parameter, magnetic number, viscosity parameter, and slip velocity parameter in the same table, but the opposite is seen for the mixed convection parameter and Browninan motion parameter.

Conclusive Remarks
To examine the impact of physical factors on the heat and mass transfer of a non-Newtonian Maxwell slippery nanofluid flow on a vertical stretched sheet, a mathematical and theoretical investigation was carried out.The stretching sheet is thought to be submerged in a porous medium.Studies are also done on the importance of Brownian motion, thermal slip, velocity, and non-uniform heat generation.In this study, the ability to vary thermal conductivity and diffusivity within the nanofluid flow is incorporated.The primary aim of this study is to investigate the influence of a magnetic field on the flow characteristics of a Maxwell nanofluid within a porous medium.Ordinary differential equations with 1 3 Journal of Nonlinear Mathematical Physics (2023) 30:1221-1238 suitable boundary conditions are created by simplifying the set of governing equations and the boundary condition using dimensionless transformations.In order to treat modelled equations numerically, the shooting method is used.The key findings of this contemporary research are listed below: 1.The wall heat transfer coefficient is shown to be significantly affected by the thermal slip parameter and the mixed convection parameter.2. The velocity, temperature, and concentration profiles are found to be significantly affected by the influences of the slip velocity parameter and magnetic parameter.3. When a porous media is available and a non-Newtonian Maxwell fluid is provided, the velocity is lower than it is when a Newtonian fluid is present and a porous medium is absent.4. It has been discovered that the parameters that control slip velocity, viscosity, and magnetic field result in a thicker thermal and concentration boundary layer.5. Skin-friction coefficients, the local Nusselt number, and the local Sherwood number are all affected by slip velocity parameter variations in a similar manner.6.In future studies, there is potential to expand this research by incorporating the Jeffery nanofluid and investigating the impact of thermal and concentration slips.

Fig. 1 3
Fig. 1 Sketch of the physical fluid flow problem