Mixed Convection Flow of Magnetic Viscoelastic Polymer from a Nonisothermal Wedge with Biot Number Effects

Magnetic polymers are finding increasing applications in diverse fields of chemical and mechanical engineering. In this paper, we investigate the nonlinear steady boundary layer flow and heat transfer of such fluids from a nonisothermal wedge. The incompressible Eyring-Powell non-Newtonian fluid model is employed and a magnetohydrodynamic body force is included in the simulation.The transformed conservation equations are solved numerically subject to physically appropriate boundary conditions using a second-order accurate implicit finite difference Keller Box technique.The numerical code is validated with previous studies. The influence of a number of emerging nondimensional parameters, namely, the Eyring-Powell rheological fluid parameter (ε), local non-Newtonian parameter based on length scale (δ), Prandtl number (Pr), Biot number (γ), pressure gradient parameter (m), magnetic parameter (M), mixed convection parameter (λ), and dimensionless tangential coordinate (ξ), on velocity and temperature evolution in the boundary layer regime is examined in detail. Furthermore, the effects of these parameters on surface heat transfer rate and local skin friction are also investigated.


Introduction
The development of modern functional materials which can be manipulated using electromagnetic fields has stimulated great attention in polymer engineering sciences in recent years.These materials include ferromagnetic polymers [1], magnetotropic fluids [2], and liquid crystalline electrically conducting polymers [3].Such fluids exhibit many complex characteristics including non-Newtonian behavior.For better understanding of the manufacture of such materials, magnetorheological fluid mechanics plays a central role.Although many different magnetic properties may arise in such materials, a fundamental methodology for their synthesis involves the application of transverse magnetic fields which alters momentum transfer and therefore also influences coupled processes such as heat transfer and mass transfer (species diffusion).Many constitutive material models have been developed to simulate the departure of such fluids from Newtonian viscous behavior.In these fluids, the constitutive relationship between stress and rate of strain is nonlinear in comparison to the Navier-Stokes equations which are generally good for Newtonian fluids.Most non-Newtonian models involve some form of modification to the momentum conservation equations.A comprehensive summary of such models which include the Maxwell model, power-law model, and also Eyring-Powell and Giesekus viscoelastic models, is provided in the treatise by Shaw [4].Heat transfer is often also frequently present in polymer dynamics [5].Furthermore, in many polymer fabrication flows, boundary layer phenomena arise.The convective boundary condition has therefore also attracted some interest and this is usually simulated via a Biot number in the wall thermal boundary condition.Recently, Ishak [6] discussed the similarity solutions for flow and heat transfer over a permeable surface with convective boundary condition.Aziz [7] provided a similarity solution for laminar thermal boundary layer over a flat surface with a convective surface boundary condition.Aziz [8] further studied hydrodynamic and thermal slip flow boundary layers with an isoflux thermal boundary condition.Buoyancy effects 2 International Journal of Engineering Mathematics on thermal boundary layers subject to a convective surface boundary condition were examined by Makinde and Olanrewaju [9].Gupta et al. [10] used a variational finite element to simulate mixed convective-radiative non-Newtonian shrinking sheet flow with a convective boundary condition and Eringen's micropolar material model.Swapna et al. [11] studied convective wall heating effects on hydromagnetic flow of a micropolar fluid.Makinde et al. [12] studied cross diffusion effects and Biot number influence on hydromagnetic Newtonian boundary layer flow with homogenous chemical reactions and MAPLE quadrature routines.Bég et al. [13] analyzed Biot number and buoyancy effects on magnetohydrodynamic thermal slip flows.Subhashini et al. [14] studied wall transpiration and cross diffusion effects on free convection boundary layers with a convective boundary condition.In many of these studies, magnetohydrodynamic effects were also considered generally via the introduction of a Lorentzian magnetic drag force and consideration of the Hartmann number.Relatively few studies however have considered magnetic field effects on Eyring-Powell viscoelastic polymer flows.This rheological model has certain advantages over the other non-Newtonian formulations, including simplicity, ease of computation, and physical robustness.Furthermore, it is deduced from kinetic theory of liquids rather than the empirical relation.Additionally, it correctly reduces to Newtonian behavior for low and high shear rates [15].It is a three-constant model which displays a nonzero bounded viscosity at both the upper and the lower limits.Distinct from the popular Ostwald-De Waele power-law model, the Eyring-Powell model does not demonstrate infinite effective viscosities for low shear rates.Investigations employing this model of relevance to polymeric transport processes include Hayat et al. [16] in radiative magnetic convection, Sirohi et al. [17] for wedge flows, and Adesanya and Gbadeyan [18] for channel flows.Very recently, Prasad et al. [19] studied the thermal convection flow in permeable materials with Biot number effects, observing that increasing rheological effect accelerates the flow and heats the boundary layer.
Very few of the above studies have considered Falkner-Skan flows [20].This family of boundary layer flows is associated with the two-dimensional wedge configuration.Non-Newtonian flows from wedge bodies arise in a number of chemical engineering systems which have been described in detail by Peddieson [21] employing the second-order Reiner-Rivlin model.The mixed convection boundary layer flow from a heated wedge plate has also drawn some interest.The combined forced and free convection flow and heat transfer about a nonisothermal wedge subject to a nonuniform free stream velocity were first considered by Sparrow et al. [22].Watanabe et al. [23] analyzed theoretically mixed convection flow over a perforated wedge with uniform suction or injection.Kafoussias and Nanousis [24] and Nanousis [25] studied the effect of suction or injection on MHD mixed convection flow past a wedge.Gorla [26] used a power-law model to study heat transfer in polymer flow past a wedge.Yih [27] evaluated radiation effects on mixed convection flow about an isothermal wedge embedded in a saturated porous medium.Rashidi et al. [28] developed homotopy solutions for third grade viscoelastic flow from a nonisothermal wedge.
Chamkha et al. [29] presented computational solutions for MHD forced convection flow from a nonisothermal wedge in the presence of a heat source or sink with a finite difference method.Hsiao [30] reported on MHD convection of viscoelastic fluid past a porous wedge, observing that the elastic effect increases the local heat transfer coefficient and heat transfer rates at the wedge surface.Ishak et al. [31] obtained a self-similar solution for a moving wedge in a micropolar fluid.Ishak et al. [32] further studied numerically steady two-dimensional laminar flow past a moving wedge in non-Newtonian fluid.
The objective of the present study is to investigate the laminar boundary layer flow and heat transfer of an Eyring-Powell non-Newtonian fluid from a nonisothermal wedge.Such a study has not appeared in the literature to the knowledge of the authors.The nondimensional equations with associated dimensionless boundary conditions are solved with the Keller implicit finite difference "Box" scheme [33].The effects of the emerging thermophysical parameters, namely, the rheological parameters (, ), Biot number (), mixed convection parameter (), pressure gradient parameter (), magnetic parameter (), and Prandtl number (Pr), on the velocity, temperature, local skin friction, and heat transfer rate (local Nusselt number) characteristics are studied.The present problem is relevant to the simulation of magnetized polymer materials fabrication processes.

Non-Newtonian Constitutive Eyring-Powell Fluid Model
In the present study, a subclass of non-Newtonian fluids known as the Eyring-Powell fluid is employed.The Cauchy stress tensor for Eyring-Powell fluids [15] takes the hyperbolic form: where  is dynamic viscosity and  and  are the rheological fluid parameters of the Eyring-Powell fluid model.Consider the second-order approximation of the sinh −1 function as follows: When the Eyring-Powell formulation is introduced into the subsequent model, the momentum conservation equation is significantly modified with numerous velocity gradient terms including mixed derivatives.Strong nonlinearity therefore results in necessitating numerical solutions.The resulting boundary value problem is found to be well posed and permits an excellent mechanism for the assessment of rheological characteristics on the flow behaviour.

Mathematical Thermofluid Polymer Boundary Layer Model
Steady, laminar, and incompressible flow and heat transfer of an Eyring-Powell polymeric fluid from a nonisothermal horizontal wedge are considered, as illustrated in Figure 1.
The -coordinate (streamwise) is measured along the wall of the wedge and -coordinate (transverse) is directed normal to it.A uniform magnetic field of strength  0 is applied parallel to the -axis, that is, transverse to the wedge surface.Magnetic induction, Hall current, ion slip, and Maxwell displacement currents are neglected.Magnetic Reynolds number is small and the electric field due to polarization of charges is negligible.The gravitational acceleration  acts downwards.We also assume that the Boussinesq approximation holds; that is, density variation is only experienced in the buoyancy term in the momentum equation.Initially, both the wedge and Eyring-Powell fluid are maintained at the same temperature.Instantaneously, they are raised to a temperature   >  ∞ , the ambient temperature of the fluid which remains unchanged.In line with the approach of Yih [27] and introducing the boundary layer approximations, the equations for mass, momentum, and energy conservation can be written as follows: Here,  and V are the velocity components in and -directions, respectively, ] = / is the kinematic viscosity of Eyring-Powell fluid (polymer),  1 is the coefficient of thermal expansion,  is the thermal diffusivity,  is the temperature, and  is the density of the fluid.The Eyring-Powell fluid model therefore introduces a mixed derivative (second order, second degree) into the momentum boundary layer equation ( 4).The non-Newtonian effects feature in the shear terms only of (4) and not the convective (acceleration) terms.
The fourth term on the right hand side of (4) represents the thermal buoyancy force on the flow field, with "+" and "−" signs referring to the buoyancy-assisting and buoyancyopposing flow scenarios, and couples the velocity field (4) with the temperature field equation ( 5).The physical boundary conditions are given as follows: At  = 0,  = 0, As  → ∞, Here,  ∞ is the free stream temperature,  is the thermal conductivity, ℎ  is the convective heat transfer coefficient,   is the convective fluid temperature,  ∞ is the free stream velocity, and  =  2 /(2 −  2 ) is the Hartree pressure gradient parameter which corresponds to  2 = Ω/ for a total angle Ω of the wedge. 1 is a positive number.The stream function  is defined by  = / and V = −/, and, therefore, the continuity equation is automatically satisfied.In order to render the governing equations and the boundary conditions in dimensionless form, the following nondimensional quantities are introduced: All terms are defined in Nomenclature.In view of the transformation defined in (7), the boundary layer equations ( 4)-( 6) are reduced to the following coupled, nonlinear, and dimensionless partial differential equations for momentum and energy for the domain: The transformed dimensionless boundary conditions are as follows: At  = 0,  = 0, As  → ∞, Here, primes denote the differentiation with respect to  and is the Biot number.The wall thermal boundary condition in (10) corresponds to convective cooling.The skin friction coefficient (shear stress at the cylinder surface) and Nusselt number (heat transfer rate) can be defined using the transformations described above with the following expressions: The location,  ∼ 0, corresponds to the vicinity of the lower stagnation point on the wedge.For this scenario, the model defined by ( 8) to ( 9) contracts to an ordinary differential boundary value problem: The general model is solved using a powerful and unconditionally stable finite difference technique introduced by Keller [33].The Keller Box method has a second-order accuracy with arbitrary spacing and attractive extrapolation features.
It remains one of the most reliable numerical methods for nonlinear boundary layer flows.

Numerical Solution with Keller Box Implicit Method
The Keller Box implicit difference method is implemented to solve the nonlinear boundary value problem defined by ( 8) and ( 9) with boundary conditions (10).This technique, despite recent developments in other numerical methods, remains a powerful and very accurate approach for boundary layer flow equation systems which are generally parabolic in nature.It is unconditionally stable and achieves exceptional accuracy.An excellent summary of this technique is given in Keller [33].Magnetohydrodynamics applications of Keller's method are reviewed in Anwar Bég [34].This method has also been applied successfully in many rheological flow problems in recent years.These include oblique micropolar stagnation flows [35], Walter's B viscoelastic flows [36], Stokesian couple stress flows [37], hyperbolic-tangent convection flows from curved bodies [38], micropolar nanofluids [39], Jeffreys elastoviscous boundary layers [40], magnetic Williamson fluids [41], and Maxwell fluids [42].The Keller Box discretization is fully coupled at each step which reflects the physics of parabolic systems, which are also fully coupled.Discrete calculus associated with the Keller Box scheme has also been shown to be fundamentally different from all other mimetic (physics capturing) numerical methods, as elaborated in Anwar Bég [34].The Keller Box scheme comprises four stages: (1) Decomposition of the th order partial differential equation system to  first-order equations.
(4) Block-tridiagonal elimination solution of the Linearized Keller Algebraic Equations.
The algebraic details for the present problem are provided in Appendix A. A typical mesh is given in Figure 2.

Validation with Nakamura Tridiagonal Method (NTM)
To verify the accuracy of the Keller Box method (KBM) code, we have also solved the seventh-order nonlinear partial differential boundary value problem defined by the parabolic equations ( 8) and ( 9) with boundary conditions (10), with an alternative finite difference procedure, namely, the Nakamura tridiagonal method (NTM) [43].NTM is also second order, accurate, and very easily implemented in boundary layer problems.It has proven to be particularly adept at resolving nonlinear and mixed derivatives which characterize non-Newtonian (e.g., viscoelastic) problems.Applications of NTM in rheological boundary layer flows include secondorder Reiner-Rivlin wedge convection in porous media [44], power-law convective flows in porous media [45], micropolar wall plume heat transfer [46], and nanofluid bioconvection in permeable materials [47] and magnetohydrodynamic micropolar enrobing flows [48].Details of the technique with solutions are provided in Appendix B and Tables 6 and 7.

Discussion and Interpretation
In Figures 3-11, we note that asymptotically smooth profiles are always achieved in the free stream.An adequate value of infinity has therefore been prescribed in the computations.Figures 3(a 4(a), although there is a very slight increase in velocity with increase in .Similarly, there is only a very slight depression in temperature magnitudes in Figure 4(b) with a rise in .This parameter is defined as ] 2 Re 3  /(2 2  4 ) and features in a single negative term in (9), namely, −(  ) 2   , unlike the  rheological parameter which also appears in the shear term, −(1 + ) 2   .This parameter also features viscosity which is effectively decreased weakly and induces a slight acceleration in the flow.The momentum boundary layer thickness is therefore slightly decreased whereas since heat diffusion is decreased, thermal boundary layer thickness is therefore marginally reduced.
Figures 5(a) and 5(b) depict the evolution of velocity (  ) and temperature () functions with a variation in Biot number, .Dimensionless velocity component (Figure 5(a)) is considerably enhanced with increasing .In Figure 5(b), an increase in Biot number is seen to considerably enhance temperatures throughout the boundary layer regime.For  < 1, that is, small Biot numbers, the regime is frequently   designated as being "thermally simple" and there is a presence of more uniform temperature fields inside the boundary layer and the wedge solid surface.For  > 1, thermal fields are anticipated to be nonuniform within the solid body.The Biot number effectively furnishes a mechanism for comparing the conduction resistance within a solid body to the convection resistance external to that body (offered by the surrounding fluid) for heat transfer.We also note that a Biot number in excess of 0.1, as studied in Figures 5(a) and 5(b), corresponds to a "thermally thick" substance whereas Biot number less than 0.1 implies a "thermally thin" material.
Since  is inversely proportional to thermal conductivity (), as  increases, thermal conductivity will be reduced at   be sustained at the surface, as witnessed in Figure 5(b).For a fixed wall convection coefficient and thermal conductivity, Biot number as defined in  = (ℎ  /)Re −1/2  is also inversely proportional to the local Reynolds (free convection) number.However, as the Biot number increases, the local Reynolds number must decrease and this will induce the opposite effect, that is, accelerating the boundary layer flow, as shown in Figure 5(a).
Figures 6(a) and 6(b) depict the profiles for velocity (  ) and temperature () for various values of pressure gradient parameter, .It is observed that increasing  significantly accelerates the flow; that is, velocity increases.This trend is familiar from other studies, for example, Nanousis [25] and Gorla [26].However, increasing  is found to decelerate the temperature.Effectively, with greater pressure gradient effect, the momentum boundary layer thickness is reduced and thermal boundary layer thickness is decreased.In (9), with  = 0, the wedge flow becomes Blasius flow from a flat plate.With  = 1, we retrieve the case of flow in the vicinity of stagnation point flow on an infinite plate.While these two flow cases are of mathematical interest, in practical polymer processing, they are not.The generalized wedge case is more relevant for which  ̸ = 0, as studied in Figures 6(a) and 6(b).Figures 7(a) and 7(b) depict the profiles for velocity (  ) and temperature () for various values of mixed convection parameter, .It is found that an increase in  increases the velocity.However, the temperature decreases significantly.The mixed convection parameter  = Gr  /Re 2  couples the momentum field with the energy field, via the buoyancy term,  sin(Ω/2).Although this is a first-order term, it has a profound influence on the interplay of velocity and temperature.With positive values, we have buoyancy-aided flow and with negative values we have buoyancy-opposed flow.Both are studied in Figures 7(a) and 7(b).Evidently when buoyancy aids the flow ( > 0), velocities are increased, with the contrary effect for buoyancy opposition.The momentum boost achieved with assistive buoyancy force results however in a depletion of thermal energy which cools the boundary layer and decreases temperatures (Figure 7(b)).The reverse effect is generated with opposing buoyancy force which heats the boundary layer and increases thermal boundary layer thickness.
Figures 8(a) and 8(b) depict the profiles for velocity (  ) and temperature () for various values of Prandtl number, Pr.It is observed that an increase in the Prandtl number significantly decelerates the flow; that is, velocity decreases.Also increasing Prandtl number is found to decelerate the temperature.Prandtl numbers of 100 are more representative of weak polymer solutions and are associated with slower and cooler flows.Momentum boundary layer thickness is greater for such fluids as is thermal boundary layer thickness.
Figures 9(a) and 9(b) depict the profiles for velocity (  ) and temperature () for various values of magnetic parameter, .It is observed that an increase in  significantly decelerates the flow; that is, velocity decreases.Conversely, increasing  is found to enhance the temperature.The parameter, , is a Hartmann number.It simulates the relative contribution of Lorentzian magnetohydrodynamic drag force relative to viscous hydrodynamic force.As  is increased, greater opposition is generated to the flow past the wedge leading to deceleration.The supplementary work expended in dragging the polymer against the imposition of the transverse magnetic field creates heating in the polymer.This dissipation of heat leads to a temperature rise and thickening of thermal boundary layers in polymers.Such phenomena are documented extensively in magnetohydrodynamic studies, for example, in Makinde et al. [12] and Chamkha et al. [29].
Figures 10(a) and 10(b) show the influence of Eyring-Powell fluid parameter, , on dimensionless skin friction coefficient ((1 + )  (, 0) − (/3) 3 (  (, 0)) 3 ) and heat transfer rate (−  (, 0)) at the wedge surface.It is observed that the dimensionless skin friction is increased with the increase in ; that is, the boundary layer flow is accelerated with decreasing viscosity effects in the non-Newtonian regime.Conversely, the surface heat transfer rate is substantially decreased with increasing  values.Decreasing viscosity of the fluid (induced by increasing the  value) reduces thermal diffusion as compared with momentum diffusion.A decrease in heat transfer rate at the wall will imply that less heat is convected from the fluid regime to the wedge, thereby heating the boundary layer and enhancing temperatures.
Figures 11(a In Table 1, we present the influence of the Eyring-Powell fluid parameter, , on the skin friction and heat transfer rate, along with a variation in Pr and .With increasing , the skin friction is enhanced.The parameter  is inversely proportional to the dynamic viscosity of the non-Newtonian fluid.There as  is elevated, viscosity will be reduced and this will induce lower resistance to the flow at the surface of the wedge, that is, accelerating the flow leading to an escalation of shear stress.Furthermore, this trend is sustained at any Pr.However, an increase in Pr markedly reduces the shear stress magnitudes.Similarly increasing  is observed to reduce heat transfer rates, again at all Pr's, whereas Pr strongly accentuates heat transfer rates.Magnitudes of shear stress are always positive indicating that flow reversal (backflow) never arises.It is also observed that increasing mixed convection parameter, , increases both skin friction and heat transfer rates.
Table 2 shows results for the influence of the local non-Newtonian parameter (based on length scale ), that is,  and the mixed convection parameter () along with the variation in  on skin friction and heat transfer rate.Skin friction is generally decreased with increasing .However, heat transfer rate is found to be enhanced with increasing . = ] 2 Re 3  /2 2  4 and inspection of this definition shows that the effectively inverse proportionality of  to kinematic viscosity (]) (with all other parameters being maintained constant) will generate a strong resistance to the flow leading to a deceleration, that is, drop in shear stresses.Conversely, the direct proportionality of  to Reynolds number (Re  ) will imply that thermal buoyancy forces are enhanced as  increases and this will cause a boost in heat transfer by convection from the wedge surface manifesting with the greater heat transfer rates observed in Table 2.This table also shows that, with an increase in , the skin friction and heat transfer rate are elevated.Table 3 presents the results for the influence of the pressure gradient parameter () and the magnetic parameter () along with the variation in  on skin friction and heat transfer rate.Skin friction and heat transfer rate are generally decreased with increasing .Flow deceleration is therefore confirmed with greater transverse magnetic field, associated with the escalation in Lorentzian magnetohydrodynamic impedance.However, skin friction and heat transfer rate are found to be enhanced with increasing .Table 4 presents the Keller Box numerical solutions of the missing condition   (, 0) (in brackets) and skin friction   for various values of  and .It is found that skin friction is reduced with increasing values of .Furthermore, the skin friction   is observed to be increased with a rise in the Eyring-Powell fluid parameter () for all values of the local non-Newtonian parameter ().
Table 5 presents the comparison results of Nu for different values of  which are found to be in good correlation with those of [49].

Conclusions
Numerical solutions have been presented for the nonsimilar, buoyancy-driven boundary layer flow and heat transfer of an Eyring-Powell fluid external to a nonisothermal wedge.This model has been developed to analyse magnetic polymer Falkner-Skan flows of relevance in chemical engineering processing.The Keller Box implicit second-order accurate finite difference numerical scheme has been utilized to efficiently solve the transformed, dimensionless velocity and thermal boundary layer equations, subject to realistic boundary conditions.A comprehensive assessment of the effects of Eyring-Powell fluid parameter (), local non-Newtonian parameter (), Biot number (), mixed convection parameter (), pressure gradient parameter (), magnetic body force parameter (), and Prandtl number (Pr) on thermofluid characteristics has been conducted.Significant flow deceleration has been demonstrated for increasing Eyring-Powell fluid parameter () and negative mixed convection parameter (), that is, buoyancy opposition and magnetic body force parameter () with generally greater temperatures observed with stronger magnetic field.Excellent correlation with Nakamura tridiagonal method computations has also been achieved, verifying the validity of the present Keller Box solutions.Blasius and stagnation point flow cases have been retrieved as special cases of the generalized wedge flow model.Generally very stable and accurate solutions are obtained with the present finite difference code.The numerical code is able to solve nonlinear boundary layer equations very efficiently and therefore shows excellent promise in simulating transport phenomena in other non-Newtonian fluids.It is therefore presently being employed to study micropolar fluids and viscoplastic fluids which also represent other chemical engineering working fluids.The present study has also neglected time effects.Future simulations will also address transient polymeric boundary layer flows and will be presented soon.

A. Keller Box Numerical Details
As mentioned earlier, the algebraic details for the Keller Box method are documented here.

Figure 1 :
Figure 1: Magnetic polymeric flow from a wedge.
) and 3(b) illustrate the effect of Eyring-Powell fluid parameter  on the velocity (  ) and temperature distributions through the boundary layer regime.Velocity is significantly decreased with increasing  at larger distance from the wedge surface.Conversely, temperature is consistently enhanced with increasing values of .The mathematical model reduces to the Newtonian viscous flow model as  → 0 and  → 0. The momentum boundary layer equation in this case contracts to the familiar equation for Newtonian mixed convection from a plate; namely,   +((1+)/2)  + (1 − (  ) 2 ) +  sin(Ω/2) −   = (1 − )(  (  /) −   (/)).The thermal boundary layer equation (9) remains unchanged.In Figure 3(b), temperatures are clearly

8 𝜉 = 1 Figure 4 :
Figure 4: (a) Influence of  on the Velocity Profiles.(b) Influence of  on the Temperature Profiles.

Figure 5 : 4 7Figure 6 :
Figure 5: (a) Influence of  on the Velocity Profiles.(b) Effect of  on the Temperature Profiles.

Figure 7 :Figure 8 :
Figure 7: (a) Influence of  on the Velocity Profiles.(b) Influence of  on the Temperature Profiles.

Figure 9 :Figure 10 :
Figure 9: (a) Influence of  on the Velocity Profiles.(b) Influence of  on the Temperature Profiles.

Figure 11 :
Figure 11: (a) Influence of  on skin friction coefficient.(b) Influence of  on Nusselt number.
minimized for the Newtonian case ( = 0) and maximized for the strongest non-Newtonian case ( = 1.0).The non-Newtonian parameter  (= 1/[]) is evidently inversely proportional to the viscosity and also to  and  (Eyring-Powell rheological fluid parameters).As  is increased, the material viscosity therefore will decrease which would aid momentum development.Similarly, the parameters  and  are simultaneously reduced (they are strongly linked to viscosity) and the overall effect is a boost in material viscosity.This results in reduction in the  magnitude which manifests in marked deceleration in the flow as shown in Figure3(a).With decreased momentum diffusion rates, the velocity boundary layer thickness is decreased.Thermal diffusion is conversely enhanced and the fluid is heated (Figure3(b)), resulting in a thickening in the thermal boundary layer.Figures 4(a) and 4(b) depict the velocity (  ) and temperature () distributions with increasing local non-Newtonian parameter, .Very little tangible effect is observed in Figure

Table 6 :
Values of and Nu for different values of , , and  (Pr = 7.0,  = 0.1,  = 0.1,  = 0.1,  = 0.3, and Ω = 30 ∘ ).the cylinder surface and this will lead to a decrease in the rate of heat transfer from the boundary layer to within the cylinder, manifesting in a rise in temperature at the cylinder surface and in the body of the fluid; the maximum effect will International Journal of Engineering Mathematics ).