# Optimal Homotopy Asymptotic Analysis of the Dynamics of Eyring-Powell Fluid due to Convection Subject to Thermal Stratification and Heat Generation Effect

Authors: Solomon Bati Kejela; Emiru Rufo Adula
Journal: Abstract and Applied Analysis
 (2022)

---

## Abstract

In the present study, the effect of thermal stratification and heat generation in the boundary layer flow of the Eyring-Powell fluid over the stratified extending surface due to convection has been investigated. The governing equations of the flow are transformed from partial differential equations into a couple of nonlinear ordinary differential equations via similarity variables. The optimal homotopy asymptotic method (OHAM) is used to acquire the approximate analytical solution to the problems. Impacts of flow regulatory parameters on temperature, velocity, skin friction coefficient, and Nusselt number are examined. It is discovered that the fluid velocity augments with a greater value of material parameterE, mixed convection parameter λ, and material fluid parameter σ. The result also revealed that with a higher value of the Prandtl number Pr and the stratified parameter ε, the temperature and the velocity profile decreases, but the opposite behavior is observed when the heat generation/absorption parameter γ increases. The results are compared with available literature and are in good harmony. The present study has substantial ramifications in industrial, engineering, and technological applications, for instance, in designing various chemical processing equipment, distribution of temperature and moisture over agricultural fields, groves of fruit trees, environmental pollution, geothermal reservoirs, thermal insulation, enhanced oil recovery, and underground energy transport.


---

## Body

## 1. Introduction

The study of the dynamics of non-Newtonian fluids has inconceivable importance on account of its numerous biological, industrial, and engineering applications. They have a vital role in material production, mixture, paper coating, bioengineering, oil repository designing, petroleum production, food engineering, chemical engineering, and the study of electroosmatic flow, blood flow, etc. [1–8]. Non-Newtonian fluid can be found in numerous materials for example in shampoos, paints, cosmetic productions, body fluids, grease, and animal blood to name a few. Several models including the Eyring-Powell model (one among the non-Newtonian fluid models) have been proposed by different scientists for the analysis of dynamics of non-Newtonian fluids. The Eyring-Powell model is derived from the kinetic molecular theory of liquids rather than empirical relation and properly shows the Newtonian behavior at low and high shear rates; owing to this, it gained the attention of many scholars. This class of fluid does not have constant viscosity with temperature alterations [9–13].The Eyring-Powell fluid model has essential applications in designing various chemical processing equipment, distribution of temperature and moisture over agricultural fields, groves of fruit trees, environmental pollution, geothermal reservoirs, thermal insulation, enhanced oil recovery, and underground energy transport [11]. The Eyring-Powell model with heat transfer plays a crucial role in diverse industrial, geophysical, and natural processes. The understanding of the heat transfer effect in a variety of materials and boundary layer flow overstretching sheets is also imperative as heat has substantial applications in diverse power station engineering, biological phenomena, and industrial branches like paper manufacture, metal extrusion, and bubble absorption [14–18].Heat and mass transfer in the Eyring-Powell liquid model undertakes an essential role in the procedures which comprise the creation and spread of haze, plotting of mixture handling instrumentation, environmental pollution, drying of porous slides, raised oil recuperation, warm protection, and underground energy transport [19–21]. Ishaq et al. [22] explored the 2D nanofluid film stream of the Eyring-Powell liquid with changeable heat conduction in the existence of MHD over the unsteady permeable extending sheet and asserted that the porosity parameter lessens the movement of the fluid.Hayat et al. [23] conducted a study on MHD nonlinear extending flow of the Eyring-Powell nanofluid, where a zero nanoparticle mass flux stipulation was applied at the sheet. Their finding illustrates that with magnetic parameters, skin friction coefficient raises though the heat transmits rate at the surface declines. The heat and mass transfer analysis of non-Newtonian fluid overstretching sheets has been done by different scholars; Megahed et al. [24] conducted a study on heat transfer for MHD fluid as a result of the unsteady stretching plate with extended heat flux by varying viscosity and thermal conductivity with temperature; Reddy et al. [25] and Megahed and Gnaneswara [26] have also investigated the heat and mass transfer of different non-Newtonian flow across a stretched sheet.Mixed convection flow plays a significant role, for example, in air limit layer flows, atomic reactors, solar gatherers, and in electronic hardware, just to name a few. Wubshet and Temesgen [19] used Cattaneo–Christov heat and mass flux model to inspect the mixed convection flow of the Eyring-Powell nanofluid through a porous medium on a linearly extending sheet in the company of permeability and Joule heating effects. They noticed that with the advanced value of the Eyring-Powell fluid parameter, the fluid velocities, local Nusselt number, and skin friction coefficient promotes. Malik et al. [27] have examined a mixed convection flow for the MHD Eyring-Powell fluid over the extending plate solving the problem numerically by the shooting method. Their outcome indicates swelling of the Eyring-Powell material parameter results in a meaningful decrease in both heat and mass transfer. Moreover, Rahimi et al. [28] presented a collocation method solution for an Eyring-Powell fluid flow over a prolonging plate. They found that the velocity enhances with mounting the Eyring-Powell fluid material parameter, but the opposite behavior is observed when the fluid material parameter is augmented.Stratification plays a vigorous role in several industrial and natural phenomena. Owing to its wide uses in fluid mechanics, several investigators studied the stratification phenomenon. Zubair et al. [29] investigated the attributes of a chemical reaction and twofold stratification in the Eyring-Powell fluid flow persuaded by a nonlinear prolonging plate with changeable thickness. They confirmed the temperature decline with thermal stratified parameters; Salahuddin et al. [30] inspected the vital aspects of the stratification procedure with linear convection. They discovered that with an enhanced value of the thermal stratification parameter the velocity, the temperature profile reduces. Hayat et al. [31] examined the MHD effect for heat and mass transfer through the double stratified sheet. They reported that thermal and solutal stratification parameters have a conflicting effect on temperature distributions. Furthermore, Ogunseye et al. [32] examined the effect of thermal radiation on the flow and heat transfer of the Eyring-Powell nanofluid. The understanding of the heat generation/absorption effect of boundary layer flow is also crucial in different fields of study and applications. Owing to this, a lot of scholars [33–35] conducted a study on the heat generation effect on boundary layer flow and heat transfer due to different non-Newtonian fluids. The study of the boundary layer flow of the Eyring-Powell fluid over a stretching surface plays a vital role in industrial applications. A stretching surface has abundant significance in numerous industrial applications including electronic chips, fiber yarn, polymer industries, and glass blowing.Motivated by the above significant studies, we consider the problem of the Eyring-Powell fluid flow over a stratified stretching surface subject to mixed convection and heat generation/absorption effect. The problem has plenty of applications in various fields of science and engineering but it has not been analytically investigated so far. Moreover, the effect of heat generation and thermal stratification on the flow of the fluid under consideration is not investigated. Hence, the novelties of the current work are the investigation of the heat generation and thermal stratification effects in the boundary layer flow of the Eyring-Powell fluid over the stratified extending surface due to convection, the use of OHAM for solving the problem, and the analysis of effects of different parameters on velocity and temperature of the considered fluid. The effectiveness and reliability of the proposed method have been proved in many research articles, and it is used for solving both nonlinear and linear problems by several authors [34, 36–38].


## 2. Mathematical Formulation

In the model, a two-dimensional incompressible steady boundary layer mixed convection flow of the Eyring-Powell fluid over the stratified extending surface has been pondered. The sheet is elongated at invariable speeduw=ax along the x-axis in a vertical direction and the fluid temperature at and far away from the plate surface is given, respectively, by TwandT∞. The influence of heat generation/absorption is considered in the energy equation. The geometric representation of the problem under consideration and the coordinate system are displayed in Figure 1.Figure 1
Physical model of the problem.Applying the boundary layer estimation the governing equation of the problem is condensed into the following form (see [21, 39, 40]):Continuity equation:(1)∂u∂x+∂v∂y=0.Momentum equation:(2)u∂u∂x+v∂u∂y=υ+1ρβc−12ρβc3∂u∂y2∂2u∂y2+gβT−T∞,where u and v, x and y, respectively, represent the horizontal and vertical components of fluid velocity and the coordinate system, β and c are the dimensional Eyring-Powell fluid material parameters, and υ,g represent kinematic viscosity and gravitational acceleration, respectively.Energy equation:(3)u∂T∂x+v∂T∂y=kρcp∂2T∂y2+QρcpT−T∞.The corresponding boundary conditions of the governing flow problem are (see [16, 31]) as follows:
(4)u=uw=ax,v=0,T=Tw=T0+b1xaty=0u⟶0,T⟶T∞=T0+b2xasy⟶∞,where a is a linear stretching rate, b1 and b2 are positive dimensional constants, and T0 is the stretching sheet reference temperature.To make the arithmetical analysis simple, we convert Equations (1)–(4) into ODEs, by introducing the stream function ψx,y, such that
(5)u=∂ψ∂y,(6)v=−∂ψ∂x,where
(7)ψ=aυxfη.The dimensionless similarity variableη, stream function fη, and temperature θη are defined as follows (see [21, 39, 40]):
(8)η=yaυ,(9)θη=T−T∞Tw−T0,where T=b1θxandT∞=T0+b2x.Now, Equation (1) is identically pleased with u and v of Equation (5). Furthermore, using Equation (5) and (8), Equations (2) and (3) become converted into the following ODEs:
(10)1+∈f′″+ff″−σ∈f″2f′″−f′2+λθ=0,(11)θ″+Prfθ′+Prγθ−Prf′θ−Pre1f′=0,where ∈=1/ρβcυ,σ=uw2a/2υc2,λ=gβb1/a2,Pr=μcp/k,γ=Q/acpρ,ande1=b2/b1 represent a material parameter, fluid parameter, mixed convection parameter, Prandtl number, heat generation/absorption parameter, and thermal stratification parameter, respectively. Moreover, the prescribed side conditions Equation (4) become reduced to
(12)f0=0,f′0=1,θ0=1−e1aty=0f′η=0,θη=0asη⟶∞.The physical quantities of interest: local skin friction coefficientCf and the Nusselt number Nu, are described as follows (see [30, 40]):
(13)Cf=τwρuw2,(14)Nu=xqwkTw−T∞,where uw represents free stream speed, τw and qw given by (see [21])
(15)τw=μ+1βc∂u∂y−16βc3∂u∂y3y=0,qw=−k∂T∂yy=0,denote skin shear stress on a surface and heat flux at the surface.From Equations (5), (8), (13), and (15), we obtain
(16)RexCf=1+∈f″0−13∈σf″03,Rex−1Nu=−θ′0,where Rex=xuw/υ=ax2 denotes the Reynolds number and x is the distance from the reference point at which a boundary layer begins to form.


## 3. Solution of the Problem via the OHAM

Employing the basic principle of OHAM on Equations (10) and (11), we construct the homotopy equation:
(17)1−p1+∈f′″=H1p1+∈f′″+ff″−σ∈f′′2f′″−f′2+λθ1−pθ″=H2pθ″+Prfθ′−Prf′θ+Prγθ−Prεf′,where p∈0,1 and f,θ,H1p,H2p are expanded as follows:
(18)f=f0+pf1+p2f2+p3f3,θ=θ0+pθ1+p2θ2+p3θ3,H1p=pc1+p2c2+p3c3,H2p=pc4+p2c5+p3c6,where H1pandH2p are nonzero ancillary functions for p≠0 and zero for p=0. Using Equation (18) in Equation (17) and collecting like terms based on powers of p, we get the next zero-, first-, second-, and third-order problems with their boundary conditions.The zeroth-order problem:(19)f0′″η=0,θ0″η=0.Boundary conditions:(20)f00=0,f0′0=1,θ00=1−e1,f0′η=0,θ0η=0asη⟶∞.The solution for Equation (19) with BCs (20) is
(21)f0η=η1−η20,θ0η=1−ε1−η10.The 1st-order problem:
(22)f1′″η,a1=a11+∈f0f0″−f0′2+λθ0,θ1″η,b1=Prb1f0θ0′−f0′θ0+γθ0−εf0′.Boundary conditions:(23)f10=0,f1′0=0,θ10=0,f1′η=0,θ1η=0asη⟶∞.Solution for Equation (22) with BC Equation (23) is
(24)f1η,a1=a11+∈−η512000+η4240−η36+158η2+λ3e1−1η480−η32+5η2,θ1η,b1=Prb12ε−1η41200−η330+η2−15η2+η2−η330−20η3γ−γε−ε.The 2nd-order problem:
(25)f2′″η,a1,a2=1+a1f1′″+11+∈a1f0f1″+f0″a2f0+a1f1−σ∈a1f0′′2f1′″−2a1f0′f1′+a2f0′2+λa1θ1+a2θ0,θ2″η,c4,c5=1+b2θ1″+Prb1f0θ1′+f1θ0′+b2f0θ0′−b1f1′θ0+f0′θ1−b2f0′θ0+γb1θ1+b2θ0−εb1f1′+b2f0′.Boundary conditions:(26)f20=0,f2′0=0,θ20=0,f2′η=0,θ2η=0asη⟶∞.The solution for the problem is(27)f2η,a1,a2=1+a1a11+∈−η512000+η4240−η36+15η28+λ1−εη32−η4240−5η23+5a121+∈2η84032000+η61200−η750400−19η5960+5η432−275η2336+λ1−ε7η5360−η61200−5η436+η784000+11+∈a210η51200−η424+25η24−a12101+∈−η84032000+η750400−η6720+η532−19η256+λ1−εη6720−η750400−η512+3125η218+−σ∈a121001+∈2−η512000+η4240−η36+15η28+λ1−εη36−η4240−5η23+λa21+∈1−εη36−η4240−5η23−2a121+∈2−η7504000+η67200−η5120+5η432−125η29+λ1−εη5120−η67200−5η436+100η29+2a12101+∈2−η8806400+η712600−η6240+η516−13625η2252+λ1−εη6240−η712600−η518+125η23−a21+∈η36−η4120+η56000−5η24+λa1Prb11+∈ε−1η7504000−η67200+η5120−5η432+125η29+γ−γε−εη5120−η6720−5η436+100η29.(28)θ2η,b1,b2=Pr1+b1b1ε−1η42400−η360+η22+25η9+γ−γε−εη22−η360−10η3+Prb12ε−1η614400−η71008000−η5600+η464−275η56+γ−γε−εη436−η5600−10η39+η636000−225η2+Pra1b1ε−11+∈η71260000−η618000+η5400−11η4192+5η38−7125η280+λ1−εη624000−η5400+η418−5η39+625η30+Prb2ε−110η36−η4240−25η2+Prb2ε−1η22−η330+η41200−5η2+γ−γε−εPrb2η22−η360−10η3+Prb12γε−1η672000−η51200+η424−5η38+250η9+γ−γε−εη424−η51200−5η39+200η9−Pra1b1ε1+∈−η672000+η51200−η424+5η38−250η9+λ1−εη424−η51200−5η39+200η9The 3rd-order problem:
(29)f3′″η,a1,a2,a3=1+a1f2′″+a2f1′″+11+∈a1f0f2″+f1″a2f0+a1f1+f0″a2f1+a3f0−σ∈f1′″a2f0′′2+2a1f0″f1″−σ∈a1f0′′2f2′″−a12f0′f2′+f1′2−2a2f0′f1′−a3f0′2+λa1θ2+a2θ1+a3θ0θ3″η,b1,b2,b3=1+b4θ2″+b2θ1″+Prb3f0θ0′+b2f0θ1′+f1θ0′+b1f0θ2′+f1θ1′+f2θ0′−b3f0′θ0+b2f0′θ1+f1′θ0+b1f2′θ0+f1′θ1+f0′θ2+γb1θ2+b2θ1+b3θ0−εb1f2′+b2f′1+b3f0′.Boundary conditions:(30)f30=0,f3′0=0,θ30=0,f3′η=0,θ3η=0asη⟶∞.Solution for Equation (29) subjected to BCs Equation (30) is achieved in the same way as the zeroth-, first-, and second-order problems. However, there are very bulky expressions to be included herein this paper. Collecting the solution of Equations (21), (24), (28), and (29), we get the four terms solution given as
(31)f∼η,a1,a2,a3=f0η+f1η,a1+f2η,a1,a2+f3η,a1,a2,a3,θ∼η,b1,b2,b3=θ0η+θ1η,b1+θ2η,b1,b2+θ3η,b1,b2,b3.


## 4. Result and Discussion

The OHAM is used to solve the transformed nonlinear ODEs describing the problem. The effect of a variety of parameters on skin friction coefficientCf, Nusselt number Nu, temperature, and velocity has been investigated. The obtained results are displayed with the help of graphical illustrations and tables. Table 1 displays the obtained value of Nu and Cf for particular values of certain parameters and a comparison has been made with published results of others’ work. From the table, it is detected that Cf increases with swelling values of the Eyring-Powell fluid material parameter E and mixed convection parameter λ, but it falls off with an augmented value of thermal stratification parameter ε and material fluid parameter σ. It is also found that Nu increases with enhancing mixed convection λ and Eyring-Powell fluid parameter σ.Table 1
Comparison of results for the skin friction coefficient and Nusselt number for various value ofε,λ,σ,andE when Pr=0.7 and γ=0.1 are fixed.
EσλεPresent studySalahuddin et al. [30]Skin friction−θ′0Skin friction−θ′00.10.10.10.1-0.99480.9224-0.99530.92010.30.10.30.3-0.92350.9448-0.92210.94500.50.10.50.5-0.86320.9661-0.86290.96570.10.30.10.1-0.99870.9194-0.99880.91960.10.50.10.1-0.99310.9187-0.99220.91920.10.10.30.1-0.89730.9482-0.89690.94780.10.10.50.1-0.80460.9689-0.80450.96950.10.10.10.3-0.99961.0221-0.99981.02160.10.10.10.5-0.87981.1223-0.00431.1216The impact of the thermal stratification parameterε and the Prandtl number Pr on Nusselt number Nu and the skin friction coefficient Cf for the fixed value of the other parameters is available in Table 2. The table illustrates Cf is a decreasing function of Prandtl number Pr and the thermal stratification parameter ε. Table 2 also indicates that the local Nusselt number Nu enhances with the increasing value of Pr and ε.Table 2
Comparison of the present result for local skin friction coefficient and Nusselt number for different values of Pr andε when γ=σ=λ=E=0.1 are fixed.
Present resultBilal and Ashbar [39]Ogunseye et al. [32]PrεSkin friction−θ′0Skin friction−θ′0Skin friction−θ′00.10.1-0.9545380.281633——-0.9545290.2816310.30.3-0.9729170.524666-0.9729210.524669-0.9729080.5246680.50.3-0.9859530.734565-0.9859510.734568-0.9859520.7345670.70.5-0.9892860.804125-0.9998811.021608-0.9892850.8041250.70.7-1.0043141.121622−1.0043131.121621-1.0043091.1216220.70.9-1.0086141.220177-1.0086131.220178——0.91.0-1.0213721.237912——-1.0213681.237908The effect of a variety of parameters on the fluid velocity profile is displayed in Figures2–7. In Figure 2, the impact of mounting the material parameter E on dimensionless fluid velocity is displayed. The figure points out that the velocity increase with augmenting value of E. Physically, E has an inverse relation with the non-Newtonian fluid’s dynamic viscosity. Hence, with increased value of E, the flow resistance reduces, and consequently, the fluid velocity intensifies. Figure 3 is sketched to illustrate the influence of mixed convection parameter λ on velocity profile. It depicts that an intensification in λ results in the rise in the velocity profile. In reality, while fluids are heated their density reduces and buoyancy increases, as a consequence their velocity boosts. Thus, by advancing value λ, one can enhance the velocity within a boundary layer when required.Figure 2
Impact ofE on velocity profile when λ=0.4,Pr=0.72,ε=σ=γ=0.1.Figure 3
Impact ofλ on velocity profile when E=0.3,Pr=0.72,ε=σ=γ=0.1.Figure 4
Impact ofε on velocity profile when λ=0.4,E=0.3,Pr=0.72,σ=γ=0.1.Figure 5
Impact ofγ on velocity profile when λ=0.4,E=0.3,Pr=0.72,ε=σ=0.1.Figure 6
Impact ofσ on velocity profile when λ=0.4,E=0.3,Pr=0.72,ε=γ=0.1.Figure 7
Impact of Pr on velocity profile whenλ=0.4,E=0.3,ε=σ=γ=0.1.The influence of thermal stratification parameterε on the fluid velocity is illustrated in Figure 4. It is witnessed that with higher value of ε, the velocity profile declines. Actually, when ε increases, the convective potential between the heated surface and ambient temperature reduces. The effect of heat generation parameter γ on velocity is exhibited through Figure 5. As can be seen from the figure, the velocity profile enlarges with the value of γ. In reality, while the heat generation parameter rises, additional heat is produced in the fluid; consequently, the fluid turn out to be hotter and less dense, and as a result, its viscosity decreases. Figure 6 displays the effect of material fluid parameter σ on the velocity profile. It is confirmed that the velocity graph increases with the mounting value of σ. In actual fact, σ is inversely related to its viscosity. So, the higher the fluid parameter σ leads to the deterioration of the viscous forces; for this reason, the fluid flows rapidly. Figure 7 describes the impact of the Pr on the velocity profile. It is detected that an increment in Pr corresponds to a decrease in the velocity profile of the fluid. Physically, Pr is directly related to the viscosity of the fluid. Because of this fact, as the Pr increases, the fluid turns out to be denser, which causes a decrease in the velocity profile. This may be accredited to the fact that convection currents became diminutive at higher Pr and leads to a decline in velocity.The effect of a variety of parametric quantities on the nondimensional temperature distribution of the fluid under consideration is displayed in Figures8–13. The characteristic of stratification parameter ε on the temperature is depicted in Figure 8. Physically, the rise in ε declines the temperature distinction between the fluid on the surface and the surrounding fluid, this in turn trim down the temperature profile. Thus, the temperature distribution in the fluid drops with an augmentation in the value of ε. Figure 9 reports the effect of the Prandtl number Pr on the temperature profile. It is revealed that when the value of Pr is augmented, the temperature graph reduces. Mathematically, Pr is inversely related to thermal conductivity; hence, the fluid with greater Pr has a low energy diffusion rate. So, the raise in the Pr produces a sturdy reduction in the temperature, which upshots the thermal boundary layer thickness. Therefore, in heat transfer, Pr is used to rheostat the thicknesses of thermal and momentum boundary layers.Figure 8
Impact ofε on temperature profile when Pr=0.72,λ=0.4,E=0.3,σ=γ=0.1.Figure 9
Impact of Pr on temperature profile whenλ=0.4,E=0.3,ε=σ=γ=0.1.Figure 10
Impact ofγ on temperature profile when Pr=0.72,λ=0.4,E=0.3,ε=σ=0.1.Figure 11
Impact ofσ on temperature profile when Pr=0.72,λ=0.4,E=0.3,ε=γ=0.1.Figure 12
Impact ofE on temperature profile when Pr=0.72,λ=0.4,ε=σ=γ=0.1.Figure 13
Impact ofλ on temperature profile when Pr=0.72,E=0.3,ε=σ=γ=0.1.Figure10 is plotted to scrutinize the impact of heat generation/absorption parameter γ on the fluid temperature. The figure confirms that when the value γ intensifies, the temperature profile also increases. Physically, higher γ generates extra heat in the fluid system, due to this the temperature profile enhanced. Thus, when the heat generation parameter upshots, considerable heat which can increase the temperature profile of the fluid will be produced. Therefore, increasing the heat generation/absorption parameter is used to enhance the temperature of the fluid, but sometimes, the augmented values of heat generation parameter results in overmounting in the temperature boundary layer and it should be controlled by using suitable mechanisms.The characteristic of material fluid parameterσ on the fluid temperature is displayed in Figure 11. It is witnessed that an increment in the value of σ corresponds to a turn down in the temperature profile. This behavior is revealed as σ is inversely associated to the viscosity of the fluid. So, when σ increases, viscosity reduces which in turn increases the heat transfer rate. Consequently, the temperature of the fluid reduces. Figure 12 portrays the behavior of material parameter E on temperature profile. As we can observe from the figure, an improvement of value of E results in a decrease in the fluid temperature graph. Physically, when E is increased, fluid resistance reduces and contributes to the reduction of temperature distribution. Figure 13 depicts the impact of varied convection parameter λ on temperature. It has been found that when the value of λ increases, the temperature profile declines. This happens attributable to the reduction of the fluid density, as hotter fluids are less dense. Thus, with enhanced value of λ heat transfer rate increases, thereby temperature profile decreases.


## 5. Conclusion

In this work, we analyze heat transfer and flow of the incompressible, steady Eyring-Powell fluid flow over a two-dimensional stratified prolonging surface with mixed convection. The nonlinear PDEs governing the flow problem and their respective prescribed side conditions have been converted into a couple of nonlinear ODEs via variable similarity transformation and then solved analytically by the OHAM. The quantities of attention in this work are velocityf′, temperature θ, skin friction coefficient Cf, and Nusselt number Nu. Thus, the effect of varied parameters embedded in the model on f′, θ, Cf, and Nu has been examined and described with the aid of tables and graphs. The main observations drawn from this finding are the following:
(i)
With an intensification in value of the thermal stratification parameterε, both the fluid velocity and the temperature profile decreases. Hence, in various manufacturing processes ε can help to control heat and velocity of fluids(ii)
As the value of the Eyring-Powell fluid material parameterE upsurges, the velocity profile boosts, while the temperature profile declines(iii)
With snowballing values of the mixed convection parameterλ and the material fluid parameter σ, the velocity profile enhances(iv)
As the value of Prandtl number Pr enlarges, the velocity profile dwindles(v)
The temperature distribution turn down with enhancing values of the Prandtl number Pr, the material fluid parameterσ, and the mixed convection parameter λ(vi)
As the value of the heat generation/absorption parameterγ amplifies, both the temperature and the velocity profile increase(vii)
As the values of the mixed convection and the Eyring-Powell fluid material parameters heighten, both the Nusselt number Nu and the skin friction coefficient augment(viii)
With mounting values of the thermal stratification and the material fluid parameters, the skin friction coefficient falloffs(ix)
The local Nusselt number Nu is an increasing function of Prandtl number Pr and the thermal stratification parameterε, the mixed convection parameter λ, and the Eyring-Powell fluid material parameter E



---

*Source: 6777513-2022-09-23.xml*