Non-Newtonian nanofluid flow across an exponentially stretching sheet with viscous dissipation: numerical study using an SCM based on Appell–Changhee polynomials

The objective of this article is to investigate how the properties of a non-Newtonian Williamson nanofluid flow, which occurs due to an exponential stretching sheet placed in a porous medium, are influenced by heat generation, viscous dissipation, and magnetic field. This study focuses on analyzing the heat transfer process by considering the impact of temperature on the thermal conductivity and viscosity of Williamson nanofluids. Additionally, the research significantly contributes by investigating the flow characteristics of these nanofluids when influenced by slip velocity. Using the spectral collocation method (SCM), the equations that describe the current problem are transformed into a collection of ordinary differential equations and then solved. The SCM proposed here basically depends on the properties of the Appell-type Changhee polynomials (ACPs). First, with the aid of ACPs, we give an approximate formula of the derivatives for the approximated functions. Through this procedure, the provided model is transformed into a nonlinear set of algebraic equations. Physical factors of interest, such as skin friction, the Nusselt number, and the Sherwood number, are explained using tabular expressions. Data are displayed as graphs for the nanofluid’s velocity, temperature, and concentration. The primary findings showed that increasing the Williamson, magnetic, thermal conductivity, and Brownian parameters significantly improves the thermal field. Finally, testing the suggested method with specific cases from some past literature-based publications reveal a good degree of agreement.


Introduction
Due to the limited applications of Newtonian fluid models, non-Newtonian fluid research has drawn interest.Non-Newtonian fluids include materials like starch, honey, ketchup, lubricating sprays, and more.Earlier, many non-Newtonian fluid models were proposed.From these a Williamson model was constructed.In 1929, Williamson [1] proposed a non-Newtonian fluid model that illustrates the rheological characteristics of these fluids.Shear thinning is a characteristic of non-Newtonian fluids like in the Williamson fluid model.In the literature, the Williamson fluid model is the description of a fluid that thins under shear.The Williamson model has been utilized by numerous researchers (see [2][3][4][5][6][7][8]) due to its applicability in illuminating fluid dynamics in the preceding decade.
Due to its applicability in the industrial sector, nanotechnology is currently attracting the attention of scientists and researchers, for instance, the application of nanofluids to cool nuclear reactors, regulate heat flow through heat valves, lower automobile radiator temperatures, cool computer processors, etc. Cancer patients are treated in the medical field using medications and radiation delivered by devices made of iron-based nanofluids.By incorporating nanoscale material particles in regular fluids, Choi [9] established the name "nanofluid" and mathematically validated the approach.The nanoparticles mainly consist of metals, nitride ceramics, oxide ceramics, carbide ceramics, and common base fluids like methanol, water, ethylene glycol, and oil.The results of experiments by the authors of [10,11] proved that nanofluids have a higher thermal conductivity than regular fluids.Mixed convection is a phenomenon where fluid flow involves both natural and forced convection mechanisms simultaneously [12].It holds significant importance in practical applications, especially in the fields of heat transfer and fluid dynamics.Its significance lies in its capability to affect fluid flow patterns, boost heat transfer rates, enhance energy efficiency, and influence thermal stratification.Precise modeling and analysis of mixed convection phenomena are crucial for optimizing the efficiency and performance of various engineering systems in diverse industries.Due to the significance of mixed convection, Abbas et al. [13] addressed this phenomenon in their study on nanofluid flow in the second grade.After that many authors (see [14][15][16][17][18]) have presented an extensive survey of the flow of nanofluids in various thermo-physical contexts.
In this paper, we will, for the first time, derive an approximate formula of the derivatives with the help of the Appell-type Changhee polynomials and apply it to solve the model under study by using the spectral collocation approaches (see [19,20]).The most famous advantage of these methods is their capability to generate accurate outcomes with a very small degree of freedom error [21].They are widely used because of their good properties in the approximation of functions.The orthogonality property of the ACPs is used to approximate functions on their domain.These polynomials have a main and important role in the methods for ODEs [22].Many researchers used and implemented these polynomials to numerically solve many problems, such as in [23] were they were used to solve the high-dimensional chaotic Lorenz system.Therefore, the primary goal of the current paper is to examine the behavior of a mixed convection boundary layer flow of a non-Newtonian nanofluid flow that includes nanoparticles towards a vertical stretching permeable sheet under slip velocity, nonuniform heat generation, and magnetic field impacts.This goal is motivated by the earlier cited references on the subject of nanofluid models.For the solution of the coupled nonlinear model equations, the spectral collocation method is applied.The motivation of this investigation is to learn more about the non-Newtonian Williamson nanofluid problem and to improve our understanding of it.By examining this research, we hope to add to the body of knowledge, fill in knowledge gaps, and perhaps open the door to breakthroughs in a variety of disciplines, including engineering, materials science, and nanotechnology.The primary objective of this study is to investigate the influence of slip velocity, magnetic field, internal heat generation, and viscous dissipation on the flow behavior of Williamson nanofluid, providing a distinctive viewpoint.The study's novelty and purpose stem from the fact that it is the first of its type to implement the proposed numerical technique to solve the proposed model.The proposed method possesses numerous advantages.The use of Appell-Changhee polynomials in the spectral collocation method results in significant advantages such as exceptional precision, extensive usability, global approximation, efficient execution, high convergence order, and numerical reliability.These advantageous characteristics establish it as a potent and dependable numerical approach suitable for tackling diverse mathematical problems.
This study is organized as follows.A full description of the problem is given in Sect. 2. In Sect.3, we give some basic concepts about Appell-Changhee polynomials, and an approximation of the D n via ACPs, with procedure solution using SCM.The code validation is given in Sect. 4. In Sect.5, we present the results and discussion.Finally, the conclusions are given at the end of the paper through Sect.6.

Basic model
Here, we consider the two-dimensional non-Newtonian Williamson flow of a nanofluid along an exponentially stretched sheet with a constant temperature T w and concentration C w .The appropriate values for ambient concentration and temperature are C ∞ and T ∞ , respectively.Also, we presume that the constant axial surface temperature T w is greater than the ambient fluid temperature T ∞ .Based on the flow model, the behavior of thermophoresis and Brownian motion are also taken into consideration.Additionally considered are the novel slip velocity, nonuniform heat generation with Brownian motion, and thermophoresis properties.The velocity is assumed to be U w = U 0 e x L and the x-axis is taken in the direction of the stretched sheet, where L is a characteristic length and U 0 is a constant velocity (Fig. 1).Furthermore, a suction velocity of v w is produced by the hypothesis that the sheet is porous.
In addition, the following fluid properties are taken into consideration as a variable besides the assumption of nonuniform heat generation.
(i) Viscosity We assume that the nonlinear exponential function of temperature and the nanofluid viscosity μ are related to one another [24], i.e., μ = μ ∞ e -α( T-T∞ Tw-T∞ ) , where μ ∞ is the nanofluid viscosity away from the sheet and α is the viscosity parameter.Also, this relation explains that the nanofluid viscosity behavior is influenced by the thermal parameter α.Usually, α is regarded as a positive value for a liquid.
(ii) Thermal conductivity Thermal conductivity κ(T) demonstrates how a material's conductivity changes.Here, we assume that the relationship between the nanofluid thermal conductivity and the distributed temperature is as follows [24]: where ε is the thermal conductivity parameter, assumed to be small, and κ ∞ is the ambient thermal conductivity.Further, if ε = 0, this relationship indicates that thermal conductivity is a constant.
(iii) Nonuniform heat generation One can observe that the nonuniform heat generation q * is directly articulated as follows in light of the multiple correlations between nonuniform heat generation and the sheet as well as the ambient temperatures [25]: where D 1 is the kinematic viscosity, ν ∞ and E 1 are constants.The conservative equations after using the boundary layer and Boussinesq estimations yield under the consideration of all the aforementioned factors [26] and the related boundary conditions are [8]: It is crucial to note that the momentum equation ( 2) introduces two significant additions: the inclusion of a variable thermal viscosity and the mixed convection term.Similarly, the energy equation ( 3) incorporates the viscous dissipation term and the internal heat generation term.Also, it is important to emphasize that the initial component of equation ( 5) represents the phenomenon of slip velocity, where the x-and y-axis coefficients of the fluid velocity vectors are u and v, respectively; D B , , B 0 represent diffusivity of the Brownian motion, the Williamson parameter, and the magnetic field intensity; τ is the ratio of the nanoparticle's effective heat capacitance to that of the base fluid, V 0 is the constant suction velocity, k 0 is the porous medium's permeability, D T is the thermophoretic diffusion coefficient, σ is electrical conductivity, β c is the coefficient of concentration expansion, β T is the temperature expansion factor, g is the acceleration due to gravity, c p is the specific heat at constant pressure, λ 1 is the factor of slip velocity, and ρ ∞ is the ambient density.Now, the following dimensionless transformation is used to convert the nonlinear partial differential equations ( 1)-( 4) into nonlinear ordinary differential equations (ODEs), assuming the apposite similarity variable as shown below [26]: where f is the dimensionless stream function, θ (η) is the dimensionless temperature, and φ(η) is the dimensionless concentration.Now, equation ( 7) enables the continuity equation (1) to be satisfied precisely.The following dimensionless system of ODEs with associated boundary constraints is obtained from the system (2)-( 4) with boundary constraints ( 5)-( 6) by invoking transformation ( 7)- (8).So, we have subject to the following restrictions on boundaries: where is the local modified Grashof number, is the thermophoresis parameter.Accordingly, the Sherwood number, Nusselt number, and skin friction are [27]: where 3 Solution procedure

Basic concepts about Appell-Changhee polynomials
It is well-known that Changhee polynomials Ch m (t) are usually defined using generating functions (see [19,22]), where Ch m = Ch m (0) are the Changhee numbers, see [19].These polynomials can also be given in the following form: where S 1 (m, ) and E (t) are Sterling numbers of the first kind and Euler polynomials, respectively.However, the Appell-type Changhee polynomials Ch * m (t) are defined by the generating function given by [28] The ACPs of degree m are defined by From formula (14), one can easily get that therefore from (15), we can confirm the following formula: It is also worth noting that Ch * 0 = 1 and 2Ch * m + mCh * m-1 = 0, ∀m ≥ 1.
In addition, we can prove that the ACPs satisfy the following identity: Let {Ch * i (t)} m i=1 ⊂ L 2 [0, 1] be the set of ACPs and suppose that For a function g(t) of L 2 [0, 1], one has a good and unique approximation of it in .If g * (t) is the unique approximation of g(t), we can write the following error estimate: But because is a closed subspace of L 2 [0, 1], according to [29], we can write L 2 [0, 1] = ⊕ ⊥ , where ⊥ denotes the orthogonal complement of , and so we have and then r(t) = g(t)h(t), which also means that g(t)g * (t) ∈ ⊥ .Therefore, this confirms the following: where •, • denotes the inner product.Since g * (t) ∈ , we can write the following: where By taking h(t) = Ch * i (t) and substituting equation ( 18) into (17), we get Also from (18), we can see that where A = Ch * (t), Ch * (t) is an N × N matrix, defined by and A can be computed by using (16).Therefore, the coefficients' vector from (19) can take the form In this subsection, we will show that the nth derivative D n of the function g * (η) given in (18) can be approximated through the following theorem.

Theorem 1
The derivative of order n > 0 for the function g * (η) given in (18) can be approximated by where Ch * i-j is the Changhee number, and χ i,j,n is given by Proof Considering the ACP, Ch * i (η) of degree i, with i = 0, 1, . . ., N , and by using (18), we can get where χ i,j,n is given in (20), and this completes the proof.

Code validation
In this section, we compare the numerical values of the skin friction coefficient against various values of the Williamson parameter W e with the data calculated by Nadeem and Hussain [30] to confirm our numerical results that were achieved using the spectral collocation method, which depends on the characteristics of the ACPs.As shown in Table 1, these findings are closely linked to those reported in the literature by Nadeem and Hussain [30].Therefore, we can confidently state that the data presented here is accurate.

Results and discussion
The findings of using the SCM to establish the significant features of the flow, heat, and mass transfer characteristics with nanoparticles are shown in Figs. 2 through 13 and in Table 2. Figure 2 depicts how a magnetic parameter M affects temperature θ (η), velocity f (η), and concentration φ(η) graphs, respectively.It can be shown in Fig. 2 that the Lorentz force becomes strong for large magnetic field parameters, increasing the fluid resistance and causing velocity to decrease.Additionally, it is observed that as the magnetic parameter increases, species concentration drops, followed by a decrease in the thickness  of the solutal boundary layer.Furthermore, it appears that a rise in the magnetic parameter improves the temperature distribution.Figure 3 displays the variations in the fields of velocity f (η), temperature θ (η), and concentration φ(η) for various values of the slip velocity parameter λ.Diminishing behavior is observed in the velocity field with rising values of the slip velocity parameter λ, whereas the temperature field exhibits the opposite trend for the same parameter.It is shown that increasing the slip velocity parameter decreases the volume friction of nanoparticles, improving the boundary layer thickness and concentration profile.Physically, the reduction in nanofluid velocity associated with the slip velocity parameter can be explained by slip flow, which refers to the situation where the fluid near a solid surface moves at a different speed compared to the surface itself.When nanoparticles are present in nanofluids near the solid surface, they form a boundary layer that influences the flow behavior of the fluid.The slip velocity parameter quantifies slip flow and represents the relative velocity between the fluid and the solid surface.
For a particular viscosity parameter α values that define the flow behavior and the heat mass transfer spread through the boundary layer, Fig. 4 illustrates the velocity, temperature, and concentration curves.It was made clear that an increase in the viscosity parameter led to a greater improvement in the panoplied's viscosity, which was related to the fact that when the viscosity parameter increases, the fluid velocity decreases.Additionally, a bigger viscosity parameter leads to stronger mass and heat transport, which enhances the temperature and concentration fields.Physically, the viscosity parameter gauges a fluid's ability to flow.Nanoparticles alter the viscosity of nanofluids when they are added.As a result, the velocity of the nanofluid decreases.The addition of nanoparticles alters the interactions between the fluid molecules, increasing the viscosity of the nanofluid as a whole.This increased viscosity causes the fluid to experience greater internal friction, which slows the flow and lowers the velocity.
Figure 5 examines the impact of the Williamson parameter W e on the distributions of f (η), θ (η), and φ(η).In contrast to concentration and temperature distribution, the velocity distribution changes.The temperature and concentration distribution advance as the Williamson parameter increases.It indicates that as nanoparticle mobility increases, the thermal conductivity of nanofluids also increases, which boosts the heat-and masstransfer mechanisms.
The velocity, concentration, and temperature profiles for several porous parameter γ values are shown in Fig. 6.In this case, raising the porous parameter γ induces the velocity profile to drastically deteriorate, although the temperature and concentration fields exhibit the opposite tendency.Additionally, the porous parameter produces a small improvement in the thickness of the thermal boundary layer.
Figure 7 shows how raising the local Grashof number Gr T improves the velocity profile f (η).Physical justification can be found in the fact that as the local Grashof number increases, the convection mechanism improves, increasing the nanofluid velocity.In contrast, the temperature θ (η) and concentration φ(η) fields exhibit different behavior.Physically, the Grashof number is a dimensionless parameter that captures the equilibrium between buoyant and viscous forces in fluid flow.The following physical explanation ex- plains why a higher Grashof number could result in a faster nanofluid: An increase in the Grashof number denotes a greater dominance of buoyant forces over viscous forces.The fluid consequently encounters a stronger pushing force as a result of temperature or density variations.This more powerful pushing force causes the nanofluid to move more quickly, which improves convective heat transfer.As a result, an increase in the Grashof number causes the nanofluid velocity to increase, encouraging more active fluid flow and better convective heat transfer.
Plotted in Figs.8(a) and 8(b) are the findings of the velocity, temperature, and concentration distributions for the local modified Grashof number Gr C operating in the laminar Williamson nanofluid flow zone.According to the findings in these graphs, the velocity field and boundary layer thickness both improve as the locally modified Grashof number increases.Additionally, the opposite is shown for the same parameter on mass and heat transport.
The relationship between the thermal conductivity parameter ε and the Williamson nanofluid flow and heat mass properties is depicted in Fig. 9.This graph demonstrates that the velocity f (η), concentration φ(η), and temperature θ (η) profiles through the boundary layer region are ascending in nature for larger values of the thermal conductivity parameter.This is because the conductivity behavior of the nanofluid advances as the thermal conductivity parameter increases.Physically, when it comes to nanofluids, the presence of nanoparticles can increase the fluid's thermal conductivity parameter, which is a measure of a material's capacity to conduct heat.And so, the nanofluid exhibits higher tempera- For various values of the thermophoresis parameter Nt, Fig. 10 depicts the variations in the concentration and temperature fields.It has been noted that as the thermophoresis parameter is elevated, both concentration and temperature profiles rise because thermophoresis is a phenomenon by which small particles are drawn from hot to cool surfaces, which causes the temperature to rise.
Figure 11 shows how the temperature θ (η) and concentration φ(η) fields for nanofluids are affected by the Brownian motion parameter Nb.This figure illustrates how raising the Brownian motion parameter causes fluid particles to move randomly, which raises the concentration and temperature profiles.
The temperature profile θ (η) within the thermal boundary layer region is affected by Eckert number Ec, as shown in Fig. 12(a), since an increase in the Eckert number speeds up the movement of kinetic energy.Fluid particles collide more regularly as a consequence, which causes kinetic energy to be converted into heat energy in the process.As a result, the temperature profile rises, as seen in the graph.Figure 12(b) shows how the spacedependent heat source parameter D 1 influences the temperature profile.The fluid temperature increases for high values of the space-dependent heat source parameter; this phenomenon is physically maintained because high values of D 1 result in increased heat generation and increased temperature distribution for a nanofluid.The effect of the temperature-dependent heat source parameter E 1 on the temperature and concentration profile is seen in Fig. 13.The graph demonstrates how increasing the temperature-dependent heat source parameter increased the capacity of heat and mass transference, leading to an increase in both the nanofluid temperature and concentration.
Table 2 shows the variations in the skin friction coefficient √ 2Cf x Re

Conclusions
In this work, non-Newtonian Williamson nanofluid flow characteristics caused by a vertical exponential stretching sheet were addressed.The spectral collocation technique which basically depends on the properties of the Appell-type Changhee polynomials was used to give a numerical solution to the aforementioned problem.The flow was thought to be in a steady state.Williamson nanofluid flows' slip velocity phenomena, the influence of magnetic fields, heat generation, and viscous dissipation all play a role in how this study is observed.Graphical analysis was used to view and analyze the impact of the embedded factors.These numerical calculations and mathematical modeling have essential applications, particularly in cooling operations.As a consequence, the present study can specify several essential aspects that are crucial in engineering applications.The key findings are as follows: 1.The skin friction decreases as the parameters for viscosity, slip velocity, and the local Grashof number are improved.2. In contrast to the effects of the magnetic field, slip velocity, and viscosity parameters on fluid motion, an increase in the local Grashof number causes fluid velocity to accelerate.3. The Nusselt number is increased by the slip velocity and local modified Grashof numbers, whereas it is decreased by the magnetic, Williamson, and porous parameters.4. The temperature of nanofluids is controlled by the magnetic field and higher levels of the local Grashof number.5.By increasing the values of the thermophoresis parameter, the Nusselt number declines, whereas the viscosity parameter exhibits the reverse tendency.6. Increased values of the variable viscosity parameter, Grashof number, modified Grashof number, and thermal conductivity parameter result in a decrease in the skin friction coefficient.7.In future endeavors, it would be advantageous to broaden the scope of analysis conducted in this study by examining the impact of different variables such as nanoparticle concentration, porosity, and the chemical reaction on the flow behavior of non-Newtonian Williamson nanofluids.

Figure 1
Figure 1 Schematic geometry of the flow problem

Table 1
Values of skin friction√ 2Cf x Re 1/2x with various W e and the results of Nadeem and Hussain 2L Sh x caused by modifying the magnetic field parameter, slip velocity parameter, local Grashof number, thermal conductivity parameter, thermophoresis parameter, porous parameter, local modified Grashof number, and Williamson parameter.It is evident that with high values of the thermal conductivity,