In the current research, the numerical solutions for heat transfer in an Eyring–Powell fluid that conducts electricity past an exponentially growing sheet with chemical reactions are examined. As the sheet is stretched in the direction, the flow occupies the region . MHD, radiation, joule heating effects, and thermal relaxation time are all used to represent the flow scenario. The emergent problem is represented using PDEs, which are then converted to ODEs using appropriate similarity transformations. The converted problem is solved numerically using the SLM method. The main goal of this paper is to compare the results of solving the velocity and temperature equations in the presence of changes through SLM, introducing it as a precise and appropriate method for solving nonlinear differential equations. Tables with the numerical results are created for comparison. This contrast is important because it shows how precisely the successive linearization method can resolve a set of nonlinear differential equations. Following that, the generated solution is studied and explained in relation to a variety of engineering parameters. Additionally, the thermal relaxation period is inversely proportional to the thickness of the thermal boundary layer and the temperature, but the Eckert number is the opposite. As grows, the temperature within the channel increases.

1. Introduction

Non-Newtonian fluids are widely encountered and are used in a wide variety of engineering applications. Some of these applications are notable and are applied in the paper, food, personal care products, textile coating, and suspension solutions industries. These fluids have mostly been divided into three categories: differential, rate, and integrals. Recent technological and engineering advancements have resulted in the development of a diverse range of non-Newtonian fluids with a number of major differences from viscous fluids. Ziegenhagen [1] explored the slow flow of a Powell–Eyring type fluid and used variation techniques to obtain results. He studied the behavior of Oldroyd and Powell–Eyring fluids and discovered that both fluids behave identically in situations involving extremely slow fluid flow. Sirohi et al. [2] studied it by observing the flow of Powell–Eyring fluid around the accelerating plate. They compared three distinct techniques. Yoon and Ghajar [3] pioneered the concept of a stretched sheet by providing a precise solution to the resulting differential system. Recent academics have investigated this topic from a variety of perspectives [412]. Mushtaq et al. [13] investigated the Powell–Eyring fluid flow and heat transport past a stretched sheet exponentially. They discovered that increasing the velocity ratio parameter results in a thinned boundary layer. Malik et al. [12] examined the Powell–Eyring fluid flow and heat transport with varying viscosity over a stretching cylinder by examining the steady condition. They concluded that as Prandtl and Reynolds numbers increase, the boundary layer shrinks. Sher Akbar et al. [14] studied the effect of magnetic factors on Eyring–Powell fluid flow past a stretched surface. They investigated flow resistance as the magnetic and hydrodynamic properties of the fluid under study increased.

Kumar and Srinivas [15] investigated the Powell–Eyring nanofluid passing via an inclined permeable sheet. They demonstrated that temperature increases as thermophoresis parameter values increase. While the contrary is true for nanoparticle concentration due to higher chemical reactions and Brownian parameters, increasing thermophoresis parameter values results in an increase in concentration. Pal and Mondal [16] demonstrated magneto-bioconvection of the Powell–Eyring nanofluid via a vertically stretched sheet that is convectively heated and also contains motile, gyrotactic microorganisms. They discovered that as the Schmidt number and chemical reaction parameters increase, the concentration of nanoparticles drops. Thermal relaxation time is the time required for fluid to return to its original temperature after being heated. It is a frequently used parameter for determining the time required for heat to leave a fluid. Hayat and Nadeem [17] investigated the effects of mass flux models on Eyring–Powell fluid flow in three dimensions. They discovered that temperature and thermal-relaxation time have an inverse relationship. Reddy et al. [18] studied the effect of chemical reactions on the activation energy of the Eyring–Powell nanofluid flow via a stretching cylinder. They concluded that as the relaxation parameter increases, the temperature curves decrease in shape. It takes a long time for an increase in the relaxation parameter assessment to transfer heat to neighboring material particles. Additionally, the Nusselt number improves behavior when nondimensional thermal relaxation calculations are performed.

Mustafa [19] researched the Maxwell fluid with a generalized heat flux model for rotating flow and heat transfer. They also discovered that the thermal relaxation period is inversely proportional to temperature and thermal boundary thickness. Ishaq et al. [20] demonstrated that the entropy production of the Eyring–Powell fluid flow with nanofluid thin film flow can be calculated by considering the heat radiation and MHD impact. They discovered that when the Brinkmann, Hartmann, and Reynolds numbers grow, so does the entropy profile. For increasing values of the Eyring–Powell and radiation parameters, the entropy profile reduces. The Eyring–Powell nanofluid flow with nonlinear mixed convection and entropy generation was explored by Alsaedi et al. [21]. They arrived at the conclusion that entropy generation showed a falling tendency for some fluid parameter values while increasing for others. Through a permeable stretching surface, Bhatti et al. [22] studied the irreversibility of the MHD Eyring–Powell nanofluid. More interesting articles can be seen in [2330] and cross references.

According to the existing literature, no attempt has been made to investigate the electrically conducting Eyring–Powell fluid with radiation, thermal relaxation time, and joule heating effects beyond an exponentially stretched sheet with chemical reaction. This work visually depicts and tabulates the impacts of various flow parameters encountered in the governing equations. The SLM technique is used to solve the issue numerically, which is more computationally efficient. The relevant results are graphed and quantitatively analyzed. This research fills a void in the literature and lays the groundwork for future researchers to contribute their perspectives to the open literature. This is structured as follows: Section 1 contains the literature survey; Section 2 contains the mathematical formulation; Section 3 contains the methodology; Section 4 has the results; and Section 5 contains the conclusion.

2. The Problem’s Formulation

Consider an incompressible Powell–Eyring fluid flowing across an exponentially stretched surface subjected to magnetic, joule heating, thermal radiation, and thermal relaxation periods, as illustrated in Figure 1. The sheet is put on the -and -axes, respectively, and the flow is restricted to Let represent the sheet velocity, represent the external fluid velocity, and represent the surface temperature, with being the ambient temperature.

The governing equations so obtained are given as (see for example, [13], [21], [22], [31]).where are kinematic viscosity, fluid density, velocities, fluid parameters, temperature, thermal conductivity, thermal radiation, specific heat at constant-pressure, strength of the magnetic field, heat flux, the consternation field, diffusion coefficient, and chemical reaction rate, respectively, which satisfy the relation

The appropriate boundary conditions are

Using the similarity transformations as follows:

The continuity equation is satisfied in the same way using (6), and (2)–(5). is transformed into the following form:

Here, and denote velocity ratio and Prandtl number, respectively. Where, Schmidt number, chemical reaction parameter, and are the dimensionless fiuid parameters. Since is a function of , therefore, we use a local similarity solution of (8)–(10) that allows us to analyze parameter behavior. For , we have the case of a Newtonian fluid. The and the local are mathematically described as follows:

Here, and are mathematically described as follows:

The mathematical form of the local Nusselt number and skin friction coefficient are given as follows:where the local Reynolds numbers are .

3. Solution Methodology

Bhatti et al. [22] solved a non-Newtonian model known as the Powell–Eyring fluid model using the collocation approach. Rahimi et al. [23] addressed this model numerically by using a sequential linearization approach and the Chebyshev spectral collocation method. Agrawal and Kaswan [24] solved the Eyring–Powell fluid model using a fourth-order precision methodology (BVP4C) and the homotopy analysis method (H.A.M). Jafarimoghaddam [25] studied the Eyring–Powell model and described fluid flow and heat transfer over a stretching sheet. He then solved the governing PDEs by using homotopy perturbation and homotopy analysis methods to convert them to ODEs. The third-order nonlinear ordinary differential equations (7) and the second-order nonlinear ordinary differential equations (8) are expressed as differential equations and solved using the successive linearization technique (SLM) [26, 31] in this article.

3.1. Procedure of Computational

SLM is used to find the numerical solutions for the nonlinear systems (8)–(10) that conform to the boundary condition (11). We choose the initial guess functions for the SLM solution, i.e., , and are in the form

Here, the two functions and are representative unknown functions. , are successive approximations, which are obtained by recursively solving the linear part of the equation that results from substituting (15) in the governing equations. The mean idea of the SLM is that the assumption of unknown function , , and are very small when becomes larger; therefore, the nonlinear terms in , , and , and their derivatives are considered to be smaller and thus neglected. The intimal guess functions , , and , which are selected to satisfy the boundary conditions

Which are taken to be in the form

3.2. Convergence Analysis

Table 1 illustrates the convergence for the numerical values of the skin friction coefficient, the local Nusselt number, and the local Sherwood number for various values of the parameters involved in using SLM, when .

3.3. Numerical Scheme Testing

Here, we test the validity of our numerical results and contrast them with those of published works as limiting examples. As a result, we compare the current results to those obtained in reference [13], and we discover that they are in reasonable agreement, as shown in Table 2.

4. Result and Discussion

The velocity ratio parameter, the fluid parameter , the magnetic parameter , the nondimensional fluid parameter, and the velocity profile are all monitored for variation. Additionally, this section discusses the influence of the Prandtl number , the velocity ratio parameter, the fluid parameter , the Eckert number , the radiation parameter , the thermal relaxation time , and the magnetic parameter on the dimensionless temperature . Lastly, this section shows the effect of the velocity ratio parameter, the fluid parameter , the magnetic parameter , the Schmidt number , and the chemical reaction parameter on the dimensionless concentration . Two types of boundary layers near the sheet have evolved in a flow with exponentially changing free stream velocity over an exponentially stretched sheet. This means that they are dependent on the velocity ratio parameter , for values of greater than or equal to one. Additionally, it's worth noting that when , no velocity boundary layer arises near the sheet. The velocity profiles for various values are depicted in Figure 2. The influence of the fluid parameter on the velocity is seen in Figure 3. A rise in can be interpreted as either a fall in viscosity or a decline in the Powell–Eyring fluid’s rheological effects. Here, we see that velocity and the thickness of the boundary layer are rising functions of when . This observation leads to the conclusion that the increase in the elastic effects of the Powell–Eyring fluid leads to a thinner momentum boundary layer. However, an opposite trend is noticed when increasing results in a drop in fluid viscosity, which results in an increase in velocity. Additionally, as increases, the viscosity of the fluid becomes lower due to which the increase in the velocity of the fluid accrues. The velocity profile declines as grows but changes toward the border, indicating that the boundary layer's thickness has decreased, which is depicted in Figure 4. As the magnetic field intensity increases, the velocity profile in Figure 5 drops. This is because an increase in the Lorentz force creates resistance to fluid flow, resulting in a drop in the velocity profile.

The fluctuation of the velocity ratio parameter on the temperature profile is depicted in Figure 6. The temperature is discovered to be a decreasing function of . This data may imply that a greater sheet velocity results in a thicker thermal boundary layer. As increases, there is a slight reduction in temperature, as illustrated in Figure 7. Due to the lack of viscous dissipation effects, the fluid parameter is not explicitly included in the energy calculation, and hence has a reduced effect on the thermal boundary layer. Figure 8 illustrates the effect of on temperature . The temperature profile falls as increases. Additionally, rising values of decreases the thickness of the thermal boundary layer. As a result, heat travels rapidly, leading to a decrease in fluid temperature.

The influence of radiation on temperature distributions can be seen in Figure 9. Increases in Rd result in an increase in heat fluxes from the sheet, which results in a rise in temperature. effect on the temperature profile is depicted in Figure 10. As the value grows, the sheet’s wall temperature increases. Due to the fact that when is high, the rate of heat transfer at the surface is low, and the thickness of the thermal boundary layer increases. Frictional heating happens at the surface, raising the fluid’s temperature. The effect of thermal relaxation time on the temperature profile is illustrated in Figure 11. Temperature and thermal relaxation time have been found to have an inverse connection. Physically, when we increase the pressure, the fluid elements have to work harder to transfer heat to their neighboring components, resulting in a temperature drop. When , heat rapidly spreads throughout the fluid. Figure 12 illustrates the effects of the magnetic parameter on the temperature profile. When increases, there is a slight reduction in concentration, as seen in Figure 13. The fluid parameter is not explicitly included in the energy calculation since there are no effects of viscous dissipation, which reduces its impact on the concentration boundary layer. Figure 14 depicts the effect of the magnetic field on dimensionless concentration. The increase in is thought to raise the concentration profile. Figure 15 shows how the velocity ratio parameter varies in relation to the concentration profile. It is shown that the concentration decreases as it increases. According to these findings, a thicker concentration boundary layer is produced by a higher sheet velocity. The effect of the Schmidt number on dimensionless concentration is shown in Figure 16. It is seen that as the Schmidt number increases, the concentration falls. Figure 17 shows how the chemical reaction affected the concentration profile. The concentration decreases as the of the chemical reaction rises.

The local Nusselt number is listed in Table 2, and was estimated using the SLM. In Table 3, the skin friction coefficient increases as increases. As a result, as increases, the coefficient of friction on the skin lowers. According to Mushtaq et al. [13], on an exponentially stretched surface, the magnitude of the skin friction coefficient decreases significantly as the velocity ratio grows. It has already been noted that when k grows, the thermal boundary layer’s thickness decreases. As a result, the heat transfer rate at the stretching sheet is increased. Additionally, as grows, the size of the local Nusselt population decreases dramatically. Additionally, it increases as the values of and increase.

5. Concluding Remarks

In this article, the numerical solution for thermal transport in the Powell–Eyring model via generalized heat flux over an exponentially stretching sheet with a chemical reaction is obtained. By resolving expressions for velocity, temperature, and concentration distributions, the SLM approach is utilized to numerically solve the flow equations. The impact of the Powell–Eyring fluid parameter , magnetic parameter , Eckert number radiation parameter , thermal relaxation time , and chemical reaction was investigated and presented in tables. The validity of the current results was tested, and they were contrasted with those that had previously been published [13]. Table 2 shows a limited example where there is strong agreement. The study’s most important features are listed as follows(i)The velocity increases as the fluid parameter is increased, while reverse behaviour is noticed for the temperature profile.(ii)For increasing values of the magnetic parameter , the velocity profile falls while the temperature rises. In addition, as the resistance to flow increases, the magnetic field intensity and increase.(iii)The temperature and thickness of the thermal boundary layer are inversely related to the thermal relaxation time , whereas the Eckert number has the opposite trend. With an increase in , the temperature within the channel rises.(iv)Increasing values of the (radiation parameter) increase the heat fluxes from the surface, which will cause an increase in the fluid’s temperature and velocity.(v)Simulations of local Nusselt number are verified with published work.


:Velocity components,
:Ratio of expansion rates
:Stretching velocity,
:Dimensionless Powell Eyring fluid parameters
:Velocity of external flow,
:Prandtl number
:Characteristic length,
:Skin friction coefficient
:Surface temperature,
:Local Nusselt number
:Ambient temperature,
:Wall sheer stress
:Stress tensor,
:Surface heat flux
:Kinematic viscosity,
:Local Reynolds numbers
:Dynamic viscosity,
:Magnitude of magnetic field vector,
:Density of the fluid,
:Radiative heat flux
:Powell–Eyring material parameter,
:Mean absorption coefficient,
:Powell–Eyring material parameter,
:Stefan Boltzmann constant,
:Specific heat,
:Dimensionless thermal relaxation time
:Temperature of the/fluid,
:Thermal relaxation time
:Thermal conductivity of the fluid,
:Eckert number
:Radiation parameter
:Magnetic parameter
:Schmidt number
:Chemical reaction parameter.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The author declares that there are no conflicts of interest.


The author acknowledge King Abdul-Aziz University, College of Science & Arts, Department of Mathematics, Rabigh 21911, Saudi Arabia.