Analysis of Convective Transport of Temperature-Dependent Viscosity for Non-Newtonian Erying Powell Fluid:A Numerical Approach

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

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

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

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

Abstract:Non-Newtonian is a type of fluid that does not comply with the viscosity under the Law of Newton and is being widely used in industrial applications.These include those related to chemical industries,cosmetics manufacturing,pharmaceutical field,food processing,as well as oil and gas activities.The inability of the conventional equations of Navier–Stokes to accurately depict rheological behavior for certain fluids led to an emergence study for non-Newtonian fluids’models.In line with this,a mathematical model of forced convective flow on non-Newtonian Eyring Powell fluid under temperature-dependent viscosity(TDV)circumstance is formulated.The fluid model is embedded with the Newtonian heating(NH)boundary condition as a heating circumstance and is assumed to move over a stretching sheet acting vertically.Using appropriate similarity variables,the respective model was converted into ordinary differential equations(ODE),which was later solved utilizing the Keller box approach.The present model is validated by comparing the existing output in literature at certain special limiting cases,where the validation results display a firm agreement.The current outputs for the proposed model are shown in tabular and graphical form for variation of skin friction plus Nusselt number,velocity and temperature distribution,respectively.

Keywords:Temperature-dependent viscosity;erying powell fluid;vertical stretching sheet;numerical solution;boundary layer

1 Introduction

Over recent decades,several wide-ranging investigations have been conducted on heat transference in a non-Newtonian fluid.Fluids with a non-Newtonian nature are partly or wholly made up of macromolecules or dual-phase materials.These fluids are typically described as those equipped with the capacity chronology of flow.As a rule,the categorization of the numerous forms of non-Newtonian type of fluid is realized on the shear stress equation.Unlike Newtonian fluid behaviors,industrial materials,such as polymeric liquids,paste,melts,emulsions,muds,foodstuffs,condensed milk,soaps,shampoos,molten plastics,etc.,are able to follow the characteristics of non-Newtonian fluid.Therefore,the fluid flow modeling and estimation have fundamental and practical importance in the applications of engineering and its businesses.The behaviors rheological of such fluid flow are substantial inside the flows of nuclear gasoline slurries,paper coating,greases and lubrication with heavy oils,plasma and mercury,polymers,etc.Besides,one of the essential characteristics of fluids under investigation herein is the pressure tensor,which is associated with the fee of deformation tensor utilizing non-linear courting.Valuable contributions to the investigation of the properties of such fluids are highlighted by following experts[1–5].

The new development underwritten in fluid mechanics has paved the way for innovative methods in investigating the two-dimensional flow,where the non-Newtonian fluids’ flow mathematical system is far complex.Most known modules are the second grade of the non-Newtonian,which contains Maxwell,Oldroyd-B,as well as the power law of mathematical model.Eyring and Powell proposed a whole new mathematical model for non-Newtonian fluids in 1944 called the Eyring Powell model.Thus,very scant research has been involved in the fluid model of Eyring Powell on a two-dimensional type even though the model is better in several ways than other forms of non-Newtonian fluid models.Eyring Powell’s model has possessed numerous benefits for two main reasons.First,it is based on the theory of kinetic towards objects instead of that experimental formulas.Apart from that,the modeling of Eyring Powell reacts accordingly to low and high behavior of Newtonian’shear rates.

Despite facing the challenges of solving the Eyring Powell model due to its complexity,Jalil et al.[6]took the challenge of solving the non-Newtonian Eyring Powell problem passing a permeable surface that is continuously moving parallel with the free stream.Besides,an analytical solution for the effects of convective boundary conditions over a surface moving with convective boundary conditions under a constant free stream movement has been analyzed by Hayat et al.[7].Moreover,the ordinary equation arising from the partial governing differential upon applying proper transformation has been solved in a numerical manner using the Keller box method(KBM).Furthermore,Khader et al.[8]employed a numerical approach to investigate an unsteady flow and thermal conduction of Eyring Powell fluids in a laminar fluid-filled on a horizontal stretch surface under the condition of internal thermal generation.Likewise,Jalil et al.[9]and Javed et al.[10]devoted their study on fluids’ flow of the non-Newtonian Eyring Powell on the thin layer moving over a horizontal stretching sheet solved utilizing a finite difference approach.Moreover,Malik et al.[11]discussed a boundary layer of Eyring Powell for variable viscosity of a stretched cylinder.

Apart from that,Akbar et al.[12]deliberated the impact on free convection flow,mass and heat of viscous dissipation employing the numerical technique for computation.Also,Roşca et al.[13]discussed the flow and heat transfer of Powell Eyring fluids on a shrinking surface aligned with the free stream movement.Also,the effects of magnetohydrodynamics of Powell Eyring fluids with heat transfer were examined by Hayat et al.[14],Akbar et al.[15],Ellahi et al.[16],Hayat et al.[17],Gireesha et al.[18].Moreover,Abbasi et al.[19]and Hina[20]concentrated on investigating the peristaltic fluid of Eyring Powell.In contrast,Hayat et al.[21]and Rauf et al.[22]solved such fluid using optimal control embedded with double-diffusive Cattaneo–Christov heat and mass flux theories,respectively.The series comparing numerical solutions of Eyring-Powell fluid flow with NH and heat generation is presented by Hayat et al.[23]and Rehman et al.[24],separately.Furthermore,the study on Prandtl-Eyring fluids forced convection flow Eyring Powell and its theoretical similarity was carried out by Shukla et al.[25].Moreover,Sajid et al.[26]and Khan et al.[27]solved the fluid flow problem of Maxwell and Erying Powel models,respectively in the existence of nanoparticles with the effect of non-linear thermal radiation on a stretched surface.

Motivated by researchers examined above,the incompressible boundary layer of forced convection in two dimensions over a stretching sheet can be considered with the NH boundary condition as a heating circumstance.This recent study is dedicated to examining the forced convective of Eyring Powell fluid’s flow incorporated with TDV along a vertical stretching sheet through embedded NH boundary conditions.Numerical solutions are worked out by employing the KBM,and the solution is expressed through the graphical and tabular arrangement.

2 Problem Formulation

This study aims to investigate the characteristic of Eyring fluid under the influence of TDV embedded with the Newtonian heating(NH)boundary condition.The flow is assumed to flow over a vertical stretching sheet.It is presumed that the sheet stretches in the direction and moves with uniform velocity,where thexaxis is defined as an upward direction alongside the sheet.At the same time,they-axis assumes its upward direction concerning the surface.The geometric configuration of the study case is shown in Fig.1.

Figure 1:Physical configuration of problem

The derivation of the Eyring Powell model is based on rate processes theory to demonstrate the sheer of non-Newtonian flow.Summing up from the assumptions mentioned above and adopting the boundary layer approximation in the Cartesian coordinate system,the equations representing the model of fluid with TDV can be conveyed as

where(u,v),T,ρ,cp;and μ respectively refer to velocity components in(x,y)directions,temperature,density,constant pressure of specific heat,and fluid’s viscosity coefficient whilec*andare the Powell Eyring parameters’ fluid model which has(time-1)dimension.

The model is referred to conditions at the boundary as

In Eq.(4),the parametersk,hf,TfandT∞correspond to thermal conductivity,heat transference coefficient,hot fluid and ambient temperature,respectively.To obtain the set of similarity equations in the form of ODE,the similarity transformation variables as in Eq.(5)are adopted and applied to the governing Eqs.(1)to(4).

To simulate temperature-dependent viscosity variation,we adopt the Reynolds exponential viscosity model which provides an accurate approach,given by:

where the inverse linear function of temperature is assumed for variation of fluid viscosity.The resulting equations are obtained as follows:

together with the reformed boundary conditions given by:

3 Computation Procedures

The obtained equations under the arrangement of a set of ODE is solved by employing the Keller box technique.Eqs.(8)and(9)subjected to the boundary conditions(10)are reduced to the system of first order.For that matter,the independent variables are demonstrated as follows

Then,the respective equations can be written as

In Eqs.(7)to(12),the prime notation(′)corresponds to the differentiation with respect to η.Additionally,fluid parameters,MandB,Prandtl number,Pr,specific heat ratio of mixture,γ,thermal diffusivity,α and fluid kinematic viscosity,νfcan be defined as follows

The net rectangle in the η level is considered as in Fig.2.The points of net are denoted as

wherehjis the Δηj-spacing andjis the order of numbers specifying the coordinate location.

Figure 2:Net rectangle for difference approximations

Eqs.(10)–(12)are predicted by expending the central difference at mid-pointof the segmentP1P2.Thus,the following set of equations are obtained

Eqs.(15)–(19)are computed forj=1,2,…,Jat the givennwhere the boundary conditions(9)become

The matrixAis acknowledged as a tridiagonal matrix having zero elements except for its main diagonal.The respective system can be elucidated by using a block elimination technique with the assumption that the matrixAis non-singular and able to be factorized in the form of

4 Validation Procedure

The numerical results for this current investigation was computed using KBM.The outcome of fluid parameters(MandB),Prandtl number(Pr),and thermal diffusivity(α)on fluids’ velocity and temperature against η are computed using MATLAB software.It is well known that this method is undeniably one of the suitable approaches for solving the flow problem based on its wide application in many previous and recent studies within this research area.By having the finite boundary layer thickness,η∞=8,the boundary conditions of this study are fully satisfied based on both velocity and temperature profiles attaining the asymptotic behavior.It is important to mention that the exact solution for skin friction coefficient for Eq.(12)without the parameter representing Eyring Powell fluid and viscosity is given by

The exact expression of Eq.(41)has also been documented by Salleh et al.[28].It is worth mentioning here that the comparative outcomes with the exact solution are necessary to claim the current model and its output are acceptable.The physical measures of skin friction as well as Nusselt number for this present mathematical model are given as

Besides that,the present model can be reduced to the existing equation,as in the literature.The following summaries are presented to show its connection to the model.

5 Results and Discussion

This present paper concentrates on the discussion of the solutions obtained(tabular and graphical form).The results were computed under the NH boundary condition for various parameter α,M,B,and Pr.The computation is performed by assigning a set of fixed values of the parameter,where the value ofandare later computed one by one.The direct comparison study with the exact expression(46)as well as the existing study reported by Salleh et al.[28],Gorla et al.[29],and Vajravelu et al.[30]carried out to corroborate the numerical solutions acquired in this study.From Tabs.1 and 2,an excellent agreement is achieved,which indicates that the current model and its findings are acceptable.

Table 1:Comparative study on f ′′(0)

Table 2:Comparative study on value θ(0)and -θ′(0)

Tab.3 depicts the variations ofCfRex1/2andNuxRex-1/2with the variation of the values of the parameter involved.For the entire computation,the value γ was taken as 1.It is noticed that for larger values of α andB,the skin friction coefficient increases while decreasing trend is observed for larger values ofMand Pr.On the other hand,the variant of Nusselt number increases only for larger values of α,while a decreasing trend for Pr,MandBis observed.The increment in skin friction is due to the enhancement of drag forces occurring on the surface,and the increase in Nusselt number attributes to the store of heat from the process of fluid movement.It can be noticed that all the parameters existing in the present model affects the fluid flow characteristics.

Table 3:Numerical results of Cf Rex1/2 and NuxRex-1/2 for numerous values of Pr,α,M and B

Furthermore,the analysis on parameters ofM,B,Pr and α on velocity and temperature of fluid,respectively,are described in Figs.3 to 10.An increased behavior of velocity can be seen in larger values of Pr andB,as captured in Figs.5 and 7,respectively.This is because the properties of those parameters characterize the ability to exchange thermal energy with its surroundings,which reduces the viscosity of fluid particles.Meanwhile,a contradict behavior on velocities of fluid can be seen for larger values ofMand α,as illustrated in Figs.3 and 9,correspondingly.The temperature distribution shows a decreasing trend for larger values ofM,Band Pr whereas an increasing trend for larger α is obtained,as captured in Figs.4,6,8 and 10,respectively.This is due to lessen heat loss and viscosity.Besides that,all the figures indicate that the velocity and temperature of the fluid asymptotically approaches zero,which fittingly fulfils the boundary condition.

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

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

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

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

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

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

Figure 9: f ′(η)at M =0,B=1 and Pr=10 for various values of α

Figure 10:θ(η)at M =0,B=1 and Pr=10 for various values of α

6 Concluding Remark

In this research,the convective force of an Eyring Powell fluid along vertical stretch plate has been discussed by highlighting the effect of TDV involved the parameters of α,M,B,and Pr.A similar trend can be noticed from the mathematical analysis in the motion and temperature distributions of fluid,respectively,when parameters are increased.The utilization of KBM has developed a numerical solution,where the given graphs and tables demonstrate the effects of parameters.Nevertheless,the deviations on the distribution of velocity and the temperature of fluid exhibit the significance of the parameters studied.The findings in this study are expected to contribute to a better understanding of single-phase fluid flow characteristics as much as the solutions of its flow problems.

Funding Statement:This project has been supported by Universiti Malaysia Pahang&Ministry of Higher Education under The Fundamental Research Grant Scheme for Research Acculturation of Early Career Researchers(FRGS-RACER)(Ref:RACER/1/2019/STG06/UMP//1)through RDU192602.

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