Effects of Variable Viscosity in Unsteady Magnetohydrodynamic Hybrid Nanoﬂuid Flow over Stretching/Shrinking Cylinder with Partial Slip and Stefan Blowing

: In the manufacturing sector, transport phenomena near the stagnation region are frequent, particularly in the polymer and extrusion processes, which require continuous improvement to raise the process’s quality standards. The aim of this study is to explore the improvement of heat and mass transmission using unsteady magnetohydrodynamic (MHD) hybrid nanoﬂuid (HNF) ﬂow over a stretching/shrinking cylinder with variable viscosity and Stefan blowing. The governed equations of heat and mass transfer processes are converted into ordinary differential equations (ODEs) using the appropriate transformations, and the resulting equations are then solved using the MATLAB package bvp4c . With an upsurge in the volume fraction of nanoparticles, the skin friction increases, but the reverse trend is detected with negative values for the unsteadiness constraint. The use of 2D graphs to show how important parameters affect the velocity, temperature, and concentration is thoroughly discussed. There is a discussion of the quantitative ﬁndings from the wall shear factor and the heat and mass transfer rates calculated for the stretching/shrinking cases.


Introduction
Scientists, researchers, and engineers are attentive to the boundary layer flow, heat transfer, and mass transfer of Newtonian fluids over a stretching/shrinking cylinder due to their many applications, including hardening, mining of metals, swelling processes, pipe manufacturing, wire diminishing, etc.Crane [1] highlighted the boundary layer flow on the stretching plate.Anuar et al. [2] explored heat transfer past an unstable stretchy sheet.Kumbhakar and Nandi [3] explored an HNF unsteady MHD radiative flow over a stretchy sheet with slip.Khashi et al. [4] examined the axisymmetric flow and heat transmission of an HNF flow past a shrinking/shrinking disc.Sreedevi et al. [5] conducted an investigation into the unsteady MHD heat and mass transfer of HNF flow through a stretching sheet and found that the temperature profile intensified as the tiny (nano) particle volume fraction rose.Raju et al. [6] conducted research on the nonlinear motions of axisymmetric ternary hybrid nanofluids in thermally radiated Darcy Walls with various shapes.Zainal et al. [7] investigated the unsteady MHD flow towards a horizontal cylinder with two different nanoparticles and found that the unsteadiness parameter improves the thermal performance.Zainal et al. [8] scrutinized the stability of an unstable HNF flow passing a permeable stretching/shrinking cylinder.Raju et al. [9] examined the role of body acceleration and gold nanoparticles in an uneven/composite inclined stenosis artery's blood flow.Raju et al. [10] studied the stagnant power-law fluid flow over a stretching spinning disc containing homogeneous-heterogeneous chemical reactions.Hussain et al. [11] conducted a heat transport investigation of HNF under the thermal radiation regime.
Applications of magnetohydrodynamics (MHD) include energy conversion, lowemission power generation systems, high-efficiency MHD accelerations and thrusters, etc. [12][13][14].Ashwinkumar et al. [15] scrutinized HNF flow over a convectively heated cylinder with MHD with different geometries.The effect of a stagnant flow of HNF with MHD was scrutinized by Mabood et al. [16].Liu et al. [17] discovered the slip flow of an unsteady nanofluid flow past a stretching plate.An investigation into an unstable MHD nanofluid stretching cylinder with a heat cause and sink was conducted by Song et al. [18].Waini et al. [19] studied unsteady HNF flow over a permeable shrinking/stretching surface and found that skin friction decreases but heat transmission rises with an increasing volume fraction of copper.Bilal et al. [20] inspected unsteady HNF flow through a horizontal channel.An unsteady stagnant flow of nanofluid stretching/shrinking over an outward surface with partial slip was investigated by Rajput et al. [21].Related work was carried out by Wang [22] and Waini et al. [23].The study of the stagnation-point flow of a viscous or non-Newtonian fluid is important in several engineering and industrial applications.Knowledge of the stagnation flow assists in the design of thrust bearings and radial diffusers, drag reduction, and transpiration cooling, as well as thermal oil recovery.
The assumption made in each of the aforementioned studies was that the fluid viscosity would remain constant.However, it is well known that temperature changes can significantly alter the fluid's physical characteristics.The viscosity change with temperature variation must be taken into account, to correctly predict the flow performance.Numerous recent studies have looked into how heat transfer past a moving sheet is impacted by fluid viscosity and thermal conductivity variations.The effect of variable viscosity on the natural convection of HNF flow because of a stretching surface was deliberated by Tulu et al. [24].The temperature-dependent viscosity variation in Maxwell nanofluid was predicted by Mabood et al. [25].
The species (concentration) field and the Stefan blowing effect are related.On an impermeable surface, the Stefan blowing effect appears.By stretching a cylinder to examine the HNF flow under Stefan blowing, Rana et al. [26] discovered that the increasing nanoparticle size and Stefan blowing have extreme effects on skin friction.Nanofluid flow towards a non-liner stretchy sheet with slip on unsteady MHD was studied by Jagan et al. [27].Hassen [28] looked at the effects of Stefan suction and blowing on an HNF based on ethylene glycol and water that exhibits partial slip and nanoparticle shape effects.Kiran et al. [29] examined the thermal transport analysis of non-uniformly heated annular geometry unsteady buoyant convective flow.In recent years, the nanofluids are used as an emerging working medium for fluid flow and heat transfer under various conditions [30][31][32].
Velocity slip is a phenomenon that occurs when fluids do not adhere to the surface of the body.One of the frequently occurring boundary conditions at which the surface velocity linearly falls on the fluid shear stress is Navier's partial slip.There are numerous fields where the slip state at the boundary surface is used, including lubrication, medicine, polishing artificial heart valves, etc.The HNF flow and heat transfer past a stretchy surface were researched by Kumar [33].MHD nanofluid flow with the effect of partial slip was scrutinized by Iskander et al. [34].Nadeem et al. [35] considered the slip effect in a cylinder of micropolar HNF.
The current research focuses on an HNF flow along a stretching/shrinking cylinder with an unsteady stagnation point that is subject to Stefan blowing and variable viscosity.In addition, the impact of the partial slip on the surface and magnetic field are taken into account.With the aid of the MATLAB bvp4c package, the resulting nonlinear dimensionless ODEs were solved, and the explanations were determined in terms of the velocity, friction drag, concentration, temperature, heat and mass transmission of the HNF, all of which were identified while they were constrained by the associated flow conditions.

Mathematical Formulation
Consider an incompressible, unsteady stagnation-point HNF flow past a stretching/shrinking cylinder with radius a = 1.Here, the magnetic field is applied vertically to the fluid flow path.The cylinder free stream velocity and surface are denoted by u e (x) = c 1 x L(1−c 3 t) and u w (x) = c 2 x L(1−c 3 t) .When c 2 = 0, the cylinder is static; when c 2 > 0, the cylinder is stretched; when c 2 < 0, the cylinder is shrunk; c 1 , c 2 , and c 3 are factors, and the characteristic length is L. The surface temperature, upstream temperature, surface concentration, and free stream concentration are symbolized as T w , T ∞ , C w , and C ∞ .Stefan blowing, stagnation point, variable viscosity, unsteady MHD, and partial slip effects are included as shown in Figure 1.When we start the stretch/shrink cylinder, the fluid starts flowing in x-direction but, in our model, we consider only the positive x-direction.The Stefan blowing effect applied is vertical to the fluid flow direction.Since water is used as the base fluid in our model, the Schmidt number is fixed.
account.With the aid of the MATLAB bvp4c package, the resulting nonlinear dime less ODEs were solved, and the explanations were determined in terms of the ve friction drag, concentration, temperature, heat and mass transmission of the HNF which were identified while they were constrained by the associated flow conditio

Mathematical Formulation
Consider an incompressible, unsteady stagnation-point HNF flow past a s ing/shrinking cylinder with radius a = 1.Here, the magnetic field is applied vertic the fluid flow path.The cylinder free stream velocity and surface are denoted by   From this model, the governing equations (referring to Zainal [7] and Wain [19]) are described as Continuity Equation From this model, the governing equations (referring to Zainal [7] and Waini et al. [19]) are described as Continuity Equation ∂ ∂x Temperature Equation The physical realistic boundary conditions are where the velocity mechanisms are u and v in the directions x and r.T represents the temperature of HNF.Moreover, lists of the thermophysical features of the base fluid and nanoparticles and the physical attributes of the HNF are described by Waini et al. [19,23].
The temperature-dependent variable viscosity is taken as B > 0 denotes an accelerating flow, while B < 0 denotes a reversing flow.In this study, we examine the fluid flow features and locate the point of separation for a decelerating stretching/shrinking cylinder (B < 0).Here, the nanoparticle volume fraction of Al 2 O 3 and Cu are designated by ϕ 1 and ϕ 2 .The nanoparticle volume fraction HNF is defined in Equation ( 6): The suitable transformations are when Equation ( 7) is substituted into (1), ( 1) is satisfied exactly.The following equations result from substituting (7) into Equations ( 2)-( 5). (1 with the boundary conditions: where Meanwhile, the coefficient of skin friction (C f ), Nusselt number (Nu), and Sherwood number (Sh) (refer to Waini et al. [23]) are expressed as follows.
whereτ w = µ hn f ∂w ∂r r=a , q w = −k hn f ∂T ∂r r=a and q m = −D ∂C ∂r r=a (13) Through the adoption of the similarity transformations Equations ( 7) and ( 13) by Equation ( 12), the following Equation is obtained.

Numerical Methods
The simplified form of the mathematical model boundary value problem solver MAT-LAB (version 7.50) bvp4c package is utilized to solve Equations ( 8)-( 11) mathematically, as designated by Waini et al. [23], and the following equations are obtained.

Results and Discussion
Different combinations of significant parameters are discussed in this study.Alumina Al 2 O 3 (ϕ 1 ) and copper Cu (ϕ 2 ) from 0 to 0.02 (2%), respectively, are used in the nanoparticle volume fraction.In Table 1, the current fallouts of f (1) are compared with those of Wang [22] and Waini et al. [23].In Table 2, the mathematical values of C f , Nu, and Sh are obtainable for different values of Sb, M, θ r , A, and ε, and it is numerically found that the skin friction, heat transfer, and mass diffusion increase with a rise in the variable viscosity parameter.From Figure 2a, increasing the Stefan blowing parameter velocity profile increases for the stretching case but decreases for the shrinking case.We can infer that, physically, mass diffusion causes a slowing of HNF as it approaches the surface due to an increase in Stefan blowing from the surface to the free stream when the cylinder is shrinking, but the opposite trend is obtained in the stretching case.Figure 2b,c show how the rising Stefan blowing parameter profile for the temperature and concentration are enhanced for both the stretching and the shrinking case.The presence of nanoparticles in the base fluid energizes species diffusion physically, which raises the temperature and concentration profiles.
Stefan blowing from the surface to the free stream when the cylinder is shrinking, but the opposite trend is obtained in the stretching case.Figure 2b,c show how the rising Stefan blowing parameter profile for the temperature and concentration are enhanced for both the stretching and the shrinking case.The presence of nanoparticles in the base fluid energizes species diffusion physically, which raises the temperature and concentration profiles.From Figure 3a,b, it can be seen that, with an increasing stretching/shrinking parameter, the velocity rises, but the temperature profile decreases.This is because the stretching/shrinking state improves the HNF velocity, which, in turn, increases the thermal boundary layer thickness and results in high temperature gradients on the surface; the temperature drops as the stretching/shrinking parameter upsurges.The velocity profile declines as the velocity slip parameter increases in the stretching case but, in the shrinking case, the velocity profile is boosted (see Figure 4a).In fact, the velocity is partially transported to the fluid with an upsurge in the slip parameter when the cylinder is stretching, whereas the opposite trend is seen in the shrinking case.The velocity decreases as the magnetic field rises for the stretching case, but it increases for the shrinking case (see Figure 4b).The retarding force increases along with the values of M, which causes the velocity to decrease.
Both the temperature and concentration profile reduce for negative values of the unsteadiness parameter for both the stretching and the shrinking case, as depicted in Figure 5a,b.When the unsteadiness parameter increases, the stretching/shrinking cylinder loses heat physically.As a result, the temperature of the HNF drops.
case, the velocity profile is boosted (see Figure 4a).In fact, the velocity is partially transported to the fluid with an upsurge in the slip parameter when the cylinder is stretching, whereas the opposite trend is seen in the shrinking case.The velocity decreases as the magnetic field rises for the stretching case, but it increases for the shrinking case (see Figure 4b).The retarding force increases along with the values of M, which causes the velocity to decrease.Both the temperature and concentration profile reduce for negative values of the unsteadiness parameter for both the stretching and the shrinking case, as depicted in Figure ported to the fluid with an upsurge in the slip parameter when the cylinder is stretching, whereas the opposite trend is seen in the shrinking case.The velocity decreases as the magnetic field rises for the stretching case, but it increases for the shrinking case (see Figure 4b).The retarding force increases along with the values of M, which causes the velocity to decrease.Both the temperature and concentration profile reduce for negative values of the unsteadiness parameter for both the stretching and the shrinking case, as depicted in Figure Figure 6a shows that, with negative values for the unsteadiness parameter in contradiction to the stretching/shrinking state, the skin friction coefficient decreases.According to the graph, parameter B is dependent on the buoyancy force, and the flow is in an upward direction because of a stronger buoyancy influence that causes the velocity field for the flow of fluid to decrease.Figure 6b demonstrates that Cf rises as a volume fraction of ϕ 1 and ϕ 2 from 0 to 2% against the stretching/shrinking case.The volume fraction of Figure 6a shows that, with negative values for the unsteadiness parameter in contradiction to the stretching/shrinking state, the skin friction coefficient decreases.According to the graph, parameter B is dependent on the buoyancy force, and the flow is in an upward direction because of a stronger buoyancy influence that causes the velocity field for the flow of fluid to decrease.Figure 6b demonstrates that C f rises as a volume fraction of ϕ 1 and ϕ 2 from 0 to 2% against the stretching/shrinking case.The volume fraction of nanoparticles served to significantly improve the skin-friction coefficient for both the stretching and contracting states.Figure 6a shows that, with negative values for the unsteadiness parameter in contradiction to the stretching/shrinking state, the skin friction coefficient decreases.According to the graph, parameter B is dependent on the buoyancy force, and the flow is in an upward direction because of a stronger buoyancy influence that causes the velocity field for the flow of fluid to decrease.Figure 6b demonstrates that Cf rises as a volume fraction of ϕ 1 and ϕ 2 from 0 to 2% against the stretching/shrinking case.The volume fraction of nanoparticles served to significantly improve the skin-friction coefficient for both the stretching and contracting states.

Conclusions
In this numerical study, the unsteady MHD transport phenomenon of HNF flow near a stagnation point over a stretching/shrinking cylinder is explored.Using suitable similarity, the governing PDEs are transformed into ODEs and solved via bvp4c code.The study's main conclusions are listed below:

Conclusions
In this numerical study, the unsteady MHD transport phenomenon of HNF flow near a stagnation point over a stretching/shrinking cylinder is explored.Using suitable similarity, the governing PDEs are transformed into ODEs and solved via bvp4c code.The study's main conclusions are listed below:

•
As the Stefan blowing parameter increases, the profiles of the temperature and concentration rise; • With a rise in the variable viscosity parameter, the velocity profile decreases, but this is more significant for skin friction;

•
The velocity, temperature, and concentration profiles decrease while there are negative values for the unsteadiness parameter; • The addition of the nanoparticle volume fraction leads to a better thermal conductivity performance, resulting in an increase in thermal performance;

•
The skin friction decreases with negative values for the unsteadiness parameter.

.
When c2 = 0, the cylinder is static; when c2 > 0, the cyli stretched; when c2 < 0, the cylinder is shrunk; c1, c2, and c3 are factors, and the charac length is L. The surface temperature, upstream temperature, surface concentratio free stream concentration are symbolized as  ,  ,  , and  .Stefan blowing, s tion point, variable viscosity, unsteady MHD, and partial slip effects are inclu shown in Figure1.When we start the stretch/shrink cylinder, the fluid starts flowin direction but, in our model, we consider only the positive x-direction.The Stefan b effect applied is vertical to the fluid flow direction.Since water is used as the base fl our model, the Schmidt number is fixed.

Figure 1 .
Figure 1.Physical model for the stretching/shrinking cylinder.

Figure 1 .
Figure 1.Physical model for the stretching/shrinking cylinder.

Figure 4 .
Figure 4. Velocity profiles for (a) different A and (b) different M values.

Figure 4 .
Figure 4. Velocity profiles for (a) different A and (b) different M values.

Figure 4 .Figure 5 .
Figure 4. Velocity profiles for (a) different A and (b) different M values.