Mixed Convection of Non-Newtonian Erying Powell Fluid with Temperature-Dependent Viscosity over a Vertically Stretched Surface

Ahlam Aljabali,Abdul Rahman Mohd Kasim,＊,Nur Syamilah Ari n and Sharena Mohamad Isa

1Centre for Mathematical Sciences,UMP,Gambang,26300,Malaysia

2Faculty Computer and Mathematical Sciences,UiTM Johor,Pasir Gudang Campus,Masai,81750,Malaysia

3Manufacturing Engineering Technology Section,UniKL,Italy Design Institute,Kuala Lumpur,56100,Malaysia

Abstract:The viscosity of a substance or material is intensely influenced by the temperature,especially in the field of lubricant engineering where the changeable temperature is well executed.In this paper,the problem of temperature-dependent viscosity on mixed convection flow of Eyring Powell fluid was studied together with Newtonian heating thermal boundary condition.The flow was assumed to move over a vertical stretching sheet.The model of the problem,which is in partial differential equations,was first transformed to ordinary differential equations using appropriate transformations.This approach was considered to reduce the complexity of the equations.Then,the transformed equations were solved using the Keller box method under the finite difference scheme approach.The validation process of the results was performed,and it was found to be in an excellent agreement.The results on the present computation are shown in tabular form and also graphical illustration.The major finding was observed where the skin friction and Nusselt number were boosted in the strong viscosity.

Keywords:Temperature-dependent viscosity;Erying Powell fluid;numerical solution;combined convection

1 Introduction

Advancements in the study of fluid mechanics have contributed to the significant innovation in engineering devices where it is widely used in the technological and industrial fields.The process in manufacturing these devices involves the heat transfers procedure.Heat transfers also can be found in the boundary layer flow,heat exchangers,and solar core receivers,which are visible to wind currents,fancooling,the manufacturing of electronic devices,cooling nuclear reactors during an emergency shutdown,low-velocity heat exchangers and electrical equipment[1–4].Currently,this research area is emerging and has begun to gain a high interest among researchers due to the considerable applications in the engineering production comprising polymer processing.The advancement in fluid mechanics as well as heat transfer,for example,has continuously been the critical subjects in stipulating the pieces of machinery of manufactured polymers.Polymer properties,on account of their complex configuration,are different compared to viscous fluids like water or oil due to their excessive viscosity and viscoelasticity and that they are commonly non-Newtonian in nature.Some vigorous features of polymers are its retardation and relaxation criterion.In many circumstances,operations involving polymers by most manufacturing processes are generally in the laminar form as a result of their high viscosity level[5].The action of natural buoyancy forces on the fluid flow is termed as free or natural convection.Predominantly,the natural buoyancy forces increase density gradients in the fluid that is induced by[6]and differences in the temperature.Therefore,previous literature has worked on different areas of this respective flow in varying their characteristics,such as geometries,boundary conditions,and fluid-based forms.Such flow patterns exist concurrently in both the external forcing system and internal volumetric forces that are significant either with the influence of the buoyancy forces in forced convection or with the influence of forced flow in free convection.Forced convection occurs when λ →0 while free convection happens when λ →∞.Meanwhile,if the free convection and forced convection flows act together,mixed convection flow will occur.This mixed convection happens due to the difference between surface temperature and the surrounding fluid that results in the power of buoyancy forces.Several studies have mentioned the mixed convection flow for various physical and geometrical aspects[7–12].For example,Umavathi et al.[13]investigated the flow under free convection by considering the combined effects of variable viscosity and thermal conductivity over a vertical channel.An effect of the temperature-dependent fluid properties on the fluid-structure and heat transfer is essential in cases of moderate temperature difference cases.Studies on viscous dissipation on a vertically stretched surface for unsteady mixed convection and heat transfer have been performed by[14–18].In[19],the unsteady natural convection in the differentially heated square cavity of a fluid of temperature-dependent viscosity has been numerically analysed.It has been noticed that thermal convection is more intensive at near the hot wall for a variable viscosity fluid,and the heat input from the hot vertical wall to the cavity exceeds the heat output to the cold vertical wall during the transient phase.An enhanced viscosity correlation has been given by[20].Mass transfer effects on the flow through an accelerated vertical plate,as well as the combined effect of high viscosity and thermal conductivity have been studied by[21]and[22].In[23],the unsteady natural convection was examined with temperature-dependent viscosity within a square porous cavity.In[24],a mathematical analysis of a continuous,incompressible flow of a temperature-dependent on nanofluid from a vertical stretching surface under the external magnetic field and gravitational body force effects has been carried out using the Reynolds exponential viscosity model.The examination of convective flow and heat transport of fluids with temperature-dependent properties is of a great significance for the functional value of modern fluid systems,as mentioned in[25]and[26].Another analysis of the flow problem associated with this particular effect on temperature-dependent viscosity moving over a vertical plate has drawn some researchers’attention[27–29].

Commonly,in industries,the non-Newtonian fluid is widely used,such as Erying Powell,Power law,Maxwell and micropolar fluid models.The model considered is complex and has a preference for other fluid types.The derivation of the Power-Law model is based on the kinetic theory of liquid rather than the empirical relation,which condenses the high and low shear levels according to properties of Newtonian significance.Previously,Erying Powell fluids’ flow has been numerically analysed by[30],which showed that the effects on Erying Powell to the magnetic field over a stretching area.The velocity profile declined as magnetic flux and Erying Powell fluid parameters rise in pressure.They noticed that for high Prandtl number values,the boundary layer would decrease.They also examined the increasing velocity profile by optimising suction parameter values.It displayed that the temperature of fluid decreases for a larger suction parameter.Apart from that,Ara et al.[31]focused on an investigation under exponential representation of the Erying Powell fluid flow over a shrinking surface.They investigated the increase in mass suction that improved the velocity profile,while the temperature profile displayed contrary behaviour.Besides,the thickness of the boundary layer decreased with increasing values of Prandtl number.Moreover,Malik et al.[32]studied the movement of Erying Powel fluid over an expanding cylinder and finally found two versions namely Reynaldo’s and Vogel’s.They found that the boundary layer decreased for high Prandtl number values.They also examined the increase in velocity profile that depended on the suction parameter.For large values of suction parameter,the temperature profile decreased.The nature of mass and heat transfer on free and forced convection flow in a stretching cone has been observed by[33].They noted that the flow parameters’tangential velocity varies in behaviour.Furthermore,the coefficient of skin friction increased due to the growth in the ratio of buoyancy forces.In[34],the relation of Erying Powell fluid series and numerical flow solutions with Newtonian heating has been reported.Meanwhile,the MHD flow of hyperbolic tangent fluid through a stretching cylinder and Keller box system has been studied by[35].In addition,the Erying Powell fluid’s model,has been examined under various physical conditions in[36–39].

According to most industrial and engineering applications,the flow produced by a stretching surface generates keen interest.The analysis showed that the incompressible boundary layer of mixed convection flow over a stretching surface in two dimensions could be considered as a heating environment with the Newtonian heating boundary condition.The present study explores the solution on mixed convective of Eyring Powell fluid’s flow combined with temperature-dependent viscosity and embedded with thermal Newtonian heating past a vertical stretching surface.The numerical computation was performed by using the Keller box approach,and the solutions are presented through graphs and tables.

2 Problem Formulation

Consider the problem of two-dimensional incompressible Eyring Powell fluid under the mechanism of temperature-dependent viscosity embedded with thermal boundary condition known as Newtonian Heating(NH).The flow occurs over a vertical stretching sheet where it stretches in a direction and moves with uniform velocity.Thex-axis is defined as an upward direction alongside the sheet,whiley-axis assumes its upright direction to the surface.The basic theory on the derivation of the Eyring Powell model is constructed on the rate processes theory,which demonstrates the sheer of non-Newtonian flow.

The constitutive equation of Cauchy stress of Eyring Powell Fluid is defined as

where

Apart from that,the second-order approximation of the sinh-1function is defined as

In the Cartesian coordinate system,the governing equations representing Eyring Powell fluid with temperature-dependent viscosity can be expressed as

Continuity equation

x-momentum equation

y-momentum equation

Eqs.(5)and(6)were once investigated by[40]without considering buoyancy and temperaturedependent viscosity.The complexity of momentum equations was reduced by utilizing the boundary layer theory pioneered by Prandtl(1904).Through this theory,the complicated Navier–Stokes equation can be simplified into a set of equations known as the boundary layer equations[41].According to the principle under boundary layer theory,the velocity and temperature differential are more significant than those in the direction at right angles to the ground.Under the explanation on this simplification,the procedure is made available in which Eqs.(5)and(6)take the form

Eq.(8)indicates that the equation is unconditional toy.If the lateral velocity is zero,then the pressure is constant.The power of volume is the forces acting on fluid flow.It is also known as the long-range arm.In general,the frequency of this force varies very slowly and works equally on all areas of fluid flow.Example of this type of forces is gravitational,given by

wheregis the centre of gravity.In two-dimensional flows,the gravitational forces are defined as,

The pressure,P,in the momentum equation,is the combination of hydrostatic pressure,Phwhile dynamic pressure in the quiescent ambient medium,Pd.Sincegis moving downwards andx-direction is upwards,the pressure gradient as described in[42]can be written as

Substituting Eq.(11)into Eq.(7)produces

When performing Boussinesq approximation,all consequences of variable properties are ignored except for density in the momentum equation.This approximation is used as a pure temperature effect to estimate the difference in density,which causes the flow as a result of an interaction between gravitational body force and hydrostatic pressure gradient.The maximum value is believed to be minimal,and the subsequent Taylor expansion becomes

Manipulating the Eq.(13),yields

Therefore,after substituting Eq.(14)into Eq.(12),the continuity and energy equations are given correspondingly,which takes the form as describe in[43]

The present flow analysis is subjected to the following boundary conditions

whereuw(x)is taken as velocity of the stretched surface withaa as positive constant,(u,v),ν,T,T∞,ρ,cp,μ,hs,P,α,?β andcare the velocity components inx-,y-directions,kinematic viscosity,temperature,ambient temperature,density,specific heat at constant pressure,viscosity coefficient,the heat transfer parameter,a pressure,thermal diffusivity and fluid factors of the Powell’ model,respectively.Meanwhile,?β andchave the dimension of(time-1).

To obtain the set f similarity equation in the form of an ordinary differential equation,the similarity transformation variables as in Eq.(19)are adopted and applied to the Eqs.(15)to(17),which yields

To feign the temperature-dependent viscosity variable,the Reynolds exponential viscosity model is taken into account,which gives a detailed approach to

where fluid viscosity is assumed to vary as a linear inverse temperature function

The resulting equations are obtained as follows

together with the transformed boundary conditions as

3 Computation Procedures

The obtained solvable equations are resolved using Keller box technique in the form of established ordinary differential equations.Eqs.(22)and(23)according to boundary conditions(24)are condensed to a system of the first order.For that matter,the independent variables are demonstrated as

Then the respective equations can be written as

In Eqs.(19),(22)to(24),the notation prime(′)corresponds to the derivative with respect to η.Additionally,where the fluid parameters,MandB,Prandtl number,Pr,conjugate parameter,γ,viscosity parameter,α,fluid kinematic viscosity,νfand mixed convection,λ,are defined as,

The steps of Keller box techniques are as follows

i)The Eqs.(22)and(23)are reduced to a first-order arrangement.

ii)The attained system in(i)is changed to a system of algebraic equations by means of the central difference formula.

iii)The resultant in(ii)is linearized using Newton’s method and prepared in matrix-vector system.

iv)The output(iii)is determined using the block-tridiagonal elimination method.

In this present mathematical model,the magnitudes of the skin friction coefficient and the Nusselt number are as follows

4 Validation Procedure

The present model is computed using the Keller box approach.The upshot of fluid parameters(MandB)representing Eyring model,Prandtl number(Pr),mixed convection parameter(λ),conjugate parameter(γ)and viscosity parameter(α)on velocity and temperature profile against η are figured via MATLAB software.The proposed method is undeniably suitable since it is proven to be unconditionally stable,even for a higher-order model.It is worth mentioning here that,the comparison between the present results with the previous existing solution for limiting case is essential to claim that the current model and its output are acceptable.

The direct comparative studies with the existing literature have been conducted to corroborate the numerical solutions acquired in this study.The Tab.1 shows the summary of the work done by[24].The momentum equation on particular works can be reduced to its limiting case by fixing theBr=0 andM=0.

Table 1:Model by[24]

The output from the comparative study on momentum equation shows a firm agreement,indicating that the current model and its findings are acceptable.A clear picture of this working is presented in Tab.2.Besides,the present works also presents a comparative study on the value of θ(0)and -θ′(0)and the results were undeniably in agreement with the reports by[44]and[45].The results on this particular matter are presented in Tab.3.

Table 2:Comparison of value(1-αθ(0 ))f′′(0)at λ=1,B=M =0

5 Results and Discussion

The numerical solutions obtained were analysed after undergoing proper validation process.The results were computed under the NH boundary condition for various parameters Pr,α,M,B,γ and λ.The computation was done by assigning a set of fixed values of the parameter followed by the calculation on the value ofand.

Table 3:Comparison on value θ(0)and -θ′(0)at B=λ=0,M =γ=1

Tab.4 depicts the variations ofandwith variation of value of pertinent parameters involved.For the entire computation,the value of γ was set as 1.It was noticed that a more substantial value of Pr,λ,γ,α,MandBboosted the.This was due to the presence on those parameters that made the convective process in the fluid flow become more active.Thewas increased in the higher values of γ and α but showed a reduction trend in the larger values of Pr,λ,MandB.The increasing trend happened due to the strengthen on drag forces,which increased the viscosity of the fluid.Meanwhile,the lessening was due to the weakness in the resistance of fluid with the surface.

Table 4:Numerical results of and various values of Pr,α,M,B,λ and γ

Table 4:Numerical results of and various values of Pr,α,M,B,λ and γ

PrαMBγλθ(0)Cf Rex1/2NuxRex-1/2 70.110.110.10.381319-1.3878820.618680 9 0.331691-1.3913230.668308 120.252203-1.3949750.747796 100.150.317005-1.3897970.682994 0.30.303985-1.3828030.696014 0.40.292714-1.3785980.707285 100.10.50.322521-1.0732250.677478 0.90.315830-1.2269560.684169 1.50.261353-1.4654830.738646 100.110.30.314332-1.3536980.685667 0.60.289651-1.3005270.710348 0.80.289522-1.2672000.710477 100.10.110.70.238458-0.9916740.533079 1.20.348424-0.9862010.781890 1.70.431920-0.9818950.965734 100.110.910.30.255546-1.2482680.744453 0.50.254600-1.2413040.745398 0.70.247357-1.2404690.752642

Furthermore,the behaviours on various values of parameters Pr,α,M,B,γ and λ on the velocity and temperature of fluid are depicted from Figs.1 to 12.An increasing trend on velocity can be seen in the larger value ofMand λ as captured in Figs.5 and 11.This behaviour shows that the properties of those parameters characterises the ability to exchange thermal energy with its surroundings.Meanwhile,a contradicting behaviour can be seen for larger values of Pr,B,α and γ as illustrated in Figs.1,3,7,and 9.The temperature distribution showed increasing trend for the parameters α,Pr and λ while the opposite trend was observed for γ,MandB.These trends happened due to the respective parameters with properties to absorb and release the energy to the fluid.

All the figures indicate,at far from the plate,that the velocity and temperature profiles of fluid asymptotically approached zero,which has fittingly fulfilled the boundary condition.

Figure 1:Variation of f ′(η)at M =γ=1,B=0.1,λ=3 and α=0.3 for various values of Pr

Figure 2:Variation of θ(η)at M =γ=1,B=0.1,λ=3 and α=0.3 for various values of Pr

Figure 3:Variation of f′(η)at B=λ=0.1 and Pr=10 for various values of α

Figure 4:Variation of θ(η)at M =γ=1,B=λ=0.1 and Pr=10 for various values of α

Figure 5:Variation of f′(η)at B=0.5,α=0.1,γ=1,λ=0.01 and Pr=10 for various values of M

Figure 6:Variation of θ(η)at B=0.5,α=0.1,γ=1,λ=0.01 and Pr=10 for various values of M

Figure 7:Variation of f′(η)at M =γ=1,α=λ=0.1 and Pr=10 for various values of B

Figure 8:Variation of θ(η)at M =γ=1,α=λ=0.1 and Pr=10 for various values of B

Figure 9:Variation of f′(η)at λ=3,Pr=10 and M =B=α=0.1 for various values of γ

Figure 10:Variation of θ(η)at λ=3,Pr=10 and M =B=α=0.1 for various values of γ

Figure 11:Variation of f′(η)at M =γ=1,B=α=0.1 and Pr=10 for various values of λ

Figure 12:Variation of θ(η)at M =γ=1,B=α=0.1 and Pr=10 for various values of λ

6 Conclusion

This present research has integrated the mixed convective flow of an Eyring Powell fluid over a vertically stretched surface that was reviewed by highlighting the effects of temperature-dependent viscosity involving the parameters of Pr,α,M,B,γ and λ.From the mathematical analysis,a similar trend can be noticed in the motion and temperature distributions of fluid,respectively,when parameters were increased.The results of the parameters are illustrated in the specified graphs and tables.Nevertheless,the variations on the velocity distribution as well as fluid temperature portrayed the importance of the parameters investigated.The findings in this study contribute to a better understanding of the characteristics of fluid flow and its advancement.

Funding Statement:This project has been supported by Ministry of Higher Education and Universiti Malaysia Pahang through RDU182307.

Conflicts of Interest:The authors declare no conflicts of interest to report regarding the present study.